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On Two-Dimensional Panel Flutter’ 


Y. C. FUNG* 
California Institute of Technology 


SUMMARY 


Theory and experiments of the flutter of a buckled plate are 
discussed. It is shown that an increase in the initial deviation 
from flatness or a static pressure differential across the plate 
raises the critical value of the ‘‘reduced velocity.”’ 

The applicability of the Galerkin method to the linearized 
problem of flutter of an unbuckled plate has been questioned by 
several authors. In this paper the flutter condition was formu- 
lated in the form of an integral equation and solved numerically 
by the method of iteration and the method of matrix approxima 
tions, thus avoiding the constraint of assumed modes. For a 
plate (with finite bending rigidity) the results confirm those 
given by the Galerkin method 

An approximate analysis of the limiting form and amplitude of 
the flutter motion for a buckled plate is presented 


SYMBOLS 


Fourier coefficients of series expansions 

of 2, 2:, respectively 
= amplitude of B,(r), B.(r), Eq. (3.4 

Fourier coefficients of s, see Eq. (1.8 

nondimensional Fourier coefficients of 
z, see Eq. (1.10) 

influence function 

dimensionless influence function, Eq 
(2.3 

normalization factor of Iguchi fune 
tions, Eq. (1.80) 

Eh'/12(1 — v?), bending rigidity of 
plate 

effective elastic modulus, see Eqs 
(1.5) and (1.6) 


Presented at the Structures Session, National Summer Meet 
ing, IAS, Los Angeles, Calif., June 17-20, 1957, under the title, 
“On Panel Flutter.’’ This work is available as IAS Preprint 
No. 749. 

Revised and received October 28, 1957 
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The author is grateful to Prof. E. E. Sechler for suggestions 
in experiments and to Mrs. D. Eaton of C.1.T. and Dr tA 
Brooks and Dr. W. Frank of The Ramo-Wooldridge Corporation 
for obtaining numerical results. 
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Iguchi functions, Eq. (1 
kernel of the integral for aerodynamk 
pressure, Eq. (2.7 
kernel of an integral equation, Eq 
2.11 
plate thickness 
wl /U, reduced frequency of flutter 
nondimensional Fourier coefficients of 
the static pressure p, see Eq. (1.9 
w L/L’, stiffness parameter 
chord length of plate 
Mach Number of flow 
stress resultant in x direction, tension 
per unit span, tension positive 
the value of N, corresponding to z 
buckling load per unit span 
static pressure differential across the 
plate, positive force downward 
aerodynamic pressure induced by the 
deflection z, acting on the top sur- 
face of the plate 
coefficient of p, see Eq. (1.9 
dimensionless static pressure loading 
coefficient, see Eq. (1.15) 
dimensionless aerodynamic pressure, 
Eq. (2.3) 
1/2)p,U? = dynamic pressure of flow 
dynamic pressure ratio, see Eq. (1.14 
A,/As, amplitude ratio, Eq. (3.17 


- the Routh discriminant 


thrust ratio = compressive stress 
buckling stress, in a state corre- 
sponding to 2» 

time 

increment of distance between span 
wise edges 

speed of flow at infinity 

streamwise coordinate 

x/L, dimensionless streamwise co- 
ordinate 

spanwise coordinate 

vertical deflection, positive up 

the initial, unstressed shape of the 
plate midsurface 

the deflection surface when thermal 
expansion or edge displacements 
have taken place, but no static or 
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aerodynamic pressure icts on the 
plate 

Pa ¢/L, dimensionless deflection 

a, A, a4 coeflicients of stability determinant, 
Iq 22 

ay coefficient of linear thermal expansion 

6B,, (m 2 assumed small perturbation of B 

AT temperature rise of plate above that of 
the supports 

1, o phase angles, Eq. (3.4 

¥y QQ} = & 

Gm (aM | kr + Q 

k Wk/(M? — 1 

nN m he as nondimensional Fourier coctlicients of 

, Eq. (1.10 

Am (mn 4 nondimensional Fourier coeflicients of 
sc, Eq. (1.10 

m = p.h/plL = mass ratio of panel to air, 
Ea. (1.12 

y = Poisson’s ratio 

py density of undisturbed air 

p density of plate material 

T dimensionless time, see Ea. (1.18 

w = circular frequency of flutter 

Ww fundamental frequeney of free vibr 
tion, rad per sec 

X membrane thrust — buckling load 


ratio, Eq. (1.17 


INTRODUCTION 


Pex FLUTTER refers to oscillations of a panel (a 
thin plate, a shell, or a membrane) in a flow, for 
which the aerodynamic force induced by the motion of 
the plate plays a decisive part. The oscillations of such 
panels excited by mechanical or aerodynamic noise, or 
buffeting, will not be considered here. 

That panel flutter can occur in an incompressible 
fluid is exemplified by the common experience of the 
flutter of a flag. But as an engineering problem panel 
flutter seems to become important only in supersonic 
flight. In particular, large amplitude panel flutter is 
likely to occur only for panels that are buckled. (The 
transient aerodynamic heating may be the cause of 
buckling.) Because of the usual restrained edge condi 
tions, the maximum amplitude of flutter is limited. 
Therefore, the apparent importance of panel flutter 
lies in the determination of the fatigue life and allowable 
stresses of the panel material, and in its influence on 
the aerodynamic drag. 

The analysis of panel flutter naturally finds its prece 
dent in the works of Kelvin'® and Rayleigh*’ on the 
problem of the generation of weter waves by wind. 
Similar analysis has been made by Miles'’ on the 
flutter of infinite by studying progressive 
waves in a plate which is immersed in a flow. Miles 
shows that in such a plate a straight-crested progressive 
wave interacts with the such a that a 
critical flow speed exists above which the wave ampli 


plates; 


flow in way 
tude will diverge with increasing time and below which 
it damps. Miles!” * 
waves in a thin-walled pressurized cylindrical shell of 
infinite length, under the assumption that the wave 
length is short in comparison with the radius of the 
Again the possibility of flutter is found 


also studies such progressive 


cylinder. 
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It is, however, difficult and unreliable to extrapolatg } 
the results obtained for infinite plates or cylinders { 
those of finite dimensions. In the following we shal 
consider panels with finite chord but infinite span. 

The first recognition of panel flutter was due 
Jordan,'* who suggested that a number of early Germa 
\-2 failures were due to panel flutter. The first wind 
tunnel testing at supersonic speeds was reported by 
Sylvester and Baker.” Analytical studies of the panel 
flutter problem have been published by Isaacs, Hayes 
Miles, Shen, Goland, Luke, Hedgepeth, 
Leonard, Jordan, Fung, Eisley, Nelson, Cunningham 
Ashley, Zartarian, Constant, St. John, and others (se 
References for articles directly related to the present 


Budiansky 


paper). The literature is marked by a high degree oj 
controversy, resulting from the subtle effects of various 
simplifying assumptions used by different authors. 
The present paper is a summary of some of th 
GALCIT. 
experiments will be discussed in order to clarify th 


results obtained at Several exploraton 


physical conditions of flutter in transonic and_ super 


sonic flows. Relevant theoretical considerations wil] 


also be discussed. 


FLUTTER OF CURVED PLATES WITH EDG! 


COMPRESSION IN 


(1) 
\ SUPERSONIC FLow 


The flutter of buckled plates has been a subject oi 
The restrictive among 
and 17 refers to the following 


many investigations. least 
references 4, 5, 6, 9, 11, 
problem. A uniform plate of finite chord but infinite 
span is perfectly flat in the initial (no-load) condition 
It is buckled by edge constraints in the chordwis 
direction. The buckled plate is then exposed to 
supersonic flow on one side and to stagnant air on the 
other side. Under the assumptions that the flow is 
two-dimensional, that the plate is uniform, isotropic 
elastic, and has either clamped or pinned edges, an 
that the acoustic pressure in the air underneath th 
plate is negligible, the theory gives the critical flutter 
speed as a function of the amplitude of the buckle. 

A comparison of the results of such a theory wit! 
experiments described in reference 4+ shows, however 
On examining the exper 
mental the 
theory, it seems apparent that the assumption of per 


a significant discrepancy. 
conditions versus assumptions of — the 
fect flatness of the plate in the initial, no-load conditio! 
is the most offensive. Since the critical flutter condi 
tion of a buckled plate is sensitive to the mean com 
pressive load (the so-called ‘‘membrane’’ stress) in th 
plate, and since the compressive load for any give! 
amplitude of deflection is reduced by any initial curva 
ture, it can be expected that the effect of the initial 
deviation from flatness is to raise the critical flutte: 
speed. 

The effect of initial curvature can be investigated by 
Some details 
A. brie 


outline of the analysis and the numerical results will be 


the same method as used in reference 6. 


of such an analysis are given in reference 7. 


wher 


he 
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presented below. The assumptions listed above art 
rapolate } : : ; A Pid 
ised un, except that the assumption of perfect 
nders t : ‘ me ‘ p 
jlatness in the no-load condition is dropped Phe g@ 
we shal . So : ; 
plate is Shown in Fig. 1. The edges normal to the flow 
yan. ; . : | 5 
ire fixed on a massive structure; the distance between Pia Pcie 7 
tf ' / at | 
due t these edges can be changed by a small amount by Pi si MP Fae 
yermay . z - U>o — y = 
Germa external means, but is assumed to be unaffected by the ee 4 L oe Ae : 
st wind variation of the thrust in the plate. Pct 4 tas Pete t 
x + . . e < yh ® s col 
rted b Let 2;(“) be the initial, unstressed shape of the mid = : Lo 
. -” : : ‘ : —T —>=" - - 
he pan surface of the plate. Let 2)() be the deflection surfac f 
_ Hayes when any thermal expansion or displacement of the , 
x @ @ 
I insky span W1S¢ edges has taken pl ice, but no load acts Fi Notate 
ne } - . 
Ingham normal to the surface of the plate. Let 2(x, ¢) be the 
eTS (See | midsurface of the plate during flutter. Considering a Ik’ hg ik, the Young's modulus 5 


present strip of unit width, and assuming that the plate is thin } ; 
Whereas, if the spanwise strain vanishes 


esree ol | and deflections small, we obtain the equations governing 
Various - and > (see Symbols) KE’ E | ; / ” I: | , 1.6 
Ors. 
of th D\o'(z ei) Ox} V;,(O°S) On 0 7 Eqs. (1.2), (1.3), and (1.4) are to be solved for the 
oratory | DId4(- 2) /Ox4] N.(O°s/Ox?2) + N.. X boundary conditions 
rify the 0°) /Ox-) / p.l(O-z Of) b Be =((). t) o(L. t 0) } a 
| super i ( is ; 4a) 
I ‘ ‘ ; ) f \ Of 
Peni Che aerodynamic load will be assumed to be of the O"2/ On = o3/0 be 
orm for pinned edges and 
py l?/V M 1); [Os(x, 2) Ox] 4 (0, t) = 2(L, t 0) { 17h 
' ) 
\/ 2)/(M 1) }(1/ U) [Os(x, t)/ot| 4 Oz Ox) (0, 7) 02 /Ox)(L, f Ot 
GI : 
11/(M 1) ](1/U’)[d2(0 Of |; 3 ¢ 
for clamped edges 
. rhis is valid for slow oscillations** (downwash velocity 
ee speed of flow), if the Mach Number is sufficiently (1) Approximate Solution by Galerkin’s Method for 
among : ; : ; ; Panel With Pinned Edges 
tlow; remote from 4/2; and also for higher Mach Numbers 
LLOW1ng . . . 
“afin; when the so-called piston theory! becomes applicable Let 
Infinit . . é , 
fit | (in which case (1/ 1) and (\/ 2 \/ 1) z 
1d1tion . ; 4 
| should be replaced by ./ and 1, respectively). \" > _ MT an (mer/] (1.8) 
rdwist bh ; aa 7 2, \a .sin (max L : 
i Phe compressive load — NV, is computed as follows }. | 
£0 . . é ~ 
If w. is the increment of the distance between the two 
on the ' : \ 
. | spanwise edges, and A/7/ is the temperature rise of the P on / I 1.9) 
fowl - p Po 2, Ry sin (max, L J! 
plate above that of the supports, then 
tropic 
Ss and i which satisfy the boundary conditions (1.7a), and 
vie \ E’h/2L) [(Oz; Ox) Oz Ox)*| & 
ith th . V D/L 
flutter dx Ee'h/L)us + E"ha;Al 14 ; ; 
™ which is the Euler buckling load of the plate. 
vid where, if the spanwise stress vanishes, Let the following dimensionless parameters be 
y witl 
ses introduced 
wever 
exper \ J ° 
- a 1" ‘2 
= th \ V 31 yds / ie — ty mth component of » Zo, 1.10) 
Ol per IB \ H }; \ M ) 
1ditior 
condi V compressive load in plate corresponding to 2 
; ; : i ; Bi} 
1 com NN.) Euler buckling load 
in the 
givel u mass ratio ph! pol 1.12) 
curva 
initial T dimensionless time Ut/L 1.13) 
flutter 
¢ dynamic pressure ratio pol/?L*/DaxtV M 1.14 
ted by . 
“1 PF; V 3(1 y)E'/E L4/a4Dh Be 
letails 
briel | 


lhe governing equations are reduced by Galerkin’s method into 
will be 
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\ 2 
B, | > n>(B,,”? — X,?) + m? — ‘| = Mn, — Pe km — } C 
21 m: ont 
QO jX, 2mn[(1 — (-1)"*"] M? — 2 dB,, dB, | Me 
4 B, + — + uv I? — 1 (m ae , NV) By 
m? (az, (m+ n)(m — n) MP -— 1 dr dr? § i 
The symbol ’ indicates that the term m = n must be excluded. | . ' 
Numerical results have been obtained for the case of a two-term approximation, N 2. In this case we notic : 
that the factor wt 
xX =A? — BY? + 4(\.? — B,?) [(—.N,) (—N,,,)] -— s 1(—N,z) — (—N,,)]/(— Neer) 1.17 | cur 
appears in all equations. Hence we write, for m = 1, 2, | a 
By(—x +1 —s) = Aw — Poki + Q) (8/3) Be — [CM — 2) CMe — 1)])(dB, dr) — up AP? — 1(d2B, dr*)} 1.18 I 
easy 
Bo(—x + 4 — s) = 4dov — (Poko 4) + (Q/4)} —(8/3)Bi - pre 
(MW? — 2) (MP — 1))(dB2/dr) — uw V MP — 1(d?B2/dr?)} (1.19) | of ( 
At static equilibrium all derivatives with respect to 7 vanish. One may then eliminate B,, B, to obtain an equatior ~ 
relating x and Q. In case Ax» = ky = 0, the resulting equation is quadratic in Q? and yields simply ‘ “ 
1602 Que — =~ | 1244 — s — x)? — x) } forn 
+ ie -s— wo +s ~ ae : I “> * | 1.20 
9 2(a° — x) L (Au — Pok1) 
An example here will clarify the nature of the solu- 
tion. In Fig. 2 is shown the relationship between the 5 ] T | ~—s . 
membrane compression ratio x and the dynamic pres- | dio?! 
sure ratio Q, for the special case of Aip = 1, Ax = Po = 0, | |B +0 
the ordinate being Q*. The solution depends on the 4 | bes rhe 
initial compression ratio s. When s is small, one of the | | I> 
roots Q* is negative and the other is positive. Both ™ | ; 
roots tend to infinity asymptotically as x approaches Oley 
Aw?/(1 — s)*, at which the plate becomes flat When s ’ 0 Ojlon / | | inde 
is a sufficiently large number, say s = 0.5, two curves Q? Vom | 
are shown in the Figure; the lower one can be shown 02 A | | 1% 
by the method below to be unstable, and the upper one 2 9.6 ] Re 
has a peak. No solution exists for Q greater than this 07. | enn t 
relative maximum. It can be shown that the equi- 19] | | - 
librium configurations corresponding to the front por- ' oo“) ‘ | | the 
tion of the curve, OA, are stable, whereas those corre- Gis the 1 
sponding to AB are unstable. // XH 7?8 | mpl 
The stability of the static equilibrium configurations ae ot Fae * a 5 the e 
will be defined against infinitesimal perturbations. OT77 77 I a > ~*~ ; = 2 differ 
Let B,, Bs, x be the values at a static equilibrium con- is é sith \ | | OAB 
figuration. Consider the variational equations by / My) a ] \ | | OAB 
imposing small perturbations 6B,(r), 6B2(7) onto By, Be, a a \ eal J also 
respectively. Keeping only linear terms in 6B,, 6Bz, Fic. 2. The square of the dynamic pressure-bending rigidity | boun 
_ he tie nent aeventiel equations for 6By, Be. Se sev geht ely dite eusian confi 
On assuming solutions of the form infins 
curyvs 
6B, ~ ee’, 6B, ~ e”’ (1.21) a = (1/Q)((AP — 2)/(1P — 1)] Xx ait 
(17 — 5s — 5x + 2B)? + 32B. ADE 
two linear homogeneous equations are obtained. For aw = (4/Q*)} (1602/9) + (lL —x —s\4—x-— 5) -]} ape 
nontrivial solutions the determinant of these equations 2(1 — x — s)(x — A? — 42") + GB LR 
must vanish, which may be written as 
i The condition of stability of the disturbed motion 4T¢ \ 
ase* + aze® + ave? + me +a =O0 (1.22) is that all the roots of Eq. (1.22) must be pseud ation 
negative (no root has a positive real part). This re-| Pec! 
where quires that a, a), ..., a; and the Routh discriminant P. 
a, = (M? — 1)p? ®, have the same sign: ) igait 
initia 
a; = 2[(( 17? — 2) V i? — 1 |u R = AyAva; — apa;” — asay" reduc 
ao = [(M? — 2)/(A? — 1) + (MM? — 1)p?ta? + 2[(I? — 2)3/(M? — 1)7/7] x} vent; 


((A7? — 1)*#?/(MP — 2) ]uay (aru) — 4[(M? — 2)/( MP — 1)Paof (1.23 Fis 
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Checking of the conditions of stability for all static 


equilibrium states can be done numerically. As 
» J? — 2), one of the conditions of stability is 
> V2. When this 


generally found that for smaller values of s and Aj the 


condition is satisfied, it is 


stabilitv is governed by the Routh discriminant, which 


becomes negative when Q exceeds a certain critical 


value. For larger values of s and Ay, when the Q vs. x 
curve has a relative maximum, the stability is governed 
by a, Which becomes negative when x exceeds the 
critical value at which Q reaches its maximum. 

The critical condition will be expressed by Q It is 
easy to show that 7/Q is proportional to the reduced 
velocity lL’ (wf), and that the inverse of the cubic root 
af O, V 


required to prevent flutter. 


, is proportional to the thickness of plate 
An extreme case is of interest. Consider a slightly 
warped plate which assumes in the unstressed state the 
form of a low sinusoidal arch with a center rise 

(1.24) 


Lay fh V30 — PEE] = O.58AA, 


We assume 


Aw| < 1, Po| <1, Aw = Re 0 (1.25) 


[hen it can be shown that, if s is not too close to 1, 


> V2, and wu is sufficiently large (say, > 20), 


0.. = (3/4), [U7 — 5s) 4P — (4-5) —s)f'? (1.26) 
independent of 7, uw, Py and Aw. Q., = 45°16 when 

0; it decreases as s increases until the expression 
1.26) becomes invalid when s ~ 0.7. 

Results of calculations based on the equations above 
can be summarized graphically as follows. In Figs. 3a 
and 3b is plotted the dynamic pressure ratio Q versus 
the arch-rise ratio \; for plates that are perfectly flat in 
the unstressed state. Remember that A; describes the 
amplitude of the buckle due solely to the movement of 
Fig. 3a gives the dividing line between two 
above the curve 


the edges. 
different states.of static equilibrium: 
OAB the plate will be blown flat by the wind, below 
OAB the plate may remain buckled, but of course may 
also become flat. In Fig. 3b are shown the stability 
boundaries of the initially flat specimens. A buckled 
configuration below the curve OAB 1s stable against 

A flat configuration in the 
OCDEAO is stable 
In the small area 


infinitesimal perturbations. 
curved-triangular region also 
against infinitesimal perturbations. 
ADE both flat and buckled configurations are stable. 
Above the curve CDB no static equilibrium is stable. 
Fig. 3e shows the stability boundaries for plates that 
are warped in the unstressed state. The initial devi- 
flatness is assumed to be sinusoidal and 


For this Figure it is assumed that 


ation from 
specified by Arp. 
F \ 0. The region above each curve is unstable 
igainst infinitesimal perturbations. It is seen that the 
initial deviation from flatness, through its effects on 
reducing the compressive stress, is beneficial to the pre- 
vention of panel flutter. 

! gives some calculated values of Q, 


Fig for several 


DIMENSIONAL 
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FIG.3b- STABILITY BOUNDARIES, INITIALLY 
FLAT PLATE 
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FIG.3c- STABILITY BOUNDARIES, INITIALLY 
WARPED PLATE 
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values of the parameter P)k;, with k» 0 and Ajo 0.5 
(except that the curve for ?,k; = 0 is drawn for Ayo = 0). 
The effectiveness of the static pressure in raising the 
flutter boundary is seen. (Note, however, from Eqs. 
(1.25) and (1.26), 0Q,,/ oP, 0 when P, 0.) 

In reference 25, Sylvester and Baker showed experi 
mentally that pressurization is an effective method of 
raising the reduced velocity of flutter. Fig. 4 is in 
agreement with this trend, even though it refers to a 


sinusoidally distributed load, 
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> 4 Exp. points for which a negative value of Pok; indicates an excess oi) hc 
*;, \4 | “Lor mplitude”fiutter = ‘ 
ile \ Ts} ; 4c aeaeemne pressure underneath the plate (pressurized chamber) | tio 
°o|* = Ay + . + — 4 la oge ° ° ~* ‘ 
el 4 AF 1Z 1% NI nee FeOhy rhe stability boundaries in Figs. 3 and 4 are pra tra 
° owe hk PP tically independent of the mass ratio uw for practic 
g a - ‘ad I Cal} ats 
© tk Oke a” 2 sot T T | values of u (say, 20 to 150). 
= a\e \ 4 .. Pree al RS hme) | - sho 
2 4 4 ~*~ el o ac, le © | 
° \ ~ — . a] + 
eo A le. [-—};—"P 1.2) Plates With Clamped Edges , 
2 jt ee ee SE seen —— : 
3 i ie — . , 
5 —— ——— ae For plates with clamped edges we use Iguchi func | edg 
a = | | | i! | | ee tions to describe the deflection surface and the externa 
: loading: 
é | | | pom 
re) 2 4 6 ~ 10 12 14 16 \*| » \ on | 
A, =V3-(initiot arch neight)/(plate thickness) 20 gy > i. Cart 1.28 | 
Fic. 5. Comparison of theory and experiment for plates with I. | m if; 7 \ | 
clamped edges 
p PoLRmCmFm(X ) 
where 
Foa(X) iY — 1)? + (-1)"X*7X¥ — 1) - 
(1 ma) sin marX 1.29) } 
i ce 
and c,, are normalizing factors so chosen that 
5 
MAX Cm F(X) | 1.31 
Os X <1 
\ The Euler buckling load corresponding to /\Y) is 
U"xt"FLUTTER PANEL | 
OMITTED —N,. 11.0651 D/L? is 








which must be used in the definition of s. Using the | 








Galerkin method with a two-term approximation, and | 
on setting Ek’ = E (1 — v*), v? = 0.1, the static equilib 





rium configurations and stability boundaries for ), ( 
0:5. 10, 3:0, 5.0: Xen = Po 0, have been computed 11 
reference 7. Some results are plotted in Fig. 5. | 


1.3) Comparison With Experiments 


In an experiment described by Eisley, a series of 1] 

by 1l-in. aluminum panels, ranging in thickness fron ; 
0.025 to 0.051 in., were mounted on the top surface of a | 
. wedge-shaped model and tested at Mach Number 2.15 ' 
Fic. 6a (top). Assembly of model and buckling amplitude con- — JPL oe a plate could tee bectill 


trol mechanism. FIG. 6b (bottom). Installation of the model in 
the 12-in. tunnel. from outside the tunnel by moving the plate trailing 





edge with respect to its leading edge. The onset ol 
co flutter was observed. The history of the motion was}  [-f 
‘ve aR ee ee recorded by strain gages fastened to the lower surface | 
of the plate. | 

The model is illustrated in Fig. 6. The panel was} {J 
| | located near the leading edge to minimize the effects ol! 
ik i toen ae ef 008 boundary layer. The leading and trailing edges of the 

4 008 | et > a panel were fastened with screws. The side edges were ad 
Tome, se free. The lower surface of the panel was vented to the 
oi free stream static pressure by slits along the free edges 
——— > ————-$ 0 008 A combination screw, wedge, and slide arrangement was Ll 








sae = a a | |_|} Pee) Pee used to vary the amplitude of the buckle by moving the 
: trailing edge. 
z —— : Two distinct types of oscillations were observed 





[| | “= ‘i 5 ppp ape} 0.009 One occurred at small amplitudes of the buckle. It was | ) 
yada , iter Tis mostly irregular and of small amplitude. Hence it wil 


= owe be referred to as “low amplitude’ type. The other 


Fic. 7 (Le ft Panel deflection history of a “‘large amplitude” the “large amplitude’ type, was nearly harmonic time 
flutter. Fic. 8 (right). Panel deflection history of a ‘low ampli- : : : ee ‘ se 
tude” flutter wise. An examination of the oscillograph records 
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TWO-DIMENSIONA 
showed that, in the large amplitude flutter, the deflec- 
tion surface changed in such a manner that it resembled 
traveling waves, whereas in the ‘‘low amplitude’”’ oscil- 
lations, it appeared to be standing waves. These are 
shown by the deflection histories in Figs. 7 and 8. 

It was most interesting to watch the behavior of the 
plate in the wind tunnel as one gradually forced the 


i 


PANEL FLUTTER 
nel pressure was set at some selected value with the 
The panel was then 
At 


some point the low amplitude oscillation was observed, 


panel in the unbuckled condition. 
buckled and the amplitude was slowly increased. 


mostly irregular but sometimes also nearly harmonic 
but generally it could be stopped with further increase 


in the amplitude of buckle. With still further increase 






















































































































































































































































































































































































edges of the panel together. In a typical test the tun of the buckle a point was reached at which large 
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amplitude flutter occurred. However, at no time was 
it possible (within mechanical limitations) to stop the 
large amplitude flutter by further increase in the 
buckle. In some cases, particularly in those which 
required very large buckling amplitudes to initiate 
flutter, no low amplitude oscillation was observed. 
Such a sequence of events is illustrated in Fig. 9 by the 
oscillograph records from strain gages attached to the 
plate. 

The test conditions described above were such that 
the assumptions of the theory presented herein are 
approximated. In Fig. 5 both the theo- 
retical and the experimental flutter boundaries. (Note 
that Aj 1 corresponds to a sinusoidal deviation from 
flatness with a 
equal to 0.58 times the thickness of the plate.) 


are shown 


maximum deviation approximately 
Un- 
fortunately, the actual warping of the test specimens 
was unknown. It is seen, however, that each experi- 
mental point could be explained by a certain amount of 
initial warping.* 
To critically compare experiments and theory, one 
would have to examine the effects of the simplification 
introduced by the two-term approximation, of the 
possible deviation of actual boundary conditions from 
being clamped, of the influence of the boundary layer, 
and of the acoustic pressure underneath the plate. 
Thus many questions remain to be investigated, but 
the principal features of the problem are now clear. 
Panel flutter of a buckled plate at a sufficiently high 
Mach Number may be roughly characterized by a 
constant critical value of Q, if one were willing to regard 
the effect of initial warping, buckling amplitude, etc., 
as secondary influences affecting the scatter of this 
constant. 
(2) LINEARIZED EQUATIONS FOR A FLAT PLATE AND 
THE GALERKIN’S METHOD OF SOLUTION 


The Galerkin method, or the Rayleigh-Ritz method, 
is widely used in aeroelasticity, but no mathematical 
foundation has been established for its application to 
non-self-adjoint eigenvalue problems such as flutter. 
Recently, doubts were cast over its applicability to the 
panel flutter problem. Miules,'* and Goland, Luke, 
and Constant,’ using the assumption of small reduced 
frequency, (k < 1), have shown that the Galerkin 
method yields a spurious flutter boundary (one that 
for a membrane’ (without bending 
Ashley and Zartarian! 


under the 


should not exist) 
stiffness) in a supersonic flow. 
have shown the same (also for a membrane), 
assumption that the Mach Number is sufficiently high 
that the ‘“‘piston theory’’ applies. Therefore, it be- 
comes a matter of great interest to see how far the 
Galerkin method fails. One way to do this is to com- 
pare its results with those obtained by a method of 


proved validity. The method of Laplace transforma- 


* Except for two “low amplitude” points. The randomness 


of the “‘low amplitude’”’ oscillations is not explained, nor is its 
stoppage. There remain some doubts as to the possibility of 
buffeting of the model in the tunnel 


implitude flutter is needed 


Further investigation of 


the low 
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tion, as outlined by Goland and Luke,* can be used fo, 
this purpose. In the following, however, the iteratio; 
method proved by Wielandt*® and the method of matr; 
approximation of the integral equation will be used. 


2.1) Statement of the Problem, Fundamental Equation; 
Consider an infinite strip of flat plate, one side ; 
which is exposed to a uniform supersonic flow in th 
direction of the chord of the strip, whereas the othe; 
When not deflected 
Flat solid walls parallel 


side is exposed to stagnant air. 
the plate lies in the x, y plane. 
to the uniform flow extend to infinity both in front an 
in rear of the panel. (See Fig. 1.) The problem is t 
find whether the plate may flutter. 

If flutter exists, the plate will oscillate harmonicall 
without any external excitation other than the inter 
action of the flow and the plate. It is assumed that 
the amplitude of oscillation is so small that the elasti 
and aerodynamic actions can be described by linearize 
equations. The acoustic pressure beneath the plat 
will be ‘neglected. The latter has been discussed by 
Nelson and Cunningham*! and is expected to be un. 
important unless the flutter frequency happens to be 
near the resonance frequency of the chamber beneatl 


the plate. 


curvature of the plate, or a static pressure acting on 


the plate, has no effect on the flutter boundary. 
Consider a strip of unit span and of chord length L 


In the linear range, the elastic behavior of the panel is 


specified by the influence function which describes the 


deflection under a unit static lateral load 


acting at an arbitrary point. 
two-dimensional, the coordinate y plays no role in the 
Let c(x, &) be the influence function describ 


pattern 


problem. 
ing the deflection z at x due to a unit load at €. 
load distribution f(é, ¢) yields a deflection surface 


z(x, t) = 


. 
| c(x, €)f(E, Ode (2.1 
de 


where S is the domain over which the lateral load acts 
The coordinate system is shown in Fig. 1. For the 


flutter problem it is assumed that the motion is suff 


ciently slow in comparison with the time of travel oi 
elastic waves in the plate across the chord so that 
c(x, &) is uniquely defined and is independent of the 
loading rate. In this formulation, the inertia fore 
must be regarded as an exterior force. Hence, 

f(é, t) = —p.h(0*s/Ot?) + pal€, t) 2.2 


where p, is the density of the plate material, / is the 
plate thickness, and p, is the aerodynamic load. For 
(x)e", the 


To reduce Eq. (2.1 


harmonic oscillation, 2(x, 1) aero 


dynamic load Pp, is also harmonic. 


into dimensionless form, let po be the density of air at 


infinity and w, be the fundamental frequency of free 
vibration of the panel in a vacuum, and define 





Since linearity is assumed, a small initial 


Since the condition is 


Then a 
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of matriy 
Pp used. 


[quation; 


e side oj 
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A@=</L, = t/L, Z 2/L 
C(X, =) ph Lay?c(x, &) | 
Palf, #) (poU?/2)P(=)e™ (2.3) 
lu pul. pi = density ratio of air to plate] 
k wh U, k wo l 
7) ad | 
Then we have tt) (Rk? ky?) | C(X, =)Z(Z)d= + (1 2wko?) | C(X, =)P(=)d= (2.4) 
The aerodynamic load is induced by the downwash distribution “ 
UW(=) = U[(O OZ) + tkJZ(=) (2.5) 
” ~ ai om d os 
P(= = +W(O)G(E) + | G(= — v) + 1k} W(e)ds 
V i? —1 nF dv f : 
a _ 2.6) 
< j °= O a | 
= W(=) + W(v) + ik} G(= — v)di 
Vire-—i1t i = f 
where G(x) e~ «Ms 7 (2x), K Mk(M? — 1) 2.7) 
|, is the Bessel function of the first kind. Note that when k 0, 
P(=) -—-(2V ‘A 1) W(E) = 2U V M? — 1)(dZ dz) 
as is expected. 
Consider now the integral in the last term of Eq. (2.4); writing the dummy variable as &, 
. *! *) 73 fa) 
f= — VM . | C(X, &)P(é)dé CCX, &) W(g) ~ wie)( + it) G(é - v)dv' dé 2.8) 
2 J J J df ! 


An interchange of the order of integration of the double integrals gives 


*! *! 7) O 
I C(X, e)W(g)dé + W(0) a c(x, ¢)( < + ieee — vats at 
J, of 


Integrating by parts, we obtain 


*! , : ? ; *! - *) ae O : at 
I Wi(v)C(X, 1)G0. — v)dz 4+ Wa) - G(é — v)| — + ik | CX, E)die di 2.9) 
J VJ, dé 


Eqs. (2.8) and (2.9) state in fact the well-known reverse-flow theorem: that the integral over an airfoil of the 
product of the pressure in the forward flow times the downwash in the reversed flow is the same as that of the product 
of the pressure in the reversed flow times the downwash in the forward flow. This is recognized if one regards 
C(X, £)e” as the downwash in the reversed flow, (X being a parameter). 

Returning to Eq. (2.9), substituting Eq. (2.5) for HV’, and then back to Eq. (2.4), we have, on integrating by parts 


several times and reducing, 


ie ~J int — “1 3) ee 
r(x) : Z(v) + C(X, 1)g(1 — v) + g(fé —v)| — + tk | C(X, é)dé + 
uke? M? — 1. | * de 
O , ae Fe rae eae 
~ + ik — pk?V M? — 1) CX, v)¢ dv — C(X, 1)Z(1) — CCX, 1)G(1)Z(0) 4 
ov y ukyeV M2 — 1 
ie 0 . 
Z(0) aa - ir | CX, gael 2.10) 
; of J 
where g(x) = [(d dx) + tk]G(x) { 
KC M51 Js(«x) T (2 M)Jo(xx)]$ (2.11) 


For plates with pinned or clamped edges, 


Z(0) Zi) = CiA. 1) = 8 


Hence, 


; ‘1 a) j 
Z2(X) = — (0) {| + ik — VM? — 1 Ke | C(X, v) 4 
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k 3 325.35 375 4 425 455556769 | 
-0.8 Symbol © OB A OVYVBADGH ADH F 
-1.0 | | ! ek err a —a eer 1 eee 
0.03 005 0.1 0.5 1.O 5.0 10.0 50.0 100.0 
Real part of eigen value 
Fic. 10. Complex eigenvalues of the flutter equation for M = 1.3, 1/u = 0.03. Each symbol represents a specific value of & 
loci of the smaller eiv’s become complicated at higher values of & and are not shown here Real valued eiv gives the flutter boundary 
ke = wo L/U = 0.0913366- +7 RI eit 


Imaginary part of eigen value 
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Flutter Modes at Low Supersonic Mach Numbers 
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Solution of the Eigenvalue Problem and Comparison 


| (2.2 
With the Galerkin Method 


It is most convenient to regard Eq. (2.12) as specify- 


| ing an eigenvalue problem for the functions ZY) and 
| : ; , 
| the eigenvalue ky*. & is real-valued since harmonic 


oscillation is assumed. For each set of specific real- 


valued k, uw, and .\/, the kernel of the integral equation 
is complex-valued. 

and the function ZY) are complex. 
ry the parameters &, uw, and .l/ systematically in 


Val 


Hence in general the eigenvalue 


[t is necessary 


a complex- 
valued k The lack of 


in existence theorem for such a real-valued eigenvalue 


der to locate a real-valued eigenvalue k 


has no physical significance 


is the basic difficulty of our problem. 

already in the most natural form for the 
iteration method. An 
We substitute Z 
right-hand side, evaluate the integral, 
ind collocate with 27 
ipproximate eigenvalue ko*. The right-hand side may 


Eq. (2.12) 1s 


arbitrary 
the 


ipplicati m of the 


function Z) is assumed. ee 
integral on the 


at some point to obtain an 
then be designated as Z;, and the process repeated until 
Wielandt*' 


has an eigenvalue, that eigenvalue can 


it converges. has proved that if the integral 
equation (2.12 
be obtained by such an iteration procedure, and that 
the iterated eigenvalue ky)? always converges to the one 


with the largest modulus. Biorthogonality relation- 
ships may be used in the iteration procedure to derive 
successively smaller eigenvalues, and modified formulas 
can be devised in the case of two or more eigenvalues 
having the same modulus. 

Numerical solution has been obtained with the help 
of Dr. J. A. Brooks and Dr. W. Frank at the Digital 


Computing Center of The Ramo-Wooldridge Corpora 





tion, on the Univae 1103 machine, for a plate (with 
inite bending stiffness) clamped at both the leading and 
ng edges. For such a clamped plate, 
Ww = [(1.d519)?/L27/V D/ph 
the dimensionless influence function is 
C(X, =) = [(1.51r)! 6]X2%(1 — =)? X 
[B(1 —- X)E-— N(1 — =)] for(¥ <= 
= [(1.51r)*/6JE*(1 — X)? X 
[301 — E|)X —-s(1 — X for (X >t 


rhe function Z(.Y) was represented first by 10 and then 
by 20 equal divisions across the chord length. It was 


found that convergence was readily obtained, but the 








uN 


possibility of having more than one eigenvalue with the 
same or nearly the same modulus was troublesome. 
It was found, however, that taking 10 divisions across 
the chord was sufficient for the problem, for the eigen- 
values determined by taking 20 divisions show little 
difference. 

As the program progressed, new machine methods of 
computation for non-Hermitian matrices became avail- 
able.2 The integral equation (2.12) was therefore 
approximated by a matrix equation, again with Z(Y 
represented at 10 equal intervals, and the ninth-order 
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matrix equation was solved. The results agree with 
those obtained by the iteration method. 

Some details about the variation of the eigenvalues 
with the parameters k, .1/ In Fig. 10 
is shown the variation of all 9 eigenvalues with the 
It is 


u are interesting 
reduced frequency &, when .1/ 1.3, 1/4 = 0.03. 


seen that when & is small, the eigenvalues occur as 
nearly conjugate pairs, which accounts for a difficulty 
with the iteration method. 

When 


largest real eigenvalue k 


the 


> can be obtained by interpola 


sufficient number of cases were solved 


tion. The comparison with Nelson and Cunningham's 
results*! obtained by the Galerkin method is shown in 
Table 1. 


In this case the Galerkin method appears satisfac 
tory. 

There are other solutions confirming the results of 
the Galerkin analysis for plates. Professor Shen has 
informed the author that a direct solution of the differ- 
ential equation for a panel with bending stiffness con- 
firms the Galerkin analysis. William 
Rodden of North UCLA, 


working under Professor Miles, has also solved the same 


Furthermore 


American Aviation and 
problem independently by means of the ‘‘aerodynamic 
influence coefficients,’ again confirming the Galerkin 
analysis by a method in which no mode functions 1s 
The special 

Mach 


assumed. case of static aerodynamic 


force (high Number) has been confirmed by 


Hedgepeth. 


2.3) Results for Clamped Plates at Low Supersonic Mach 
Numbers 


The iteration method vields eigenfunctions Z(x 


which resemble the observed flutter modes* at low 
supersonic Mach Numbers. In Table 2, several 
solutions on the flutter boundary are listed. The 
eigenfunctions are listed for Y = 0.1, 0.2, 0.9; 
Z being zero at XY = 0, and 1. The physical displace- 
ment is 
s(x, t) = RI LZ(x/L)e™ = 

LZpr(x L) cos (RU L)t — LZ,(x L) sin (RU Lt 


A few typical flutter modes are shown in Fig. 11. It is 
seen that the phase shift across the chord is very small. 

It is of practical interest to summarize the theoretical 
results of the linearized problem in the form of mini- 


mum panel thickness required to prevent flutter. 
SinceT 

‘ss 7 1.51 m)- p E 

ut Va — ») WN, g 


each panel niaterial and flight condition specifies a 
hyperbola which intersects the flutter boundaries on 


* The experiments will be discussed in another paper 
t It is of interest to note that for a clamped plate 
S(V E/ps/U)(h 


ko = woL/U = 6 


which shows clearly the dependence on the velocity of sound in 


the plate 
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I T T T T T Let us consider the case of plates with pinned edgy 
Edge fixity : s : , whe 
—— pinned and again use the two-term approximation together 
, es clomped 1 ‘ ° om ° ;‘ " . 
\ = Present-oncifsis - clamped (g=0.05 with the simplified aerodynamic expression, Eq. (1.3 
4 waa - : 3 : a 
‘ \ J rhe governing equations are therefore |see Eqs. lj 
, —--— -—_ 7-1 (1.19)]: 
> 3 Buckied sd )] 
¢ 6! a R= 00, Ref 3 j Sin 
re M4 J 4 ) ) ) ) 
: a a ow uv JM? — 1(d?B, dr?) +[CW? — 2)/(W?- 1)]x 
$ $=00 ity : | { 
$ 0004 _ ta, ria, (dB./dv) + (B/OK— x + 1 — 2) = 
- (8 3)B. — (1 Q)(Aw — Pok) 0 whe 
: 3. 
uV AM? — 1(d?B./dr*?) +[(.W? — 2) /(W? -— 1)]X 
L ! l ! ! ] \ | (dB. dr) + (4B. Q)(— x +4 — 5s) + 
3 4 5 6 ? 8 2 
Moen number, M ; (8 3)B, — (4/Q) [40 — (Po 4k] = 0 
Fic. 12. Theoretical minimum panel thickness required to pre : Phe 
vent flutter as a function of Mach Number for aluminum panels in x =A — BY? + 4A: B 5 ee er 
sea-level air; infinite span; two-dimensional flow, except for ‘ 
finite aspect ratio panels which are based on aerodynamic “strip” with the motion described by He 
theory from reference 27 : 
zs = [h/V30 — v2)E’ E] X | 
the 1 u, ko plane. This intersection determines the [Bi(7) sin (wx L) + B2(r) sin (24x/L)]{ 
panel thickness ratio. Fig. 12 shows such a summary th W3(1 2)! /E] X 9 
ri F ; ; : Zo = 1 al — »)'/L 1 
for aluminum panels in sea-level air. Much extrapola- . sj I ee, ey | io 
: . ; : - A; sin (mwxv/L) + Ao sin (27x L) 
tion of available data is needed to sketch such a Figure; ie ; 
hence, the numerical values are rather crude. The Even these drastically simplified equations are to 
region around |/ = 7/2 requires further investigation. complicated for a complete solution. We shall at : 
All curves refer to undamped panels with no membrane tempt only a first approximation under the assumptior \ 
stress (.V, = 0), except the one labeled “buckled,” that flutter motion exists and is essentially a sinusoidal 
which refers to \,, = 0, and ; > 1.5; and another function of time (as suggested by experiments). The \ 
labeled g = 0.05, where g denotes the usual structural following analysis, first outlined in reference 6, follows 
; ee 5 i : bia wii 
damping coeflicient. the well-known Kryloff and Bogoliuboff!* procedure 
Some numerical results obtained recently will be dis Th 
(3) LIMITING AMPLITUDES OF FLUTTER cussed below. is t 
= , , ' Let us assume that 
rhe previous studies are concerned with the deter af 
mination of flutter conditions under an arbitrary B, = A, sin (kr + ¢)) = A1 Sin 4), | cas 
infinitesimal disturbance. For a nonlinear system, $=kr+o ae 
it is possible that a structure stable with respect to B. = Ao sin (kr + gd) = Asp Sin 6. | a 
infinitesimal disturbances may still flutter under some 6. = kr + 
finite disturbances. Furthermore, experiments show 
; : ; where A,(7), Ax(7), O(7), Go(7) are slowly varying 
that the amplitude of panel flutter does not grow in- ; , Lae = es mas 
= aisehd : functions of time. The second-order differential equa 
definitely, but rather a limiting cycle is developed. - ; as 
a ener “ae tions (3.1) can be transformed into a set of first-order 
Phe possible occurrence of limiting cycles is suggested lif ‘al ; : 14 Diff Rages wh 
. . , c ae differential equations tor A), Ae, di, de. jilierentiation 
by Eqs. (1.2) and (1.4) in which the term .V,0*: Ox See ‘a ie se iti 
; , ; é ; of Eqs. (3.4) gives 
represents the elastic restoring force contributed by the 
membrane tension .V,. If .V,; is positive it tends to (dB, dr) = A,k cos 4, 4 In: 
flatten the plate, if .V, is negative it tends to buckle the (dA, dr) sin 0; + A,(d¢d, dr) cos 8 of « 
plate. From Eq. (1.4) it is clear that under suitable ; ; a \ 
Si ; ’ ' If we impose the condition 
conditions .V, becomes negative when /z! is small, but resi 
N, always becomes positive when 2 is sufficiently (dA, dr) sin 0; + (dd, dr)A, cos 8; 0 3.0 tion 
large. Thus, the influence of the nonlinear term - 

, z : ae, aaa ; then dB,/dr = Aj,k cos A; 5.0 
becomes increasingly strong in limiting the amplitude 3.1 
of flutter. d*B, dr? = —Aj,k* sin 6, + (dA, dr)k X | 

cos 6, — (dg, dr)A,k sin @, (3.7) i). ' 
| co 
The first of Eqs. (3.1) then becomes Fig 
Nu 
/ at 
uv M? — 1 [k(dA, dr) cos 6; — (dg,/dr)Aik sin 6; — A,k? sin 6;] + [(.17? — 2) (M2 1) ]A,k cos #; 4 _ 
(A1/Q) sin 0(—x + 1 — s) — (8/3)Aoe sin 6 — (1/Q)(Ay P ok) 0 (3.8 
' 
From Eqs. (3.5) and (3.8), we obtain _ 
fo= oth 
uv AM — | k(dA, dr) = cos 6; hy A, Ao, 91, P2, kr)! ( | 
(3.9) | ca 


uV M? — 1 kAj(do,/dr) = — sin & fy oo \ ass 
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TWO-DIMENSIONAL 
where 
V M? 1 uk*A, sin 6; — Ww- — 2 VW 1)] 
Similarly, uv MP? — 1 k(dAs dr 
uv MP — 1 kA2(dgd». dr) 
where 
V M2? — 1 wk?A,s sin 6 — MM? 2 \f* l 
[he right-hand sides of Eqs. (3.9) and (3.11 


period, we obtain the average values dA, dr, ete. 


admit, with respect to 6, a period 27. 
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Aik cos 4 (A; Q) sin 0) (—x + | s 
(S'3)Ay sin 2 + (1 Q)(, Pok 3.10) 
cos 6 ) / 
3.11) 
sin 42 foi \ 
|Aok cos 6 i4 x | QA», sin 6 
S 3)A, sin @ + (4 Q)|4A Po 4)k: 3.12) 


If they were averaged over a 


For steady oscillations, these average values must all vanish, 


Hence, we obtain the following necessary conditions for steady oscillations 


COS 
4,d6 


(i 


‘ sin 


On substituting Eqs. (3.4) into Eqs. (3.2) and (3.15) 


V/ 2) (MW? — 1))kA1 + (1 Q)AjA,?? sin 2¢ 
-((/ 2) (MW? — 1) JRA» — (1/Q)A2A;? sin Ze — | 
VAP? — 1 kA, + (A1/O) PD? + 4? — 145 
Vii | uk*A» + (Ao Q)[4A1? + 16.7 — 16 + 4s - 
where Y or 


i be solved for A, A», 
Q. 


These equations are ti 
is the trivial one A; = A 


a flutter mode of finite amplitude is never a purely first or second mode of free vibration in vacua. 


case A, A» + 0 needs consideration. 


QO — (3/4) cos ¢ A,?[(r? + 1)/r(r? — 1)] 
au? — 3(1 — s) = 12 + Ay2f(5/4) + fl — (1 
V M? — 1 pk? = (1/0) {-AP +1 -—s5 + Ar{l 
(\/? 2)/(M? — 1)] —(16 3)(sin ¢ k)|r/(r’ 
where, for simplicity, we have set 
1. = 0. r= A,/A>- (3.17 


In the form (3.16), it is easy to solve the inverse problem 
of calculating \,, k, Q, and ./ from assigned values of 
VM? — 1 A), A», 


and ¢. From a survey of such 


Kt, 
results the direct problem can be solved by interpola 
tion. 

3.1) Buckled Flat Plates 


The flutter motion of initially flat plates has been 
computed, and the results at 1/ = 2 are presented in 
13, 14, 15. for other Mach 
Numbers and mass ratios are quite similar—an example 
at 1/ In these Figures it is 


assumed that 


and Diagrams 


Figs. 
1.3 is shown in Fig. 16. 


Aw dai Po =0, s l (3.18 


In other words, no initial warping and lateral load 
The 


by 


other than the aerodynamic pressure is allowed. 


calculations were made on a digital computer 


assuming various values of ¢, A;, and r and computing 


0, 


yg, and k for a given set of constants VW, uw, Q, and A). 


It is also easy to see that if A, 


When this is the case, Eqs. 


vin cos 
f ) 6.40: 0 ».13) 


‘ sin 


and performing the integration, we obtain 


(S 3)Ao sin ¢ 0 
S$ 3)A1 sin ¢ = 0 | 
>. 14) 
$)A, = 2A 7 | cos o | T S/ <2 | COS ¢ ) 
2A _ 1, cos 2¢ — 12A> } S83 1, cos ¢ () 
— FF 0) @o 5.15) 


An obvious solution 


0, then A» Q, and vice versa. In other words, 
Hence only the 


3.14) may be rearranged as 


cos Pe 4 107 )} + (SO 3 
3.16) 
+ (2 ) + (cos 2¢e 7 i — S'35)(cos ¢ r 
1)] 
Q, A°, Rk, AW? according to Eqs. (3.16). Since only 


real, positive values of Q, \,°. k® have physical meaning, 
not admissible. 
A, > 0, A» 
r>0O. In order that Q may be positive, we must have 
> ifr <1. Note that 


many combinations of ¢, Ay, 7” are 


Without loss of generality we may let > 0, 


1 > w/2 lL, le < 2/2 


¢| > 2 iar 
V \/° — lu enters only in the step of obtaining & and 
MJ, and that the sign of ¢ affects only the Mach Num 
ber. 

Figs. 13 and 16 show the contours of the phase angle 
difference ¢ between the two assumed modes of Eqs. 
(3.3). [see also Eqs. (3.4) and (3.15)], and the contours 
of the amplitude of 
The plot is made on the (Q, A,) plane with ), 

W733 & 
when the assumption of vanishing spanwise 

when /’ E/(l — v*). It 
the M l. appears 
Here we emphasize that 


the antisymmetric mode, A». 
writ 
ten as (initial arch height) (plate thick- 
ness) 
1.€., 
is remarkable that 
similar to that for M/ 


the \/ - 


strain 1s made 


map for 
2.0. 


1.3 map is valid only in the region where k < 1. 
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The mass ratio slightly affects the quantitative dis- Fig. 15 shows the contours for the reduced frequenc 
position of the contours, but not their general charac- k. It is seen that for any given value of Q, Rk increases | it 
ter. as \; increases. pl 

Fig. 14 shows the contours of the amplitude ratio From these Figures we may expect the following 
A, ‘A, of the two assumed modes. It is seen that over behavior if a plate is forced to buckle in a supersoni 
most of the (Q, \;) plane the ratio A, ‘A, is large. For flow. As X, gradually increases, the flutter motion | th 
sufficiently large values of Q, the larger the ), the if it occurs at all, will first have a small reduced fre. | bi 
smaller is the second mode in comparison with the quency, an amplitude ratio A; A» not much greate; | dv 
first. Only in the smaller \, region is A» approaching than 1, and a phase shift ¢ not far from —180°, 4x 
A, in magnitude. the edges are further compressed, so that \; becomes I 
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Contours of the reduced frequency k = w L/U. Mach shift g and amplitude of the second mode As M = 13, ASSU 


Number = 2.0, mass ratio 1/M? — 1] uw = 100. VM? — 14 = 100. 
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larger, the flutter motion will have higher reduced 


frequecy, more predominant first mode, and a larger 
phase shift toward 90° between the two modes. 
- These Figures also throw some light on the flutter in 
the Mach Number range from | to y 2. In this range 
the linearized supersonic aerodynamics dictates insta- 
bility at small reduced frequency due to negative aero- 
dynamic damping. The limiting cycle analysis shows 
that if flutter occurs at small values of Q, the flutter 
motion will have a large A; A» ratio—t.e., the deflection 
will consist of a predominant first mode. 
Unfortunately, at the present time a dearth of ex- 
perimental data precludes any attempt at an evaluation 


of the accuracy of the foregoing analysis. 


3.2) Relation With Stability Boundary 


If we compare the stability boundary labeled \,, = 0 
in Fig. 3 with the flutter maps of Figs. 13 through 15, 
we see that the stability boundary divides each flutter 
map into two regions. In the region below the stability 
boundary, there exists a stable static equilibrium con- 
figuration; in the region above that boundary, no such 
stable static equilibrium configuration exists. It is 
plausible that, in the latter case, any disturbance away 
from the condition B,,» = dB,,» dr = O will cause 
a motion which ultimately will approach the limit 
cycle. On the other hand, when a stable equilibrium 
configuration does exist, a small disturbance about the 
static equilibrium will cause a damped motion converg- 
ing toward the equilibrium point, but a larger disturb- 
ance may cause certain large deviation from the equt- 
librium point and a motion which ultimately converges 
onto a limit cycle. It appears that the wind-tunnel 
experiments described earlier, in which no external dis- 
turbances other than those inherent in the flow were 
imposed, satisfied the qualification of “‘small’’ disturb 
thus revealed the infinitesimal stability 


ances and 


boundary. 


3.3) Numerical Integration of the Differential Equations 


In order to verify that the approximation (3.4) is a 
reasonable one for the differential equations (3.1 
numerical integrations of (3.1) were carried out for 
several values of Q, taking as initial conditions the 
values of A;, Ao, ¢, k given by the approximate solution 


3.16) and shown in the Figures above. Two typical 


sets of curves are given in Fig. 17. The solutions given 


by Eqs. (3.16) for these two cases are 
a QO = 2.03 A, = 2.93 
1, A 14 A 0.6 
¢ = 140 k = 0.378 
M = 1,95 
b Q = 0.524 4, = 2.99 
1,/A2 = 11.0 A 0.6 
a) 100° k = 0.679 
y= 206 


From Fig. 17 it is seen that the sinusoidal-variation 


assumption is reasonable for B,, but not for B.; and 
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Solutions from numerical integration (tep—-case a 


¢? 


om—case b). 
that in both examples B, predominates. Difficulties 
were encountered when the numerical integration was 
continued further; it seems that the rounding-off error 
and the stability of the solution tend to disturb the 
cyclic motion into a transient one. It seems remark 
able that nearly sinusoidal oscillation exists even for 
Q = 0.524, which lies well below the stability boundary 
i.e., in the stable regime—with respect to infinitesi 
mal disturbances 


CONCLUDING REMARKS 


In the present article we have shown some finer 
details about the nature of flutter of buckled plates in a 
two-dimensional flow. It appears that any relief of the 
compressive load in the plate, whether it be due to 
initial curvature, or pressurization, or support elasticity, 
or temperature differential or edge movements, will be 
beneficial to the prevention of flutter. An approximate 
method is demonstrated for the calculation of the 
limiting form and amplitude of the flutter motion, 
although the accuracy of this method remains to be 
evaluated. 

The Galerkin method is shown to give satisfactory 
The 
difficulty of the Galerkin method with the membrane is 


results for plates with finite bending stiffness. 


of course a difficulty in the question of convergence and 
of relying upon a crude approximation to settle the 
question of existence of a real eigenvalue. We may 
recall a famous remark by Richard Courant’ concerning 
the spectacular success of Walter Ritz in inventing a 
method which now bears his name: ‘‘The first success 
attained by Ritz depends largely on his good fortune in 
attacking the seemingly more difficult problem (of the 
free vibration in vacua) of the plate rather than that of 
the membrane.’ Courant has suggested a method to 
“sensitize’’ the variational integral to improve the 
convergence, through the addition of terms of higher 
vanish for the actual solution, to the 


order, which 
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original variational expression. We find that a similar 
improvement can be achieved by deriving the vari- 
ational equation through the complementary energy 
But then in both 


cases the calculation required for the membrane be 


theorem as outlined by Reissner.*° 


comes even more laborious than that for the plate. 
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| Structu ‘al and Impact Loads for the Flexible 


Airplane During Water Landings 


EDWARD WIDMAYER, JR.* ann ROBERT H. SCHWAB* 
The Martin Company 


» October SUMMARY 
: | \{ method for calculating structural loads due to water landing 
nN lane a . ° 
SAS fo impact has been developed The results of this method have, for 
hil oe the most part, been substantiated by the results of an experi 
re mental program in which both the hull bottom loads and the 
“i wis ieceleration response of the M 200 airplane were measured 
Past research has shown that the impact loads are dependent on 
She F rate-of descent, trim angle, landing speed, hull shape, and degre 
Il. Py f flexibility of the entire airplane. The method presented 
M : erein inc ow s the above factors and emphasizes the dependency 
M5 12 f the hydrodynamic loads on the flexible airplane response. 
Ow ‘ Calculations are made for both a rigid and a flexible airplane 
a 2 ind are compared with the experimental results. In_ using 
ie. jodern computing methods it was found that continuous trim 
riations and load input at several stations were possible 
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SYMBOLS 


i = acceleration 
= coefficient of mth vibratory mode 
distance in x direction from mass center of air 
plane to mass station, positive aft 
( = wetted width of hull 
( = pressure coefficient 
D viscous damping coeftlicient 
base of Naperion logarithm 


/ = hydrodynamic force acting through step of air- 


plane normal to keel, positive down 
( = structural damping coefficient 
= acceleration due to gravity 
= draft of hull after chine immersion 
draft of hull at chine immersion 
wetted length of hull measured from step 
= mass considered to act at a given 
station 
Vf = hydrodynamic moment about step, 
nose up 
= virtual mass of water 
»Pmar = Nydrodynamic pressures 
= time 
1" = forward speed 
speed normal to keel 
\ = coordinate in fore and aft direction of airplane 


coordinate in span direction 
= coordinate in direction normal to keel 


A hydrodynamic force associated with dynamic 


pressure, positive down 
Wagner's function for spray thickness 
= density of water 
¢ = slope of deformed structure 
w = circular frequency of mth mode 
T = trim angle, positive nose up 


Superscripts and Subscripts 
(" = time derivations 
= at chine or chine immersion 
Presented at the Aeroelasticity II Session, 
\nnual Meeting, IAS, New York, January 28-31, 1957 
*Group Engineer, Dynamics Section 
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positive 
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average 
may maximum 

at step 
0 = initial condition 
n,? modal designation 
n = normal 


horizontal 


INTRODUCTION 


lpn SUBJECT MATTER of this paper forms a small 
portion of the results of a research program ini 
tiated by the Bureau of Aeronautics to study loads de 
veloped during impact of water landings. (See refer 
It is felt that the results of this work 
represent a significant advance in the treatment of the 


ences | and 2. 


landing impact problem and form a base from which 
analytical treatment may be extended to cover land- 
ings in any type of sea. Historically, the water land- 
ing impact problem has been separated into two dis- 
tinct parts: the determination of the maximum water 
loads, and the determination of the transient structural 
loads arising from the structural response to these 
water loads.’ This separation of the generation of 
water loads from the structural response has been long 
recognized as an artificial device, and a method of in- 
cluding the influence of wing structural flexibility on 
the water loads has been given in reference 4. 


For relatively rigid seaplanes, the transient struc- 


tural loads determined from the maximum water 
loads generally have been sufficient for engineering 
purposes. However, for more flexible seaplanes a 


formulation of the impact problem that allows the full 
interaction of the hydrodynamic forces and the dy- 
namic response of the elastic seaplane is desired. Since 
the dynamic response of an elastic seaplane is a func 
tion of rate-of-load application as well as the peak 
value of the load, the formulation of the impact prob- 
lem should permit rapid solution for a wide parametric 
variation of initial landing conditions 

An analytical method of treating the step landing 
impact problem in smooth water is developed. It 
gives in ‘‘closed’’ form the simultaneous determination 
of both the water loads and the transient structural 
loads. In solving the resulting equations of motion, a 
procedure using an electronic analog computer has 
been developed which permits the rapid solution of the 
landing problem for a large range of initial landing 
conditions. Solutions of the impact problem obtained 
for the M-270 seaplane, considered as being com- 


pletely rigid and as being completely flexible, for a 
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Zz 
Water Level . 


Step Keel Point 


iG 


\irplane coordinate system 


range of landing conditions are presented and dis 


cussed. 


The experimental program for the M-270 airplane 
includes 


is briefly discussed. While this program 


extensive experimental measurements of hull bottom 


pressures, loads, and structural response, only those 
results directly applicable to this paper are presented. 
These experimental results are compared with the re- 
sults of the theoretical analysis. 
are discussed with respect to both the time variation 
and the maximum values of the various parameters. 


THEORETICAL LANDING LOADS ANALYSIS 


The airplane coordinate system is shown in Fig. 1. 
The displacement of any point in the elastic structure 
is considered to be represented by the sum of the dis- 
placements of that point in each mode of vibration in 
cluding the zero frequency modes. This may be ex 


pressed as follows: 


Z(x, y) = d(x) + a,Z,(x, | 
The selection of the modal shapes necessary to approxi- 
mate the structure during impact accurately is gov- 
erned largely by the vibratory characteristics of the 
structure and by the time variation of the applied 
forces. Bases for the inclusion of a sufficient number 
of modes have been discussed in previous literature 
on dynamic response and, therefore, will not be dis- 
cussed here. For the purpose of this paper the dis- 
placement of an airplane element in the elastic case is 
represented by the contributions of the first four 
flexural airplane modes in conjunction with rigid body 
translation and rotation. The displacement of an 
airplane element in the rigid body case is represented 
by rigid body translation and rotation. 


The hydrodynamic forces and moments generated 


during a step landing are given by Eqs. (2) and (3). 
Fe, = —A — M,,(Z, — 0.65/,7) ey 
M, = —F,(0.65/;) 3 


The theoretical developynent of Eqs. (2) and (3) is dis- 
cussed in detail in references 1 and 2 and will not be 
treated here. However, the results of the theory 
necessary for a solution to the equations of motion are 
discussed below. 

The term A is associated with effects of the dynamic 


pressure and is given in Eq. (4). 


These comparisons 


where 


VY Z t 17, sin T 0.65/,.7 


In this expression the influence of the vertical velocity 
of the wetted surface is obtained from contributions 
Z and 7. 


location of the center of pressure of the 


The distance 0.65/, is assumed to be th 
wetted surfac 
from the step of the keel. 
is reflected in the contributions of 1), sin 7 

The functions C, and C have discontinuities at th 
point of chine immersion. 
itil the chines are immersed, and for fur 
draft, C 


Prior to chine immersion, C 


with draft 


ther increase in is assumed to be constant 


equation : 
p 
r sin? 7 + uw? cos? r + (Z, C)? sin? + 


After chine immersion, C, is given by 


Pp 
Pyar sin? c{1 + (Z./C)? + (Z,/C 


In these expressions for C 


function of the draft. Knowing the variation C,, C 
as a function of Z, for a given value of +r. the 
. 


integral C,C(dl,) could be evaluated. 
0 


and /, 


) 


The second term in Eq. (2) 1s associated with the 


forces given by 


well-known “‘apparent mass’”’ effects 


These are the 


With the variation of C with Z, known, Cc 


may be evaluated for a given 7 
The moment about the step is taken as that due to 
the water force acting through the center of pressure. 
The hydrodynamic forces acting on a section of the 
hull are complicated by variations in trim, complex 
effects, 
The treatment of these parameters 


form factors, chine-immersions and three 


dimensional flow. 
may be made by any rationally acceptable methods of 


, It is important that 


analysis (cf. reference 2 or 5). 
rational expressions for force and moment in terms of 
the parameters of the dynamic problem be available 

In the actual evaluation of the above integrals, it 
was convenient to utilize digital machine computa- 
tional methods. The integrals were evaluated for 7 


equal to 3°, 6°, and 9° for the range of Z, expected to 


airplane, and the results 


It was then found to 


be encountered by the 
were shown in graphical form. 
be possible to approximate the curves for constant + 
with analytic expressions in which Z, and 7 appeared 
as simple functions. Further, it was found that the 
family of curves was adequately given when the value 


The effect of planning forces 


The wetted width, C, varies 


is given bv the following 


the quantity Pore/Pmer i8 a 








STRUCTURAL AND IMPACT LOADS 163 
{7 Was permitted to vary. Consequently, when the values of trim is no longer necessary to the solution 
simple analytic expressions for the integrals are sub- of the landing impact problem. 
stituted into the equation for the hydrodynamic force, TI ' pene ee ' et 
ie ; Se DTS le expressions for the C,Cd/i, anc *dl, 
the restriction of the landing conditions to constant Pt me 5 " 
were found to be as follows for the M-270 airplane 
vel . 
_— ; . 40 ~ 80. 
utions of | O< 2, <= O22 C,Cd! x 
J T 
) be the 
1 surfa : ee 1.3! — 2.27 
UTlace 02 < 7 P Z C,Cd , (e22 — 0.82 5 
12g lorces J0 T 
atin ieee S7.1 + 94.5 + 40.1 — 34.2¢7!-5 
S at the Ze . é C,Cd 
. . / 0 T 
Varies 
for fur- ind 
MStant 
lowine r r : 7 _ r r - 
Ing <6 ey Tp Cdl (17.8/7r)Z (Z? + 0.1435 
J 
ya) 
ne 4 ro | C*d 16.8 7) (Z — 0.735 | 
a JV 
In the present problem, these expressions gave excellent 
approximations of the integrals for values of Z, up to obtained 
36 in. and gave reasonable agreement for values of Z y >> 
ae » —L9 Tr (f m C 
in excess of 36 in. It should be noted that for that ont 
portion of the landing most important to the airplane n 2, 3, 4 
is a response—i.e., during the build-up to maximum F.b.+M 
i force —these analytical expressions are most valid. -—T to 0 6 
: : Yom b 
r the The forces arising from the explicit variation of trim ns 
in Eq. (2) include a component due to the dynamic pres- y 5 dn)s Fy + bs M 
i : —Z — 9 i — W, Z T > 0) 
sure and a component due to the virtual mass. In the Dan)? Mm 
h force due to dynamic pressure, |’, has components of ree _— . 
the Spano : oe when the initial values of Z,, 7, and |’, are given, Eqs 
velocity arising from both the trim and the rate of pitch , , : a 
: “sit ’ : ay (6) may then be solved for the Z(t), r(t), F(t), and 
of the airplane. During the early portions of impact i ; 
pedi: ; \J.(t). In addition, the acceleration at any mass 
the contribution due to the instantaneous trim does : oe ; 
i x station may be obtained by the use of the following 
not vary appreciably. However, the contribution due tists . 
‘ , iw : : elation: 
fect to the rate of change in trim, 7, in those cases where 
ects ° 
there is large wetted length, may become important. ; F - 
7 tgs ne ; ‘alia Ziey = AZo + bay 7 + DY a(my),Z, 
C-di, The component of force due to the acceleration of the 
“virtual” mass of water is affected by the angular ate ee ; 
° ss With these quantities known, the impact forces and 
acceleration. In those cases where the wetted length ; ales 
1e to , ‘ i : : structural response are determined, and the distribution 
is long and the pitching acceleration is large, the effect ; ‘ 
re. sepia , ; : ’ i of load over the structure may be studied. 
of this acceleration on the forces is appreciable. . as —m rs 
' the , Brace : ae From an examination of Eqs. (6), it can be seen that 
n permitting trim variation, 1t was necessary to de- ; . . ° ws , ¢ 
plex I ere Species. the equations are nonlinear functions of Z, Z, Z, and 
termine the centroid of the pressure distribution acting ne : ee 
Iree- | : ! r. Previously, when it was sufficient to treat the 
| onthe bottom. For the high-beam loaded hull where ‘ : ~ 
‘ters eee ae ‘ landing problem as the rigid body, the translatory re 
, chine immersion occurs early in impact, the center ot wee . . 
Is of | ; oe ; 4 : sponse to the hydrodynamic forcing function solutions 
pressure is taken at 65 per cent of the instantaneous : : : , : 
that ; ae a , was obtained either by a step-by-step integration or by 
value of the wetted length. The factors on which the , - mn : 
sof | reat ; : use of an analog computer. The inclusion of addi- 
selection of 0.65/, is based are discussed in reference 2. é tie 
ble. es ‘ ; capiiga te tional degrees of freedom has increased the complexity 
rhe location of the pressure centroid is of importance ; . 
, it : ; aed essai of the problem and has made the analog solution 
for landings at low initial values of trim where the 
ita : : attractive. 
wetted length may be expected to exceed four times the es ; . . . 
Es ae ; rhe solution to Eqs. (6) has been obtained readily 
beam of the hull. ee ‘ 
to by using an electronic analog computer. The left 
ilts ; hand side of Eqs. (6), giving the forces arising within 
EQUATIONS OF MOTION oar ° ° 
to the structure, presents no difficulty in machine pro- 
r The equations of motion relating the behavior of the graming. The use of a simple analytical expression for 
ed airplane to the hydrodynamic forces encountered dur- the hydrodynamic force and moments, while requiring 
he ing impact may be derived with any of usual procedures additional machines, may also be readily programed. 


ue of this type of analysis. The following equations are Continuous records of hydrodynamic forces and mo 
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ments, displacement and structural accelerations, 
shears, and bending moments may then be obtained. : Tm 
EXPERIMENTAL PROGRAM j 
A test program was initiated to obtain experimental A ti 
data for correlation with theoretical hull bottom water { fH 
; : . , : atti iti eat fs Jee + o Hy 
loads and airplane accelerations during landing impact | aan oe mneeer.” i 
: ares ra at 
of the M-270 airplane. 
° ° : . ° Anant fo i Nt yae*,. ] 
Instrumentation for the tests consisted of a total of ai fate | NMS NA f tO eg NT i 
a i ‘ , > = ae | \Y (eee ‘ Noe ee >it 
$49 transducers, which included 13 accelerometers, 20 aaass i 
1g S08 { 
‘ a Ai +t toy SOna = 286; a 8 
_ T ose | ae TY { 
— aa Fe sos O17 —— mes L fe 
—{ is eS >, - > 
f ‘ rf o i 
| Ls | bet Hi 
Ci sig cme iy 
~ solutions to impact. problem —flexible case 
2g per sec., m = 3°, Z) = 4 ft. per sec 
e Accelerometers (vertical 
ere 
Strain Gauges 
flush-type pressure transducers, and 16 strain gages 
sid The distribution of transducers in the airplane is shown 
in Fig. 2. The output of the transducers was recorded 
? on oscillographs. 
q The most difficult quantity to measure—the rate of 
-_ é descent —proved to be, perhaps, the most important 
. : ua 
__ Sy a _ _ % Attempts to use doppler radar were unreliable over 
— ; at, Oo : 
=, ot ate a anything but a glassy smooth sea. Consequently, a 
SSS ae ae — photo-theodolite system was installed at the landing 
Station -65 0 100 200 300 400 500 600 700 800 900 1000 1101 range. The rate of descent could be obtained with 
reasonable accuracy from motion picture film of each 
a landing trial. 
a ae i me = ' , “ae 
“LOR 7 —_ The range of parameters investigated was limited 
—— a by the structural strength of portions of the airframe 
Fic. 2. Location of instrumentation on aircraft structure and by the handling characteristics of the aircraft 
Canis The strength of the engine mount structure (limit load 
No Item Location factor of 4.5) placed a limitation on the maximum rate 
2 acceleromete hull sta. 18, wl. 4 of ship ; . . . «pe 
| ene = Ha ate. 18, wh. 09/s, Bol hi of descent obtainable. Also, it was found difficult to 
22 accelerometer hull sta. 396, w.1. 130, € of ship ; i 
23. accelerometer _ hull sta. 885, w.1. 160, € of ship obtain high trim angles at high rates of descent and to 
24 acceleromete ] g sta. 147.5, front sp: : : > 
oe waeace eee ie hatte EE hs ah teed obtain low trim angles at low rates of descent. How- 
25 accelerometer 1. wing sta. 147.5, rear spar ’ : E : : , : 
26 accelerometer _1. wing sta. 600, front spar ever, the flights were performed within these limits and 
27 Ss « Kage nye > o + ~~ - o be . ra . 
; train yage engine mount tube eegeed f tube all possibilities were exploited. 
28 strain gage vertical shear, hull sta. 328, web of floor frame : : : 
29 accelerometer engine mount, I.h. forward Solutions of the equations of motion for hydrody- 
30 strain gage ‘ngine mount tube, center of tube —_— . : 
ah teats ee namic impact have been obtained and are presented in 
31 strain gage sta. 236, center of longeron, r.h. side 
32 strain gage sta. 236, walkway stringer, rh. side Tables 1 and 2 and Figs. 3a to 3e. The solutions, as 
33 strain gage hull sta. 328, € top chord floor frame . . . : : 
revedsanecs coeds <i gil tin alata obtained from the analog computer, are in the form of 
34 strain gage hull sta. 386, center of longeron, I.h. side ‘ ; 
35 strain gage hull sta. 386, walkway stringer, r.h. side the time variation of the quantities being studied. 
36 strain gage hull sta. 386, outboard face of longeron, rt rT . . 
ee a eee eer PR oe These data have been reduced to the form of peak 
37 strain gage hull sta. 386, inboard face of longeron, r.h. side values and are reported in terms of the conditions exist- 
38 strain gage top of keel former, sta. 321 . : ° ° ° 
39 strain gage hull bottom midway between S-0 and 5-1, sta Ing at the time ol peak value =. scillograms ol land 
338 ings are also presented, which representatively illus- 
40 strain gage top of S-O at € of vertical stiffener to stringer eo : 4 ‘ y : 7 
, — trate variations of hydrodynamic force and moment 
attached at sta. 328 - 2 
41 strain gage top of S-0 at € of vertical stiffener between sta motions of the aircraft, and the acceleration response of 
328 and 350 ° ° ° : ome ‘ $ . : ’ 
ee ik a Ae eh lis eatin a ts, 208 the aircraft having significance in the structural dy- 
43 strain gage top of S-2 at € of vertical stiffener between sta namics problem 
328 and 350 
44 accelerometer € hull, sta. 180 
45 accelerometer € hull, sta. 574 DISCUSSION 
46 accelerometer !. wing, sta. 300, front spar 
7 ee 6 ee Oe ee ae It is a well-known result that the response of a simple 
48 accelerometer 1. wing, sta. 450, front spar 


49 accelerometer engine mount, I|.h. forward 


mechanical system to suddenly applied forces is a func- 
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tion of the time rate of force application and the indicial 


admittance of the system; thus, 
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ABRERORAUTICAL 


x(t) = FyAqg + | (dF dt) A(t — t)dt 
JV 


The response of a mechanical system also may be ex- 


pressed in terms of the response of each normal vibra- 
In either ap- 


tory mode to the harmonics of the force. 


proach to the problem, the rate of force application or 
the harmonic content of the force, in conjunction with 


TABLE 
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958 


bearing on the response of the system. 


of high rigidity, that is, having 
quency well above | cycles per sec. (this is the | 
frequency of a harmonic function which will peak at | 
tmar) Of the impact forcing function, the most sever 
distributed structural loads are generally obtained ; 
For the flexible aircraft 


Baas 


the hardest impact condition. 


having a frequency close to or within the frequency oj 
1/4tnar, it is possible for maximum local loads to be ob. | 






the frequency spectrum of the system, has an important ] 
For aircrait 
a fundamental fre. } 























Analytical iia tcp ameiiili \irplane 
V, 1% % F(lb) M(in.-lb) 2Z(in)  1(rad) ¢ < Ay (8) Ags Ap, Aqy Ain 2aeFuax Tr Landing 

maxX No. 

185 3 & 25,000 3.83x10° 12.67 2 .10 .32 31 40 8.95). TTS -37 12.0 3h 1 
8 64,000 13.0x10° 19.67 .20 32 61.28 875 1.2 2.75 3.15 95 12.0 .20 2 

12 110,000 23.33x10° 28.33 .20 55 2.4 1.00 1.78 5.38 5.00 1.50 13.33 .125 3 

6 & 431,000 2.13x10° 10.16 09 =. 031 -06 .68 55 <9 207 -76 10.0 .325 4 

8 90,000 9.33x10° 15.33 oe -06 48 21.38 «611.6 2.40 4.0 1.63 15.33 225 5 

12 150,000 18.0ox10° 19.66 .198 18 1.16 2.60 2.45 5.0 6.88 2.95 16.67 15 6 

160 3 4 23,500 5 .83x10° 13.73 138 838.115 36 .28 +42 .90 .688 33 13.33 45 7 
8 60,000 12.0x10° 21.73 198 «3461.2 -875 1.10 2.56 2.90 1.00 15.33 .20 8 

12 104,000 21.33x10° 32.0 215 .565 2.26 1.05 .170 5.20 4.50 1.00 16.67 .125 9 

= 16 153,000 35.33x10° 4.33.20 1665 3.5 3.46 2.40 8.50 6.50 4.85 16.67 .10 10 
6 4 28,500 53.33x10° 33.55 09 .017 .048 Sh 49 6635.95 65 11.00. .325 u 

8 83,000 9.13x10° 17.07 48.084 52 1.32 1.48 2.40 3.66 1.60 15.33 20 12 

12 130,000 16.33x10° 21.7 21 22 «61.16 2.06 2.30 4.15 6.0 2.60 18.0 15 13 

16 180,000 25 .07x10° 27.5 21.374 = 2.0 9.5 3.2 7.65 7.5 5.10 20.83 .125 14 

9 & 31,000 1.47108 10.27 12 225 As! .80 55 .56 1.05 85 10.0 .325 15 

8 83,000 5 .93x10° 16.27 135 .038 16 1.70 1.50 1.75 3.55 1.87 16.0 .225 16 

12 135,000 11.83x10° 20.67 Be i -O4 43 2.94 2.35 3.60 6.20 3.20 17.33 16 17 

16 185,000 18.53x10° 26.67 .20 1h gh %&.88 3.0h 6.5 7.0 6.0 26.0 .15 18 

135 3 4& 23,400 4.1x20° 14.67 15 = .133 36 2h 39 ©.65— 69 28 14.0 40 19 
8 55,000 11.13x10° 25.33 .21 36 1.12 80-90 -2.25 2.80 90 16.0 .20 20 

2 95,000 20.33x10° 39.17 20 «6.525 2.05 55 1.45 4.6 4.15 -75 14.17.1225 21 

16 140,000 31.67x10° 53.33 21 60 3.25 3.5 2.20 7.68 6.0 4.5 16.67 -10 22 

é 4 26,000 2.13x10° 12.27 .105 =.009—s« «072 50 45 .75 83 -57 12.00 35 23 

8 72,000 8.53x10° 18.67 “AT Bs S52 3435 1.50 2.2 3:35 2.35 27550 125 24 

12 120,000 16.4 x10° 23.33 .20 25 1.22 2.30 1.96 4.15 5.25 2.90 20.00 .175 25 

16 160,000 23.67x0° 33.33 .20 40 1.90 3.06 2.70 6.9 6.65 4.25 23.33 15 26 

9 4 29,200 1.60x10° 11.87 23 .051 .086 65 Sl .525 95 69 11.8 .375 27 

8 80,000 6.13x20° 17.60 IS 019 <2h 1.72 «230 «176 3.2 .196 16.67 .225 28 

12 120,000 1.2 x0® 23.33 .194 .075 .472 2.8 1.92 3.4 5.25 3.10 20.00 15 29 

16 162,000 17.330° 30.0 20 .188 96 3.8 2.64 5.75 6.5 4.70 23.33 15 30 
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‘ntal fre. 
iS is the 

peak at 
st 


tained jr 





Severs 


> aircraft 


uency of 
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STRUCTURAL 
tained for impacts other than those having the greatest 
hydrodynamic force. This condition requires that a 
range or envelope of impact conditions be examined in 
order to insure that the most severe loads are utilized 
in the design of the aircraft structure. 

Table 


with variations in the rate of descent. 


1 shows that a variation of f was obtained 
The range of 
_,r Was from approximately 0.10 sec. for the high rates 


of descent to approximately 0.40 sec. for the low rates 


AND 


[MPACT LOARPS 


The spikes appearing in these curves resulted from 
switching equations in the REAC and may be disre 
It may be noted that the general shape of 
similar for different 
Also, it 


may be seen that some delay in the generation of the 


garded. 
the forces and of the moments are 
impact conditions with only ¢,,.; changing. 


moment with respect to the generation of the force was 
experienced. This delay is caused by the hydrodynamic 


center of pressure being close to the airplane center ol 












































of descent. Figs. 3a to 3e present the shape of the force gravity for short wetted lengths. As the wetted 
and moment generated in representative landings. lengths become greater, the moment increases at a 
TABLE 2 
\nalytical Data for the Rigid Airplane 
Y, 1% Z, F(lb) M(in.-lb) Zin) t(rad) t 7 Aj, (8) Ay Aga Any, Ain ZatFmex oo imei 
185 3 4 25,600 3.93x10° 12.93 12 .099 32 20 45 (BS 45 20 12.0 35 31 
8 68,000 13.67x10° 20.0 205 33 1.30% Te 1.225 2.7% 1.25 -75 16.33 20 32 
12 115,000 23 .6x10° 30.0 205 .575 2.40 52 1.97 %&.75 1.97 -52 18.33 1 33 
6 4 32,600 2.27%10° 10.33 .09 .026 06 595 58 .60 .575 60 10.0 30 34 
8 98,000 10.0x10° 16.0 14 .063 -504 1.21 1.70 .565 1.70 1.25 15.35 20 35 
12 157,000 18 .67x10° 21.0 : 18 .188 1.184 1.65 2.76 4.00 2.8 1.60 18.33 15 36 
160 3 4 23,800 4 .08x10° 14.33 13 =3=.5 376 .166 45 87 Kb 18 13.87 4o 37 
8 61,200 12.33x10° 22.67 -205 34 1.20 Teo 6.1.30 O.46 .122 72 16.00 20 38 
12 102,000 21.33x10° 33.33 195 55 2.20 45 1.8 4.2 1.82 45 18.33 15 39 
16 148,000 32 .67x10° 46.67 .20 665 3.36 63 26 64 2.7% 63 18.33 10 4O 
6 4 29,800 2.27x10° 11.33 .20 012s .O48 50 53 59 «55 -51 11.33 35 4} 
8 88,000 9.33x10° 17.60 15 .083 = .52k 1.03 1.54 2.18 1.52 1.05 17.00 225 42 
12 138,000 17.07x10° 22.67 21 218 1.15 1.48 2.4 3.66 2.46 1.48 18.66 15 43 
16 190,000 25 .07x10° 29.0 -222 36 1.9 1.90 3.36 5.50 3.36 1.90 17.66 -10 4, 
9 4 33,000 1.53x10° 10.4 .124h 061 nh 73 60 .495 .60 Te 10.33 30 4 
8 91,000 6.53x10° 16.4 14 033 16 1.60 1.59 1.64 .163 1.62 16.0 25 46 
12 143,000 12 .67x10° 21.2 177 ~—- .OK6 452 2.30 2.56 2.94 2.60 2.32 19.73 175 47 
16 190,000 18 .67x10° 26 .67 196 «6.147.912 3.00 3.44 4.35 3.48 2.96 21.67 .125 48 
135 3 k 24,000 4 .23x10° 15.73 122 13.384 16 4O =. BAS -16 14.67 4c 49 
8 56,000 11..33x10° 26.67 20 «3 1.12 675 95 2.22 1.03 -70 16.67 20. 50 
«ga «ga,000 19.73x10° 2.0 20 ~5¢ 2.0 4) 861.60 3.64 1.72 4h 18.33 9.125 51 
16 138,000 30 .00x10° 56 .C 195 60 3.06 60 2.34 5.85 2.46 .60 16.67 10 52 
f , 27,400 2.33x10° 12.67 l= .009 072 42s 485 .57  .50 425 12.33 35 53 
8 74,000 8.67x10° 19.07 155 mt 5k 85 1.35 .19% .138 85 18.0  .225 54 
12 120,000 14 .93x10° 25.33 195 246 1.088 1.2 2.10 3.2 2.20 1.25 6.3535 .125 55 
16 165,000 22.4 x10° 53.355 21 40 =1.76 1.66 2.96 4.75 3.04 1.74 18,33 10 56 
9 4 30,600 1.57x10° 11.60 me 05 084 -64 54 0475) W555 65 11.67 35 57 
8 80,000 6 .067x10° 18.27 15 018-~—s«w.26 1.55 1.42 1.52 1.42 1.36 17.73 25 58 
12 120,000 11..2x10° 23.73 .185 .078 .496 1.90 2.20 2.60 2.26 1.90 19.33 15 59 
16 164,000 16 .67x10° 30.00 .197 19 936 2.44 2.96 3.75 3.00 2.42 22.67 .125 60 
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greater rate and reaches a maximum at approximately 
the the Consequently, the 
loads induced by both the force and the moment must 


same instant as force. 
be considered in the landing analysis 

Figs. 3a to 3c show that the frequencies excited in the 
structure are a function of f¢,,,, and, hence, the rate of 
loading. In case 3a, the hull may be seen to respond at 
both 140 cycles per min. and at 390 cycles per min., 
the first and third modes of vibration. In case 3b it 
may be noted for the shorter ¢,,,, that the hull responds 
principally in the second vibratory mode at 354 cycles 
per min. In case 3e having ¢,,,; of 0.1 sec. the hull re- 
sponse was at 510 cycles per min. which is between the 
third and fourth modal frequencies. 

The maximum accelerations at various locations on 
the airplane are given in Table | for each landing con- 
dition, while actual shears and moments were not ob- 
tained in this analysis. (The shear and moment dis- 
tributions as a function of time could have been ob- 
tained with the equipment available, but the data ob- 
tained illustrate the point.) It may be seen that the 
maximum positive or negative accelerations were not 
necessarily obtained at the impact having the maximum 
force. For instance, the maximum positive bow and 
wing tip acceleration were obtained for landing 2, while 
the maximum positive stern and stabilizer tip acceler- 
ations and the maximum negative bow acceleration 
were obtained for landing 10. In landing 18, having 
the maximum force only, maximum negative acceler- 
ations were obtained at the wing tip, the stern, and the 
center of gravity. These results indicate that even for 
the rather rigid M-270 airplane an envelope of landing 
impacts would be necessary for accurate estimation of 


the landing loads. 


IC 
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\ comparison of the maximum forces obtained 


the flexible airplane analysis with the maximum forces 


obtained for the rigid airpl ine analysis affords an ¢ por 


tunity to study the influence of flexibility on the hydr 


dynamic loads. Comparing column 7 of Tabl vit 
column 7 of Table 2, it may be seen that in most i 
stances the maximum forces for the flexible airplar 
were less than the maximum force for the rigi ir 


plane for the same initial conditions For the amount 
M-270 airplane the 


mun reduction in force was approximately 7 per cent 


of flexibility present in the 


For more flexible airplanes, particularly if considerab] 
weight is carried in the wings, it is conceivable that 
larger reductions in applied force might be 
Chis interaction between the response of the airplan 
and the hydrodynamic force, and the resulting reduc- 
tion in applied load may have an important effect o1 
the distributed loading experienced by the aircraft 

In regard to the influence of flexibility on the hydro 
dynamic moment, no definite trend is exhibited by thy 
data. 
on the force and wetted length, there is the possibility 


Since the moment depends among other things 


of compensating variations in both factors for the rigid 
more flexible airplanes, appreciabk 
When 
such information becomes available more definite con 


airplane case. Foi 


changes in pitching moment may be obtained. 


clusions may be made on the influence of flexibility o1 
the pitching moment. 

At this point, some comment on the assumptions and 
neglected quantities is in order. First, consider th 
assumption that the distributed hydrodynamic forces 
are given entirely in terms of the kinematics of th 


step. For the M-270 airplane this assumption appears 
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valid, since the rigidity of the hull and vibratory modal 
shapes were such that adequate average values of co- 
efficients could be used. However, for the more flexible 
hull lav be necessary to describe the forces in terms 
f th inematics of more than one hull station Lhe 
use of many hull stations, analogous to aerodynamic 
strip theory, would permit more accurate application 

~ e factors and of local kinematics. However, 
the increase in number of hull stations is accompanied 
by an increase in computational difficulty and for ana 
log methods may exceed the available machine cap iCi- 
ies 

Second, in the formulation of the hydrodynamic 
forces and moments tivo quantities have been considered 


; negligible. ind 


moment and pitching moment due to water resistance 


These quantities are buoyant force 


to forward motion of the airplane. In addition, it has 
been assumed that the weight of the airplane is equal 
to the lift for the duration of impact. Some consider- 
ition of the effect of loss of lift in landing impact has 
been given in reference 6. The work of this reference, 
while treating the problem only approximately, indi 
cates that some account for the effects of lift variation 
may be required. For impacts of long duration or for 
flexible highly coupled lifting surfaces, the lift variation 
might enter appreciably into the dynamics of the 
landing. 

As landing impact becomes more critical in deter- 
mining the design loads for aircraft structure (in addi- 
tion to the hull bottom) and as more powerful compu- 
tational methods are used, these forces may be included 


in the analysis. 


COMPARISON OF EXPERIMENTAL AND THEORETICAL 


DATA 


Briefly reviewing, the landing impact analysis was 


performed for symmetric step landings in smooth 


water. Among the major assumptions of the analysis 
is that the only hydrodynamic forces acting on the air 
plane would be approximated by an analytical expres- 
sion which is a function of the forebody shape and the 
kinematics of the step. Also, it was assumed that the 
dynamical response of the airplane is given by the 
rigid-body response and the superposition of the re- 
sponse of the first four flexible-airplane modes of vibra- 
tions. 


The vibratory analysis of the airplane was 


performed for a gross weight of 55,000 Ibs. It was found 
iter completion of the vibratory analysis that certain 
changes in 
landing tests to be performed at gross weights of 64,000 
Ibs. and 71,000 Ibs. As a 


is to be made in interpreting the comparison of 


the experimental program required the 


consequence, some correc 
tion 
the analytical results with the experimental results. 
A suggested basis for this correction is in terms of the 
peak loads that might be expected for the rigid airplane 
as given in reference 6. The loads in accordance with 
this correction would vary as the ratio of the gross 
weights taken to the two-thirds power. 


The experimental program was formulated for land- 
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Fic. 6. Experimental vs. REAC 


and landing 


ings at discrete rates of descent, trim angle, 
speed. However, in actual practice, it was found that 
descent and trim angle 


REAC 
there is 


only broad ranges of rate oi 


could be obtained. The analysis is a result of 


discrete initial conditions; no direct corre 
spondence between the analytical and nominal experi 
mental initial conditions. One exception to this was 
landing 4 of flight 29, where all the initial conditions 
were duplicated in the REAC solution. ol 
scent (Z to the earth and not 


the water surface, so that even | 2 to |-it. 


Rate de- 
was measured relative 
waves or 
swells could increase or decrease effective Zp. Analyti- 
cal (REAC) data are compared to experimental trends 
rather than individual test points for this reason. 


The comparison of the time variation of the analyti 
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cally determined quantities with the experimentally 
determined quantities for two landings may be made 
in Figs. 4 through 7. These Figures show the hydro- 
dynamic force, pitching moment, trim pitching velocity, 
normal accelerations, and, in one instance, the pitching 
acceleration. In Fig. 4, the analysis was performed 
for the initial conditions of landing 4, flight 29. This 
landing represents a rather hard impact, with the vari- 
ables being of sufficient magnitude that they could be 
determined with some precision. It is seen that the 
analytical results are in excellent agreement with the 
experimental results. The two curves are almost coin- 
cident up to fn, and are reasonably close after fya:. 
Excellent agreement was also obtained between ana- 
lytical and experimental results for 7 and 7 over the 
greater portion of the impact; some divergence was 
observed after ¢,,., had occurred. In regard to the 
pitching moment, it may be seen that there is a con- 
siderable difference between the curves over the range 
225 sec. No apparent 


of time from 0.06 sec. to 0. 
reason has been found for this discrepancy. The agree- 
ment for the pitching acceleration is considered to be 
good in light of the manner of measurement. The 
shape of the analytical curve coincides reasonably with 
the experimental curve and the maximum values are of 
the same order. A comparison of the analytical acceler- 
ation response of the structure as depicted at three 
stations with the measured response at those stations 
is shown in Fig. 5. The experimental acceleration re- 
sponse is seen to be predicted quite well by the analytical 
response. The absence of the higher harmonic con- 
tent in the analytical response may be attributed to 
using only four vibratory modes in the structural rep- 
resentation. As may be seen, this harmonic contribu- 
tion while adding some distortion to the response does 
not materially affect the correlation between experi- 
ment and analysis. 

Landing S of flight 31 for a gross weight of 72,000 
Ibs. also offers an opportunity to compare analytical 
and experimental results. Referring to Fig. 6, again it 
may be seen that the shape of the experimental curves 
for both the force and moment is given by the analytical 
curves. Applying the previously suggested correction 
for gross weight, the comparable magnitude for the ex- 
perimental force for the 72,000-1b. airplane corrected to 
59,000 Ibs. is given by S4 per cent of the magnitude. 
This reduction brings the experimental force magni- 
tude into excellent agreement with the analytical force 
magnitude, particularly at the maximum values. The 
application of the correction to the moments also tends 
to bring the moment magnitudes into agreement with 
the analytical moment magnitude. The 7 variation 
for the analytical impact is in agreement also with 
that of the experiment. Comparison of the pitching 
velocities showed little or no agreement, and no com- 
parison of pitching accelerations was attempted. 

The structural acceleration response for the 72,000- 
lb. airplane may be seen in Fig. 7. Recalling that the 
difference in gross weights of the airplane is achieved 
by the addition of concentrated weights in the hull 
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about the airplane center of gravity, it is to be expected 
that the vibratory characteristics of the hull would be 
greatly influenced, and the wing to a lesser degree, by 
the presence of these weights. Comparing the experi 
mental acceleration response of the 72,000-Ib. airplane 
with the calculated response for the 55,000-lb. airplane 
it may be seen that this occurs. The shape and magni- 
tude of the calculated center of gravity and wing ti 
accelerations are in reasonable agreement with the ex- 
perimental accelerations. The experimental bow ac- 
celerations, apparently showing the influence of the 
behave somewhat 
In light 


added hull weight, are seen to 
differently from the calculated accelerations. 
of the above discussion, it is felt that where some basis 
for comparison with the heavy airplane existed, the 
analytical results are in agreement with the experi- 
mental results. 

Figs. 3a to 3d present the maximum values of force 
and of structural acceleration from the analysis as a 
function of rate of descent for the different initial 
trim angles and forward speeds. For purposes of com- 
parison, experimental values having approximately the 
same initial trim angles and forward speeds are also 
shown against rate of descent. The analytical results 
for the structural response of the structure are in good 
agreement with the experimental data. Reduction of 
experimental data for the hydrodynamic force was not 
done for all landings. However, at a forward speed 
of 160 ft. per sec., the agreement between the analyti- 
cally predicted force and the experimentally determined 
force is good in trend and magnitude. 

The above comparisons of experimental data with 
analytical data have shown that the analytically deter- 
mined results have described accurately the landing 
impact conditions. The correlation of hydrodynami 
forces and moments has indicated that the representa- 
tion of these forces and moments in the equations of 
motion has been substantiated. The structural re 
sponse as given by the equation of motion agrees very 
well with the actual structural response of the airplane, 
and indicates that the flexibility of the airplane was 
adequately represented by the flexible modes of the 
analysis. The response of the rigid airplane in pitch 
given in the analysis was in agreement with the expert- 
mental observed response. In light of this, it is be 
lieved that the hydrodynamic step impact and the dy- 
namic response for the flexible seaplane may be accu- 
rately given by the method of analysis presented in this 


paper. 


CONCLUSIONS AND RECOMMENDATIONS 


A theoretical method of computing the dynamical 
impact landing, including the effects of structural 
flexibility and continuous trim variation, has been 
developed in ‘‘closed’’ form. This method, utilizing 
an electronic analog computer, permits rapid study of a 
wide range of initial conditions for the impact problem. 


The dynamic acceleration response of the M-270 


(Continued on page 216 
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The Compressibility Rule for Drag of Airfoil 
Noses 


R. T. JONES* ann M. D. VAN DYKE* 
Ames Aeronautical Laboratory, NACA 


SUMMARY 


iirfoil section im 
the Prandtl- 


It is shown that the drag of any semi-infinite 


purely subsonic inviscid flow follows precisely 


result for the parabola has 


Glauert compressibility rule. The 
ipplication to leading-edge corrections in thin airfoil theory 


INTRODUCTION 


S MI-INFINITE BODIES have long been used to repre 
sent the noses of elongated aerodynamic shapes. 
This role assumes increased importance in compressible 
flow, because the approximate linearized theory must 
usually be employed, and it breaks down in the neigh- 
borhood of stagnation points.’ It can be corrected, 
however, if the flow is known past a simple semi 
infinite shape (such as a parabola) that matches the 
actual body near its nose.’ Except inl incompressible 
flow, the semi-infinite shape itself cannot be treated 
exactly, but only by another approximation such as 
the Janzen-Rayleigh method. However, it 1s shown 
in the present paper that for the drag itself the effect 
of compressibility is known exactly. 


DRAG OF A PARABOLA IN SUBSONIC FLOW 


As discussed by Prandtl and Tietjens,* the drag of a 
semi-infinite body in inviscid flow is understood to be 
the limit of the pressure drag on a finite section isolated 
by a transverse slit into which the surrounding pres- 
sure penetrates (see Fig. 1). With this definition, the 
drag (per unit span) of an infinite parabolic cylinder in 
incompressible flow along its axis is tR times the dy- 
namic pressure, R being the nose radius. 

As the free-stream Mach Number increases to any 
value short of unity, the entire flow field remains sub- 
sonic and hence free of shock pa- 
a retarding influence on the flow. 


waves, because a 
rabola exerts only 
Local pressures do not, of course, follow the Prandtl 
Glauert compressibility rule even approximately near 
the stagnation the particular 
weighted average of surface pressures that yields the 
Thus, at any free 


point. Nevertheless, 


drag does obey that rule exactly. 
stream Mach Number .J/ less than unity, the drag is 


D (1/2)pl? (#R/V 1 — M 1) 


The proof can be based on consideration of any thin 
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airfoil (say, a symmetrical Joukowski section) of thick- 
If the leading edge is analytic, it has a nose 
Now the airfoil 


divided into a short nose section of length 7 and a re 


ness T. 


radius proportional to Tr suppose 


maining tail section (see Fig. 2), and consider the con 
tribution of each to the total drag. The drag vanishes 
as 7 approaches zero, but the drag coeflicient remains 
as the nose 


finite if referred to a length of order 7*, such 


radius. The nose section disappears in the same limit; 
but if it is magnified so that its nose radius remains 
constant, it approaches an infinite parabola 

Over the rear section the slope is uniformly small, so 
that the linearized theory becomes exact in the limit 
giving a drag coefficient proportional to (1 \/*) : 
Over the nose section linearized theory cannot be ap 
plied. However, if the flow is purely subsonic, the 
total drag is zero according to Theodorsen’s extension 
d’Alembert’s paradox Hence 
as (1 \/*) ae 


The same compressibility rule holds for the well 


to subsonic flow of 


the drag of the nose also varies 


known leading-edge suction force associated with angle 
of attack. This is clear from a repetition of the pre 
ceding argument with the airfoil considered to have an 
angle of attack a that tends to zero in proportion to the 
obtain asymmetric flow 


thickness 7. In the limit we 


= 
10 
] 1g. 


Qne cannot 


past an infinite parabola (see 
speak of the angle of attack, because the slope of all 
streamlines tends toward zero far upstream, decaying 
finally like [s, 
length. However, 
degree of asymmetry of the flow field far upstream by 


1/2 . . . 
—x)]~ where sy 1s some characteristic 


one can instead characterize the 


specifying this parameter so in the dimensionless form 


so R. 


(In the limiting process just envisioned, the 














airfoil 


Flow past a semi-infinite 
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Fic. 2 Division of airfoil into nose and tail sections 


ratio a 7 that was held fixed is proportional to V so R.) 
As indicated in Fig. 3, the stagnation point lies a dis 


tance sy downstream of the vertex in incompressible 


flow, but moves nearer the vertex as the Mlach Number 


is increased. Then for a fixed flow pattern far up- 


stream, corresponding to constant so, the drag of a 


parabola is given by 
D 1 2)pU?[r(R — so) V1 — JZ] 2 


This rule now holds only below the critical Mach Num- 
ber, which is less than unity for asymmetric flow. At 
speeds above the critical Mach Number, local super 
sonic zones appear which undoubtedly terminate in 
shock waves, and these invalidate the extended d’ Alem- 
bert paradox on which the argument rests. It may be 
noted in Eq. (2) that even in nonlinear flow the drag 
terms are additive. 

The rule provides an integral check on detailed ap- 
proximate solutions. For example, Imai has calcu- 
lated the Janzen-Rayleigh expansion in powers of 1/° 
for the surface speed on a parabola in symmetric flow, 
retaining terms in ./*. Carrying out the complicated 
integration for the drag gives 


D/(1/2)pU°xrR 1 + (1/2) M? + (3/8) Mt +... (3) 


This is just the expansion to order .\/* of the Prandtl- 
M7?) 


with the 


1/9 


Glauert factor (1 — ~, so that Imai’s approxi- 


mation is in accord present rule. Again, 
Kaplan® and others have treated the inclined ellipse by 
the Janzen-Rayleigh expansion to order \/*, and from 
that has been extracted the series to ./° for asymmetric 
flow past the parabola.’ Integrating for the drag gives 
again the compressibility factor [1 + (1 2)1/? +... ]. 
To obtain sufficient accuracy for purposes of correcting 
compressible thin-airfoil theory for lifting round-nosed 
wings,” this last solution should be extended to include 
terms in .J/*; when that is accomplished the present 
rule will provide a useful partial check on the formidable 
calculations involved 


DRAG OF OTHER SEMI-INFINITE BODIES 


It is clear that the proof given for the parabola will 
not be invalidated by local changes from parabolic 
shape near the vertex. Hence any semi-infinite cylin 
der that approaches a parabola far downstream has 
(below its critical Mach Number) the drag of that 
parabola. For example, in Helmholtz’ solution for 
incompressible free-streamline flow past a plate of unit 
height normal to the stream, the dead-water region 
far downstream approaches a parabola of nose radius 


2/(r + 4). According to Eq. (1), the dead-water re 
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gion replaced by a solid body would have a drag 
rol (mr + 4), and this is in fact Helmholtz’ value for 
the plate. For the solid body this value would ris 
with Mach Number as (1 — _l/*) 

of sound was attained on the surface. 


until the speed 
Furthermore, the proof indicates that the rule applies 
as well to a semi-infinite body that ultimately grows 
more slowly or more rapidly than a parabola. How 
ever, a body has zero drag if it grows more slowly and 
infinite drag if it grows more rapidly. Hence only the 
dividing case of a body asymptotic to a parabola pr 
vides a significant result. 

Similar considerations could be applied to a thre 
dimensional semi-infinite body —for example, a body of 
revolution. However, the significant shape that sep 
arates axisymmetric bodies of zero and infinite drag 
has been found by Gurevich? to be a peculiar one whos 
radius varies asymptotically as x! *(In x) ?) As in 
plane flow, its ultimate form must be the same as that 
of the dead-water region in axisymmetric free-stream- 
Levinson’ has indeed found the 


line flow, for which 


same form. For a body of noncircular cross section 
since the drag depends only on the shape far down 
stream, where slender body theory becomes accuraté 
the area rule shows that the drag depends only on the 
cross-sectional area. Hence, three-dimensional 
semi-infinite body having finite drag in subsonic flow 
nal 


lidl 


any 


grows with cross-sectional area ultimately proportio 


to x(In x) 


APPLICATION TO SUBSONIC LEADING EDGES 


Thin airfoil theory breaks down near the leading edge 
of a wing in subsonic flow or in supersonic flow with 
the edge swept behind the Mach cone. For round 
edges, the error in pressure is such that the drag can 
be found correct only by including a leading-edge force 
associated with the singularity.!_ This force is just the 
drag of a parabola having the same nose radius, given 






by Eq. (1). (On swept edges, U’ and lJ are to be re- 
placed by their components normal to the edge.) ‘The 
Continued on page 180 
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1) INTRODUCTION 


- ESTIGATIONS of the heat transfer in incompressible, 
laminar boundary layers must usually proceed by 
means of numerical solutions, due to the complexity 


of the governing differential equation. Exact solu 
tions in closed form have generally been obtained only 
when the Prandtl Number o is equal to one. Analyti 


cal results of an asymptotic nature have been obtained 
for large the approximation being based on the fact 
that the 
than the viscous boundary layer 
In the 


problem for plane boundary-layer flows when 


} 


thermal boundary layer is much narrower 
when o 1s large. 


present paper we consider the heat-transfer 
is Very 
small and obtain the first two terms of a series repre 


sentation for the quantities of interest. The procedure 
is based on the fact that 


narrower than the thermal layer when o is small, so 


the velocity laver is much 


that the velocity field is given essentially by the po 
tential flow solution over the major portion of the 
thermal layer. This permits a great simplification oi 
the coeflicients of the differential equati n. 


\ general solution is obtained for arbitrary potential 


flow and arbitrary wall temperature in the cooling 
Received April 19, 1957 

results presented in this paper were obtained in the 

f research sponso! the Office f Naval Research 


t\ 





Contract Nonr 56 with Brown Universi 
ciate of Applied Mathem 
Research Assistant in Applied Mathematics 


Applic d 


Professor tics 


Mathem 
nd Mate 


Formerly, Research Associate in ities 


N Assistant Professor, Department of Mechanics a 


rials, University of Minnesota 


arbitrary potential 
The 


of integrals involving the known 


problem and tort 
lated wall problem. solutions are given in 
velocities and, in the 
the prescribed wall tem 
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form results are determined tor wedge 
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temperature power of the distance fron 


Varying as a 
the leading edge, or with insulated wall. 
Prandtl Numbers are, for the 


The in portance ol the heat 


Fluids having small 
most part, liquid metals. 
transfer problem for such liquids is growing rapidly, 


nuclear reactor 


especially due to the development oO 


technology. The problem presents many hitherto un 


} 


solved questions concerning the physical properties of 


the flow, such as the extent to which the liquid wets 


the bounding surface and the appropriate boundary 


condition for the temperature. These questions art 


not considered here, the investigation being based on 


boundary conditions 


fluid 


usual 
that the 


the assumption that the 


also assumed properties 


apply. It is 
mav be considered to be constant 


2) GENERAL THEORY 


Che dimensionless energy equation for steady plane 
a fluid with constant properties expressed 1n 


parallel and 


flow of 
terms of dimensionless coordinates, (A 


al to the boundary, respectively, is 


norn 
oO] O7 [ ar. 7 
Ky ou - & Rol (1+ Ky)? O 
A O l \ O / ] 
A Uv A J A 1 R 
r] A \ / O \ 
2 wae a oe 
fy 6) ra j/ YU 
/ { 
| O! O \ 
LO | K O } 7 
wilere¢ 
l , a 1 1 ] A7 
i y hits Ul 


U and 7. being the constant free-stream velocity and 
temperature, ;, 7 
and perpendicular to the wall, 7) = 


AT 


the velocity components parallel 
the temperature, 


a characteristic temperature difference, L a 
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characteristic length, R = the Reynolds Number, ¢ = 
the Prandtl the kinematic viscosity, 
and A(x 

This equation, together with appropriate boundary 


Number, v 
the curvature of the wall. 


conditions, will give rise to a boundary-layer problem 
provided that the coefficient Ro 
thermal layer is thicker than the viscous layer, then u 


is large. If the 


will be of order one in the thermal layer and a stretching 
of the normal coordinate y by the introduction of a new 
variable & (Ro)' *y will lead to an equation in which 
both the convection term uw O7'/Ox and the leading 
Hence, 


conduction term 0O°7 O& are of order one. 


using the usual boundary-layer arguments, the thick 

° . - 1/2 
ness of the thermal layer will be of order (Ro ’ 
'?) the order of thick- 


Thus, the 


which is large compared to R 
ness of the viscous layer, whenever ¢ < 1. 
velocity field will be essentially potential throughout 


most of the thermal layer, and we therefore write 

Ue, V Vy A Y, 2.2 
where i), UV, are the velocity components of the poten 
tial flow solution and u,, v, are the corrections which 
must be added to give the components for the viscous 
dealing with a boundary-layer 


flow. Since we are 


flow, u,, v7, are taken to be determined approximately 
by boundary-layer theory. 

Since Ro 
pared to a characteristic geometric dimension of the 


1, the thermal layer is very thin com- 


problem and hence the potential velocity components 
vary only slightly across the width of the thermal layer. 


We expand them in power series of y about y 0; 
Uy, = U(x) + yu) (x) + a a 
: (23 
Up yo (ve) + ...J 


v(x) being zero since the potential flow satisfies the 
condition of zero normal velocity at the boundary. 
Since the potential flow satisfies the continuity equa- 


tion 
(Ou,/Ox) -- (Ov,/Oy) = 0 


we have u\°’’ -7'") so that 


—yuO’4t... (2.3a 


It will be convenient to introduce the variable 7, 
which is the stretched normal coordinate appropriate 


to the viscous boundary layer; if we set 


1/2 
n Ky, ule, 9 u(x, ¥)q a 
ie o> 
1/2 2. 
i(x, 7 R’'“v(x, 
then “#,, 3, and 7 are all of order one in the viscous layer, 
and a, vanishes as one leaves the viscous layer. 
, 1/2 92 
We now set &£ = (Ro) ‘“y (2.6 
substitute Eqs. (2.3), (2.3a), (2.5), and (2.6) in Eq. 


2.1), and neglect all terms which have a coefficient 


(Ro) ‘'~ or a higher power thereof. This yields 


; a” : pa bere oi. 2002 
(u + ,.) (O72) Ox) + (—£u +oa‘d,.) X 


= (0°7/O&) + (U?/c, AT) (O%,/On)? (2.7) 
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We note that to the above order of approximation the 
potential velocity component parallel to the wall is 
approximated throughout the thermal layer by its 
value at the wall (which is also the approximation 
which arises in the solution of the flow problem), and 
that the leading term of the potential flow portion of ; 
is taken into account. 

Eq. (2.7) together with boundary conditions on 7 
constitute the heat-transfer problem in the boundary 
layer, the velocity components being assumed known 
The problem is prohibitively difficult since the coeffi 
cients of the convection terms as well as the ronhomo 
geneous dissipation term are complicated functions 
To simplify the equation we now make use of the fact 
that the viscous layer is very thin compared to the 
vanishes and 


thermal layer, and that approaches 


a limiting value outside the viscous layer. As usual 
because of the linearity of the energy equation, the 
solution to the cooling problem with dissipation taken 
into account may be considered to be the sum of two 
solutions, one being the solution to a cooling problem 
with dissipation neglected and the other being the solu- 


at 


tion to a problem in which there is no heat transfer 
the wall. It 
two problems separately, and this will be done in the 


is therefore suflicient to consider those 


following to permit a clearer exposition. 


(2a) The Cooling Problem Without Dissipation 


The velocity components being simple functions of ¢ 
outside the viscous layer, but complicated functions 
within this layer, it is convenient to separate the solu 
Outside the viscous layer 
To obtain the 


tion into two components. 
let the solution be represented by 7 
equation for 7; we make use of the fact that 7,(x, 
vanishes and @,(x, 9) is essentially equal to @,(x, 


outside the viscous layer. Hence (omitting the dis 


sipation term) we may write Eq. (2.7) as 


u” (OT) Ox) + [—Eu 
a “8.(x, ©)] (OT;,/0E 0°7T% Oé 2.8 


Let the complete solution be represented as 


1 T(x, & T(x, 9 9 
so that 7’, is the correction to 7) necessary in the vis 
cous layer. Using Eqs. (2.5), (2.7), (2.8), and (2.9), 
the equation for 7, is 
(1 Tu OT ,/Ox) + ti(OTo/Ox 

—nu’ + 8.) (OT./0n) + ao *[8-(x, n) — 
v(x, ©)] (O7,/0€ o 0°7'./ On 2.10 


The boundary conditions on 7° are that 7 assumes a 
prescribed value 7),(x) on the wall and that the tem 
perature approaches that of the free stream as the 
normal distance from the wall, &, increases. By as- 
sumption, 7, vanishes as the viscous layer coordinate 


nincreases. Hence, 
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HEAT TRANSFER IN LAMINAR 

7 0 Ty(x, 0) + T(x, 0 Tula) % 21] 

] = T(x, o = 0, TAx, ©) = OF ans 

Since o' * is very small, we expand 7) and 7, in a 
series of powers of o ~, expecting that the first one or 
two terms will constitute a good approximation for 
sufficiently small o. We note that if terms whose 
magnitude is of order a'* relative to the leading terms 
are to be meaningful while terms of order (Ro 
have been neglected, we must have 1 R <a’. Let 

es : 2.12 
1 TT + o 7 t { 


substitute in Eq. (2.8) and equate coefficients of like 


powers of @ This gives 


yi Ox) — Eu'’ & 
oT; OF) = (oT, Of 2.te 
u or Ov) — Eu’ (OT; Of) + 3.(x, 4 
oT Of oT es n L.2 2.14 


may be expanded in powers of &. 


In: Bq: (2.10), 7: 
If it is assumed that the series converges in some in 
0, &,), it will converge throughout the viscous 


©1 


terval 


layer when o approaches zero, since the value of 7 cor 
responding to & is m og “é. Using & o' "yn in 
the series, we obtain the explicit dependence of 7») on 
sand 7 in the viscous layer. 

The expansions (2.12) are also introduced. Equat- 


ing coefficients of like powers of o!/*, we obtain 


O7 or _. Ol < 
(Ro + R Ro a 
0 é 0 On 0 


It is worth notice that the functional dependence of the 
result on o has been determined without actual solution 
of the equations, and that the leading term is the ther 
mal layer solution for inviscid flow. 

A general solution to the order indicated will be 


given in Section (3). 


2b) The Insulated Wall, or Thermometer Problem 


The procedure here differs only in detail from that 
followed in Section (2a). The differences are that the 
dissipation term is included in the equation for 7, 
where AT is set equal to l C,, and that the wall con- 
dition is 

oT 
= (0) (2.18) 


It is found that 7) vanishes, so that 7 is of order 


OWs 175 


BOUNDARY-LAYER FI 
"ee oy? = 0, OFZ On 0 
0°71 On- u® + 7,.)(O7 Ox) 4 9.15 


(O7 ; On T @ ol Ox é ,| 


nus?!’ + 


and more complicated expressions for O 1 On*,n >? 


Since 7 vanishes as 7 approaches infinity, so does 
rey On. These conditions are sufficient to allow 
direct integration for 7." when 7 and 7 are 
known. We see immediately that 7 and 7 are 
identically zero. 

Once we know the value of 7 at the wall, the con 
ditions on 7 are given, using Eq. (2.1] 

1 vr, O T(x / v, O TA) 

T rv, O 7 v, O 0 | 

2.16 
7 yr, O 7 r, O iW | 
7 \ 0) NA () 


Phe equations for 7 and 7 may be solved in se 


quence. At each stage the nonhomogeneous terms ot 
the equations and the boundary conditions are com 


pletely determined in terms of the solutions already 


found. 

In the following we shall be concerned only with 
finding the first two terms of tne solution, and there 
fore do not need to be concerned with any 7 equa 


tions beyond those given in Eq. (2.15). The equations 
for Ty) and 7) are (2.13) and (2.14), and the bound 
ary conditions are given in Eq. (2.16). The required 
result is found from Eqs. (2.9), (2.12), and (2.15 
with & Ro)''*y and 7 elo 

oT, (~ 

On 0 Of 0 

oT, ) | [ o7 
i: + Ole Ro | : 

On |,=0 Of ( 


o'/?. It is thus convenient to introduce new variables 


70, 7 


by the relations 


219 


i 0 T 1 es. 


so that 7) and 7, are of order one 

The details of the procedure will be omitted; only 
the final formulation of the problem and the result will 
be given. 


We assume expansions 


™j™ = T T 6 T T { 990 

T= 9 + o's 1 — 
The equations for 79°° and 7‘" are the same as those 
for 7) and 7, in the cooling problem—namely, 
Eqs. (2.13) and (2.14). The boundary conditions on 
ro‘”) are 
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Or Or > [ 
Y T . ) Zz 
O& 0 Ox Jo 
, Or Or \ 
and, in general . a Y. 0 
O& 0 On 


T(x. 0 oI r(x, 0) + 7.(x, 0 


As in the cooling problem, the nature of the dependence 











of the result on the parameter o is obtained without © Iq : [ ' 
ms ; ; L | z, y c ( p 
solving the equation. The leading term is the solution bP 
for inviscid flow, the dissipation appearing as a heat me 
ir ; fo transform the equation we shall require to know 
source distribution on the wall. BG r . Ty hit 4] } 
. : : @r'(0, ¢ O, &). he conditions of the problen 
A general solution for the thermometer problem, for ; 1 i r . i “ I % 
. ’ é' : 3 demand that QO, ¢& QO tor all € except possibly & U) 
arbitrary potential flow u(x), will be given to the rage il - 
refi : oe so that both 7 0, €) and 7 0, £) are zero. At 
indicated order in Section (+). 2 
£ 0, the wall temperature may or may not have 
=e . discontinuity depending on whether or not 77,.(0 () 
3) SOLUTION OF THE COOLING PROBLEM aria eva cee 
We therefore have 
According to our method of approximation, the cool- 
\ cp > - 
ing problem in the absence of dissipation is defined by ¢ $ 
> 2.13), (2.14), (2.16 (2.17 \ solutio , a 
Eqs. 2.13). (2.14), (2.16), and (2.17). A solution fot Che equations and boundary conditions become 
arbitrary u(x) and 7;,(x) will now be determined. 
Eq. (2.13) may be converted to the one-dimensional Pla oe” 0 3.6 
heat equation by introduction of new independent HW 7 
variables. Let 
, where 
2 tV(eldx, § fu 
a L |/ I ( I] ( 
7 \ rab) ( Q, | ) ie 0) ¢ () \ 
- 
] \ Pp } "lay 0) () a] 
df \ 
where . I 
The equations and boundary conditions become 
Phe solution of Eqs. (5.6) and (3.8 
I ? 2 
? p / p Ls 
Inversion by convolution vields 
p(s, O P,(2), PM (s, 0 3.4 ; 
¢ @ 
PY (2 O P(g, 0) 5 d Drll2. 
where subscripts 2, ¢ denote partial differentiation. We are interested in 4 >, 0). If an attempt is made 
Phe desired result can be obtained from Eq. (2.17) by to obtain this function by differentiating Eq. 
noting that and then setting ¢ = O an integral will be obtained which 
diverges as a result of the singularity / it f 0 
rey v, O)/O€ it v) ® gz, U asp _— ‘ — Fis : ; 
rhis is because the process of integration and of letting 
lo find ®”’ we define the Laplace transform with ¢ => O may not be reversed. <A different formula may 
respect to 2: easily be obtained, however, which contains an in- | 
tegral with an integrable singularity. We note that 
Le [b:'(z, 0) ] 9;((s, 0) = —s'/2g(s, 0) = —(1/s'/?) [se,(s)] 


—(] 8 }®,,( +0) + L[@,,’]} 


where we have written ®,,(+-0) to avoid an error in case ®, has a discontinuity at z = 0. Inverting by convolution, 
we have 





lade 


on, 


it Del 
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HEAT TI 


ng understood that the term ® 


Phe solution of Eqs. (3.7) and (5.9) for g'!1 


2 ae sf ie 2 Th 
a ¢ 4 é 
. 25 . 2 
To obtain ®, c, 0) we note that 
ra ‘gy I] 4 ( >. 14 
Now 
i 0 Ee 5 (0) 
x eg | Q i ; Uy 
If we assume that it is permissible to place the inversion 
perator inside the integral, then, inverting by con 
volution and recalling the definition of /7 in Eq. (3.7 
we obtain 
I 0 h(z I 2 \) XX 
\ IV \ 
We now substitute for ®& 2 ¢) using Eq. (3.11)- 
Then, again reversing the order of integration so as to 


integrate first over \ from zero to infinity, and finally 


rung the integral over ¢ from zero to 2 w by 


tranct< 


integration bv parts, we find 


| () 
/ 
LJ, (4) 
hie p,’(< ; 
) ») 
where it is again understood that the first term will 
nish whenever ®,,( +0 0. 


he functions ®,. and /; can be determined from their 


for any given flow and prescribed wall 


definitions (3.1 
temperature 7° Che integrals in Eqs. (3.12), (3.15 
be evaluated in 
numerically. The 
depends only on the flow, not on the wall temperature 
but will generally have to be evaluated numerically. 
s, O) and ®, >, 0) are then used 


; 


closed form, if possible, o1 


can then Ul 
integral appearing in Eq. (2.1% 


rhe results for ®, 
in conjunction with the relation 


oT v, O)/0€ ? 2, O 


and the definition (3.1) of 2 the desired 
wall temperature gradient from the fori-ula (2.17 


to determine 
3a) The Cooling Problem for u©) = x” and T,, = x6 


As an example of the application of the general solu 
tion obtained above, we compute the wall temperature 


flows having the potential velocity 


gradient for wedge 
component 
) {  * 1 


it the wall, and a wall temperature given by 


] 1 Tt A7 | () 


Hence 


Observe that there is no characteristic length this 
problem, but we may, for convenience, retain the 
analysis in the previous form, with / From Eq 
aus 
g } D,.(2 
1 ri 

Che zero-order solution may be computed imn ediately 
from Eq. (3.12), the integral representing Bet ul 
tion. Expressing the latter in terms Gal ul 
tions, we find, when 3 0, and hence, when 4 
() 
Pp () 

I r'; } 

~\ 

™ rs » 2 
When ) \) q Val ishes ind the integt il 1! Ee 
vanishes, while ® 0 f1IVIng 

] () Tr { 4 

As 1s to be expects the result LS roaches 
when 3 — O, the result (3.19) for 0 lo the zero 


order ot « 


O7 0 ; 
OV 


Specializing to the flat plate with constant 
perature, we set 777 OQ and 3 0 
OT (x, 0 Ro 
1 i Q (3.21 
Ov \ 71 
or O7 (x, 0) Ov Uo vrx, 


To determine the first-order approximation we need 
information about the velocity distribution within the 
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boundary layer. The boundary-layer solution for the From Eq. (3.22), we find for y > 
velocity components may be written in terms of the : 
ga. é : D(x, ©) = [alm + 1)/2] x%-Y/? 3.24 
similarity variable + 
, mer and, hence, 
uu + u, = xf! (y; m) 
h(z) = d(x, ~)/u (x) = alm + 1)/?2/22'/? (3.95 
vy —yu + R-O2DG, . ae 
(soe) 
(m 2 , on ‘he j OTs i : 3.15 7 » ¢ ute 
~(1/2)R-OlD~e™ — Ym + 1)f(y; m) 4 lhe integrals in Eq. (3.15) must now be computed, 
When 8 = O the second integral vanishes since ®,, is a 
(m — 1) vf (vy; m constant, and ®,,(+0 1. This gives 
where @-(2, 0 a(m + 1)'/?/qr2'?, B= 0 (3.26 
a >1/24.(m—1)/2, = .(m—1)/$ : 
y= Rx y, OF ¥ A n When 6 + 0, ®,(+0) = O and only the second term 


a ; . in Eq. (3.15) remains. It has been possible to evaluate 
and f(y; mz) satisfies the equation and boundary condi 


ti this integral in the special case when 28/(m + 1 n. 
ions Sep) 5S é 
where 7 is a positive integer. The calculation must be 
m(f’)? — [(m + 1)/2] ff” = m + f’” omitted to conserve space. The result is 
f(0) = f’(0) = 0, f’(~) = 1 2s, 9) = ae > Tern’ x 
ny. ( T?[1 + (”/2)] n | or 
The asymptotic behavior of f for large y is f~ y — a \T2[(1/2) + (n/2)] aa il ~ (9.27 
where a depends only on m. 
Now z,;:= x" (f’ — 1]j-so'that We have not succeeded in establishing that this relation 
, ’ also holds for noninteger values of m—1.e., for arbi 
u(x, n)dn yam ‘| uly trary 26 (mm + 1). It may be noted, however, that 
70 ss if we formally let 6 ~ 0—.e., n —~ O—in Eq. (3.27) the 
— y 523 or : , . = ae 
7 —- limiting result agrees with that given by Eq. (3.26) for 
3 0. The wall temperature gradient now becomes 
OT (x, 0) ” blag - a n/2)] 
—(Ro)! (m + 1) : 
O74 \ r[(1/2) + (n/2)] 
{1 + (n/2)] n a 
ats as(m + ] _— + O(a*'-) 3.28 
| r((1/2) + (m/2)] 2 J 


For the special case of the flat plate with constant wall temperature we have 
OT (x, 0)/Oy = —[(Ro)"2/(ax) 2] [1 — 62a’ er! 2)] + Oo 3.29 
Inserting the value 1.721 for a (see reference 3, p. 136), we have 
—(x/R)'? [OT (x, 0)/Ov] = 0.5640!/2 — 0.5470 + 0(0*”). (3.30 


The right-hand side of Eq. (3.30) is plotted as a func- 
tion of g in Fig. 1. The results computed numerically by E. Pohlhausen* are shown, and E. M. Sparrow’s ap- 
proximation® based on the Karman-Pohlhausen method 
is also given. The present results are in general agree- 
ment with Sparrow's at very low o, with Sparrow's 
values slightly smaller. Since Sparrow's values for 
ae larger o are also slightly smaller than Pohlhausen’s, we 


: a 
—T may expect that the exact curve should lie higher than 
iB Sparrow's throughout the interval. In this special case, 
he a the available comparison with Pohlhausen’s and Spar- 
a a Se 2 row's computations indicates that our zero-order ap- 
n 1 ‘mee | proximate solution is quite good up to, say, ¢ = 0.03 
— 7 a 
RE ca , 1 and the first-order solution up to, say, ¢ = 0.1. Since 
os no comparisons appear to be available for other cases, 
pa ae eo ‘ ‘ . 4 . 
as | | ia 5 the above limits might be taken as a first guess for the 
5 i ° o ge ° . 
f ace | range of validity of our general solution. 
- a : : 7 : It may be of interest to note that the results derived 
Fic. 1. Wall temperature gradient vs. Prandtl Number in in this Section from the general solution may also be 
forced convection over a flat plate at constant temperature 7\,. : é 4 j‘ : 
T = (T; — Tow)/(To — Ta), where 7; is the temperature and obtained by looking for solutions to the energy equa- 
T.. the temperature of the free stream; A = first term of series tion in terms of the similarity variable 7 i.e. by set- 


solutions, B = first two terms of series solution, C = Sparrow’s : ‘ 8 : 3 : é 
result, © = Pohlhausen’s exact results ting 7) = x’x(y) to obtain an ordinary differential 
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equation for x. The boundary conditions to be applied A 
are the same as those for 7), and the solution when 28 i 
; : - ia ; ‘ > : 7) 
m + is a nonnegative integer is found in terms ol & 
. ° _— ° o 
Hermite functions. The temperature gradient may ys | T 
; ; es ~ ‘ | | | | — — -. 
then be determined from Eq. (2.17). The results de — ys 
rived by this method agree with those obtained by Pe T 
; 3 + + 3s -F 
means of the general solution. ZA 
ie 
we 
{) SOLUTION OF THE INSULATED WALL PROBLEM a 
Phe differential equations governing this problem are f. 
the same as those for the cooling problem, viz., Eqs. 
913) and (2.14). The boundary conditions are 3 - 
2.21), and the principal result, the recovery temper Fic. 2. Recovery factor 7 ] Tad Prandt 
:ture, is to be determined from Eq. (2.22). The method Rramber sn forced convection over & pate thermometer. @ 
: ; is the wall temperature, 7. the free-stream temperature; A 
f solution 1s completely analogous to that employed first term of series solution, B first two terms of seri lu 
- . ° nr a smooth « ve hroncl ohlhaus ’ esults id low-e 
in the cooling problem. We introduce the new inde tion, smooth curve through Pohlhausen’s results at 
; : p ings results, and Pohlhausen’s exact results 
pendent variables (z, ¢) by the transformation (3.1) and 
defin« 
@'(z, ¢ r ,¢ j 0, | those in Section (3)—that is, Eqs. (3.2) and (3.5 Phe 
> boundary conditions are 
his g n 
{ J) 1.1 ®, e () D Pp ~ a) 1.2 
\ 


Ol ad ~ ~\d - ) ad > . 
D(z ( ) dy | | 5, O [(z)4 z, | | z, 0) Ea 
i v) JO On 


Again solving by using Laplace transforms, we find 


| % D(z tt 
ri - 0) 1.4 


he equations to be satisfied by ®’' are the same as 


oe ut wf | h(w)D(t &,(t, O 
P - () dtdwad al to) 
2rJo0 Jo JO (A t g \ Varo Ve / 
4a) The Insulated Wall Problem for uu’ = x 
As an example we apply the general solution to the case vy”, the wedge flow. In terms of the functions de 
fined in Section (3a) and Eq. (4.1), we have 
h(z 0(X d \ avm+12y 
2 udn - a |(m+1)2] 
“tt i 6 
l © (dit, , ; 
D(z - dn [(m + ha pis f"(y; m) }*dy 
i Y) #0 On J 0 
Substitution in Eqs. (4.4) and (4.5) and integration vields 
lr} (2m/(m + 1)] + (1/2); F , > 
f(z, 0 m-—+sl)s\ if (y; m ay 4.7 
Vim 11} [2m/(m + 1)J +1; Yo 
PY(s O al(m + 1)s] [f"(y: m) ]*dy7 1.8 
e U0 
The recovery temperature is given by Eq. (2.22), using Eq. (3.1 and writing B 2m 'm + 1 for the wedge half- 
angle in right angles 
T(x, 0 rmB+(/2)1 ey | Pe 
= g’'* LJ y;m |? Gy T Ci @ i} y; m dy 
v VUm-+i1Tr(B+ 1)J0 J0 
vLf"(y; m ray | t+ O(a*!? 1.9 
“i 
The results derived here have been checked at B = 0, 1/2, 1 by looking for solutions of the energy equation of the 


Ir 


form Ty = x°"x(4 


In the special case of the flat plate, B = m = 0, and 
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o fe ' od Lf’(y; 0 "dy 


“) 


a Lol ee Y 0) dr 
JV 
we have bie 0.4620 
Phe right-hand side of this equation is plotted in Fig. 2 
and the numerical results by Pohlhausen* for 0.6 


o < 1.0 are also shown. The smooth curve through 
Pohlhausen’s values and the present low-o approxi 
mation is expected to be an accurate representation of 
the recovery temperature for the whole interval (0, | 
The apparent range of validity of our zero-order ap 
proximation is roughly o 0.15, while the first-order 
solution is good in the range o@ < 0.50. As in the cool 
ing problem, we take these limits of validity for the 
flat plate as a first guess at the limits of validity for the 
general solution. 

The approximation here appears to be valid for much 
larger o than the corresponding approximation in the 
cooling 


s 


lowing 


This may be explained by the fol- 
the 
greater, the smaller the ratio of viscous layer thick 


problein. 


reasoning: accuracy of our method is 


ness to thermal layer thickness. In the cooling prob 
lem, heat enters the fluid only from the wall, producing 
a thermal layer of some thickness. In the thermom 
eter problem, heat is generated as far away from the 
wall as the outer edge of the viscous layer, hence the 
thermal layer is thicker than in the cooling problem 
The 
thermometer 


by roughly the thickness of the viscous layer. 
thus the 
and the approximation is more 


thickness ratio is smaller in 
problem (given R, o 


accurate. 


derivation of these corrections on the basis of thin air 
foil theory was, of course, inadequate since the flow is 
clearly nonlinear in local regions, but when the deriva 
tion is based on the extended d'Alembert paradox the 
corrections are seen to apply to more realistic fluids 
provided only that the flow remains of the reversible 


subsonic type. 
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7, 0) |°d> vif" (y; O | + Ole 
ei) =) 
vy: ( Us g 
0.2607, v[f"(y; 0) ] 0.3515 
ei) 
0.097¢ + Ola 9 


5) CONCLUSION 


We have presented a method of solution of the energ 
equation for laminar forced convection problems with 
Prandtl Number. 
Prandtl Number were obtained, 


tions were applied to give specific results for flow over 


small General solutions for very 


low and these solu 


wedges. The results for the flat plate were evaluated 
numerically and were found to be consistent with exist- 
ing results, the comparison yielding an indication of th 


range of validity of our solution. 
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SUMMARY 


ne stress problem for subjected to an arbitrary 


trical 


i plate, 


temperature distribution, is solved for the case 


he edges or the issumed to remain 


buckled 


terms of 


plate are straight 


flection function for the plate is then approxi 


the first four double infinite harmoni 


used to determine 


Rayleigh-Ritz procedure is 
when initial buckling occurs. 
SYMBOLS 


effective areas of supporting structure in the 





nd J directions, respectively 
2 length of plate 
2 width of plate 
D flexural modulus of rigidity 
D Ee /12(1 — pv? 
iverage modulus of elasticity of pl ite material 
Ky, A extensional rigidity parameters in the x and 
b directions, respectivel K ] T 1/@ 
K 1 + 1/9, 
thermal buckling coefficient 
\ tie stress resultants 
difference between the average temperature of 
the supporting structure and room tempera- 
ture 
difference between the temperature at the edge 
of the plate and the average temperature of 
the supporting structure 
difference between the temperature at the 
center of the plate ind the plate edges 
thickness of plat 
coefficient of linear expansion at average 
temperature of the pl it¢ 
Poisson’s ratio 
b,, 4 ratios of extensional rigidities of effective area 
of supporting structure to the skin in the 
vr and y directions, respectively (4 
1-1 /btE, ®, 1, /atk 
Bars over letters refer to supporting structur 

















INTRODUCTION 


— WINGs of an airplane flying at supersonic 
speeds are aerodynamically heated to high tem 
peratures. In the initial stages of heating the tempera 
of the skin is nonuniform and higher than the 
the 


Buckling of the skin may follow since, for example, 


ture 


temperature of supporting interior structures. 
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the cooler interior structure will prevent the skin from 


expanding freely. 

It is assumed that the supports remain straight, do 
not deflect transversely, and offer no rotational stiffness 
In general, the skin and supporting structure may be 
the result that 


composed of different materials, with 


thermal stresses will be induced when the temperature 


is raised uniformly. The values of the physical con 
the 


effective 


stants are taken to be equal to those existing at 
the 


1 


average temperature of skin and of the 


supporting structure, respectively. 


THEORY 


The equilibrium conditions of a plane stress problem 
are satisfied if an Airy stress function, ¥, exists such 
that the stresses are defined as follows 

oy O*y 
oy" O. 


O°y 
Ovdy 


For a simply-connected region the compatibility 


expression! for the thermoelastic plane stress problem is 


The temperature distribution 7 acting on the plate 
is assumed to be 


of the plate 


symmetrical about the centerline 


Fig. 1) and can be expressed as 


] T T ] T COS ad COS T 


. S's 


0 0 


where 7 is the difference between the average tempera 
ture of the supporting structure and room temperature 
and 7 the 


the edge of the plate and the average temperature of 


is the difference between temperature at 


the supporting structure. 
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The stress function 
y (Kiy?/2) + (Kox? 2) + 35 tw cos (sax /a) 4 

> fo, cos (try h) + > > t,, cos (sm/a)x X 
t=1 1 t=] 
cos (fr b)y (4) 
will be shown to satisfy the following conditions: 
0 K, _ 24 tr b) cos (fmv bd) 
t=] 
o Ks — bo (sw a)- cos (s7x @) 
l 
T — SY t(sx/a)(tr, b) sin 
1 ¢ ] 
where 
t Kea (S77 a) 1 
to, [Ka/ (tw/b)?|Te, (6) 
bse EaT ,,/|(sx/a)? + (tw/b)?] 
The linear thermal stress-strain relationship is 
€ l E)lo, — vo tal 
(7) 
€ | Ek 0 ae PCs) al 

Integration of the stress-strain relation results in the 
following normal displacement at the edge x 2a: 

ii (Qa/k)(A, — vAe + Ea(T + To + Tw) S) 
where 7°, 7 when s t 0 and 0 is im- 
posed. Similarly, 

(2b/E) [Ke — vA, Eal(lT + To + To (9) 


The resultant forces on any cross section of the plate 
parallel to the x and y axes are, respectively, 


Vg t | a,dy 


2K bt 


10) 


Ps t : o,dx 


2 Kyat 


Since there are no external loads acting on the system 
(plate and supports), the resultant forces in the plate 
the 


Thus the loads acting on the 


must be equilibrated by the forces acting on 
supporting members. 
stringers in the v and y directions are 


P, —K,bt, P, — Koat (11) 


From Eqs. (8) and (9) it is clear that the displace- 
ment of the plate edges depend only upon the constant 
stresses, A, and A, and upon the average temperature 
of the plate, T + 7) + Tw. The resultant force on 
any cross section is due to A; and Kz only. Therefore, 
the nonconstant temperature terms produce stresses 
which are self-equilibrating. The strain in the plate 
consists of two parts: (1) a constant strain, due to 
the average temperature on the plate; and (2) a strain 


which varies over the plate, but with the restriction 


LACAL SCIENCES 


(a) compatibility, [see Eq. 


MARCH 


(9)] 


) I, 


1958 


and (b) boundary, in 


which (1) the edges remain straight, and (2) the average 


strain along each edge equals the strain in the support— 


1.e. 


Substituting Eqs. (43) and 


(4) 


into Eqs. (2) and 


results in the following expressions for the stresses 


= »_ t.(tm/0)" cos (sax 
1 t l 
=: Sy 
me *" ( $7 /a)- COS (S#X/a) COS 


(sm@x,/a) sin (tary b 


that the resultant displacements « 


plate are zero. 


The 


boundary condition 


for 


the 


f the edges of 


strain compatibility 


along each edge of the plate is approximated by equat- 


ing the strain in the supporting members to the average 


strain of the plate edges. 


Substituting Eqs. (3), (5), 
yields 
[Ea(T) + Tw) + E 
K, 
: [Ha(7 1 / 
Ky 





‘hus 


€,a \ Z 
20 J0 
nd (11) into Eq. (12 
a)T \(v Pa 4 
- aT \(p ee 
14 








---- APPROXIMATE SOLUTION ( REF 2) 
— SOLUTION OF EQUATION (22) 
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BUCKLING OF SIMPLY SUPPORTED 
where K, = 1+ (1/®,), A, l + (1/9,) 15) A, is equal to the stringer area for a single panel and 
} LE bE. © 1 FE /atk (16) is equal to 1/2 the area for an interior panel 
Substituting Eqs. (13), (14), and (6) into Eq. (5 
results in the following stresses: 
Ea(T + 7 + Kila — @& T i fs =" . try 
7 r < Ea ia 7, cos . = 
v> — KK t=1 h 
~ ~ 1 TA trv 
Ea > + . : COs cos : 
1 1 |(s/t)*(0/a)* + | a D 
Ea(l T 1 \+ kia — a T L K = S7X Re 
© a ke Law 7 COS ge ‘ 
1 KA a 
. = D's: = LP, 
ka haut Lal ; cos COs 
l ] + t \-(a fp) a ) 
~ ~ ‘f Sx lr 
; -ka +2 _ : —~ ant sin 
lf s [)(0b,a) + (ft s)(a 0) a h 
BUCKLING OF PLATE 
The deflection function is approximated by four nee 
; Shanes gage 2 ; munimize the load necessary to cause the deformation 
terms of an infinite series in the form : pls 
(1S) for each aspect ratio a/b, and each temperature 
> Dd an, sin (27/Ap)x sin (27/u,)y (18) ratio 7) 7). 
l 
where Rayleigh-Ritz Procedure 
\ la'mu, 4 ih /es 19) The total change of energy of the system at buckling 
is due to (1) the strain energy of bending of the plate, 


u, v are arbitrary positive integers. 

[he deflection mode in the x direction, for example, 
will be symmetrical or antisymmetrical depending upon 
whether u is odd or even. The problem then is to 
determine the parameters “, v in such a way as to 


; D [ [ | (2= oe) or 
2 J0 JO Orv Ov fi 





jow Ow 
(Ox? Oy 


and (2) the total work done by the forces acting in 
the middle plane of the plate during buckling. The 
assumption is made that the stresses acting in the 
middle plane of the plate remain constant during the 
deflection, and thus the stresses in the supports must 
also remain constant. Thus the total change of energy 


( O7w yy | 
- dydx + 
OxOy f ‘ 


18 


I — . " Ow \* . Ww \- _ OW OU 
] ] 7 ( , " j ( : ” - dyes . 
2 Jo JO OX '\ Oy Ox Oy J ~ 


Since the total change of energy at buckling must be stationary, then, 


OU /Oa», 0 


Substituting Eqs. (17) and (18) into Eqs. 


equations. 


Ky K;; 
Li Ly 
My Mhs 
Ni N 


where 


l ( b ) fa l l b 
K 11 u- + v- ( : + ul 6, 
192 \a F D k, t 7r° a 


K L - oe 4 
. bor? [u?(b/a) + 4v7(a/b) 


l u7v"t ou l u-v~ 


Ky, My 


to? [4u?(b a) + v*(a bd) 


N 


m,n 1, 3) 21) 


(20) and (21) leads to four linear homogeneous simultaneous algebraic 


The nontrivial solution of these equations leads to the following determinantal equation: 


Ks, Kz 
me bi . 
0 22 
Ms; Ms 
N33 
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; 3 uv" tio, 2] 
Kz Nu . i: 

Sr? [4u?(b/a) + v?(a/b) Sa? [> 


l b 
it lj 
1 ar* a 
5 u-v-t (2 l u-v- 
Zz VW, ; 
Sar? [4u7(b a) + v*(a bd) nw” |u-(b/a 
; 9) u-v-t 8) 
i N 13 : + ‘ T ‘ < 
bor? [4u7(b/a) + 9v7(a/b) mr” [u 
l h \ [ a | l {) 
Ws) | Qu Vv" + u 
192 ad L b 4 k, lr 
9 u-v-t 3s #8 
V/ \ 
tr? [9u7(b/ a) ly?(a/b tr? a 


where 
j/ 1 [ a tin 
6 | — 
' 7 | a} T7,! 
vy + K;) 
v, 4 2 
| (v? — K,A 
and 
: . . 
} l y} x 
a6 JO J70 
mutT Nv : 7 
co v COs vy dyvdx pas 
a ) 
h V h t) al; R 26) 
CONCLUSIONS 
The determinantal equation (22) is easily solved for 


the unknown k, where w and v may be any positive 


integers. The smallest value for the temperature is the 
critical value. The critical value for the temperature 
corresponds to i v 1 in all cases for which the 
stiffnesses of the supporting structure are the same in 
the vx and y directions. This result is consistent with 
the exact solution for the case of a plate under a uni- 
form temperature (see reference 2). 

For the purposes of determining the effect of the 


nonuniformity of temperature, calculations were per- 
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u-v-t 
b/a) + 4v72(a/b) 
g a l [ b fu a [i 
} Z 6, — =e + \] | 
bar D 27 ri La z b\u/ | 
vt 2 5 u-v-l 
) v*(a/b) Sa? [u*(b/a) tv?(a b 





formed for the case of a parabolic temperature varia 


tion on the plate —1.e 
ii Fa T,;1 - a i431 
where ‘i 0, A K 


The results are presented in Fig. 2. Comparing the 
Z. that accurate 
the plat 


a uniform temperature equal to the 


two sets of curves of Fig. we observe 
results can be obtained if it is assumed that 
is subjected to 
average temperature of the plate, only for 79/7; 2 
For 7/7; - 


distribution is 


», the effect of the self equilibrating stress 
quite important and the results givet 


in this note are significant. It should be noted that 


as the stiffnesses of the supports decrease, the greatet 


will be the discrepancy between the two solutions 


Of course, in the most severe case where there aré 


supports (A K ) then the use of the averagt 


value has no significance. 
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SUMMARY 





e approximate method is advanced for the determ1 
upersonmic mVvi cid xisymmetric flow field ver 
lid body of revolution and over its wake his 
equal ipplicable to flow fields with diverging stream 
ines in the front end of the body and to flow fields with converg 
g iline t the rear end of the body ( er the wake 
nd gives good results for problems wit values of 
ic parameter Ve t larger than O(1 rhe effect 
I ul lit 1 includes ecurac I thi } i t¢ 
t is tested I give cvlinders 
( te calculation, the t get ¢ ppl 1 
pted in the nose region qd continued it the present 
Phi rhac res re coelhcient ver the ft t ( I 
ve as calculated tf 1) ‘‘ecomplete characteristic caleul 
t 2)t present approximation 1¢ }) other appt 
in Fig. 3 for comparison. The present result agree 
rly we ith t ( plete characteristic uti It 
ist dt t pres¢ ( iit i ensitive t 
tain factors in the approximate cal ti I time 
for culatin pressure curve s in Fig 1 
1 t e hours at t 
i iplete « cteristic s 1 lil the g1 
\ i ilable for pari either 1 
] \ l te scheme | ) I 1 gi 
| t Is thu t teste 1 it cc ( uci 
I" 1 1 ic foundation, the present t lis « 
( 1 in flow gions with either « erging 
ging str line Di t compari if is cer 
esirabl 
g ze 1 Savin‘ ed that 1 ISV1 et pel i 
the local variation is essentially plar rhe resent 
indicates that is the case witl pe i 
I et ¢ ] ut for the wer hyper rT ge { 
ay ic pat eter ¢ > 1, the deviatu ] the 
f 
| inplicity of the present app i cheme is preset 
i shapes involving discontinuous slopes 


INTRODUCTION 


, | {We METHOD of characteristics for the determination 
of the inviscid supersonic flow field over a solid body 


has been extensively developed and is well known. It is 


usually carried out stepwise partly graphically and 
partly numerically. It is capable of providing an ac 
curate numerical solution of a specific problem under 


usual circumstances and is often considered as the 
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however, tedious and 


Wi 
the de 


“exact” solution. The method is 


not convenient for many purposes. would like to 
flow 


the 


have a fast method of calculating ired 


property in a restricted region, appropriate for 
specific purpose, even if it is less accurat 
he small perturbation theory with or without higher 
order corrections will often serve this purpose. This is 
especially true for the planar supersonic flow over a thin 
The most convement aspect of the 
is that the 
function only of the local 


Higher 


forms are comparative ly 


airfoil. ipproximate 


result static pressure at a given point 1s a 
streamline inclination ¢ up to 
corrections in different 


the order of ¢ order 


simple. For the planar flow 
problems, the differential form of the Prandtl-Meyer re 


the local 
The 


lation can be extensively used to evaluate 


changes of flow properties error. 


situation of supersonic flow over a slender body of 
revolution is unfortunately considerably more comph 
cated. The static pressure at a given poimt on the sur 


face of a smooth pointed body is no longer determined 
by the local 
but is dependent extensively 1] 


Neither is the local change of the stati pre 


treamline first order 
geometry. 
ure simply 
determined by the local change ol unline cde 
flection. 

If there is a sudden change of the slope of the meridian 
curve of the body, the linearized solution of the flow 
field the slope 
quite complicated. Higher order appr 


Che dependences of the flow field at a 


downstream of discontinuity becomes 
Ximations are 
even more so. 
given point on the details of the upstream geometry 1s an 
This 


an isolated portion of 


to the de 


essential feature of the axisymmetric problem. 
aspect makes it difficult to study 
the downstream flow field without referring 
tails of the upstream conditions 

There are several other a methods relat 


pproximate 


ing the local static pressure and the local 


treamline de 
xisymmetrie prob 
the 
those of 


shock ex 
these 


flection which may be applied to 


lems under certain circumstances—for example 


various forms of Newtonian approximation 
the tangent-cone approximations, and the 
The 
methods has been evaluated.’ 
methods take into account the effect of the upstream 
A weak 
feature common to all these approximate methods 1s 


There 


pansion methods. relative accuracy of 


The shock expansion 
geometry, while the others are purely local. 


that they are applicable only in the nose region. 
is, as yet, no approximate method available for the de- 
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termination of the inviscid supersonic flow field with 
converging streamlines downstream of the maximum 
cross section of the body. The shock expansion method 
will give an answer, but its accuracy is doubtful. If 
one is not willing to use the method of characteristics 
and or the approximate method of Lighthill, etc., one 
must develop some new approximate method for such 
problems as the determination of the static pressure 
over the effective wake boundary of the body. It can- 
not be expected that the approximate method can do 
away with the overall effect of the upstream geometry 
in general, like the tangent-cone approximation for the 
flow in nose region. It is hoped, however, that some 
overall parameter may be introduced to represent the 
integrated effect. To have some feeling for such a 
parameter, the problem with nonuniform inflow over a 
slender body was first investigated with the method of 
singularities. With this in mind, the method of charac- 
teristics is linearized and then developed into a simple 
treatment of the supersonic flow field over axisymmetric 
body and wake. For details of the development and 


the sample computations, see reference 9. 


(2) THE APPROXIMATE METHOD—-FORMULATION 
AND DISCUSSION 


Let us consider the elementary characteristic net 
1—3-—2 in a meridian section as shown in Fig. 1 where x is 
the distance along the axis of symmetry and 1 is the 
radial distance from the axis. 1-3 is the positive family 
characteristics or Mach line, and 3-2 is the negative 
family characteristics. gq is the local flow speed and 
¢ is the angle between the local velocity vector and the 
axis, measured positive if the streamline is directed 
away from the axis (negative as shown in Fig. 1). De- 
note the Mach angle by u, the flow Mach Number by ./ 
with 6? = \/* — 1 = cot? uw. Also let v be the angle 
through which a sonic stream will have to turn in a 
planar isentropic flow to reach the local Mach Number 
\J—1.e., the Prandtl-Meyer angle given as 


v [((y + 1)/(y — 1)]*'? tan—' } [(y — 1) + 

(y + 1) — 1h”? — tan- (IP? — 1)? 
for gases with constant specific heat ratio y. The 
functions v, uw, 6 are all tabulated.® 


Let Ags 3 — ¢, etc., then the characteristic re 
lations valid along 13 and 32, respectively, may by 


written as 


91 ( Agai 41) Ags + [sin by sill nA = 
Cos (sh -F 1) | (Axs) i ie sin KM COS udSa| 
B2( Ades gz) —Ago; + [sin uw sin g - 
COS (wo + ge)| (AXNo3; %) — sin uw cos wAS 


provided that 7 is not zero and that neither of the tw 
Mach lines are normal to the axis. Here Sis the entropy 
of the gas normalized in proper units. 

Let us consider the lecal flow field in which the flow 
deflection angle varies along the positive family Mach 
line according to the power law 


o(r)/eolte) ofa 9 
where s is some constant which will vary throughout 
the flow field. Large values of s mean small variation 
of ¢. This representation of the variation of ¢ is ap 
propriate only when ¢(™) # 0. Assuming Eq. (2) is 


locally valid, we have 
Agsi —(v/s)(Ars, 7) 3 


Through some algebraic manipulations with the charac- 
teristic relations (1) and the relation (3) and some 
geometric identities, assuming locally straight Mach 
lines for a small net 1—5—2, there is obtained 


Agni Ax», | sinusin ¢cos¢ 
Avy; ro) =i Ag Thee: 
q r | COS Lb 
¢ sin (u + g) sin \4t— ¢g) Ara ¢ 
» COS uU sin Lu 7 Y 
sin (u + 2) COS (tn — ©) ; 
— Sin pcos uwS Ay 1 } 


sin MCOS LL 


where Sy is the local entropy gradient across the stream- 
line in the units of yF& for ideal gas, Y being the stream 
function. Quantities without subscripts mean the 
average value over the net. 

We shall investigate the possible simplifications of 
Eq. (4) under varying circumstances. 

For ordinary supersonic flows over an unyawed 
slender body of revolution, the following conditions are 


satisfied: 


g/ul << 1, sng2vo, cosgo=]!1 


The first condition in Eq. (5) is the same as | .\/¢ l, 
which is the well-known basic requirement of the 
linearized supersonic flow analysis. There is some dif- 
ference, however; the .1/y here is the local value, not 
the product of the undisturbed stream Mach Number 
and the semi-nose angle of the body. This local condition 
is usually less restrictive. Under conditions (5), the re 
lation (4) becomes 


Ave, —Agor — [1 —_ (1 S) Jo( Axo, Br) = 
(gy §)( Aro, fj — [3 (1 — 87) |Sy Ayo ()) 


which will enable us to determine v2 if g. is known at 
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provided also that we know the values of s and 


point ~ 
the flow is not 


§, at point | and that 8 is not too small 
) the upper transonic range). 

The physical origin of the four terms on the right- 
the middle two are 


‘ , 
hand 


side of Eq. (6) is quite clear; 


the ixisvinmetrie corrections. It also be noted 
that with S, = Oforisentropic flow, the relation (6) with 


reduces to the result obtained by a source dis 


may 


tribution along the axis with nonuniform inflow. 
For flows with large Mach Numbers ¢ and uw may be 
f comparable magnitude and the local value of the hy 


mic parameter l\Jo becomes of the order of unity 


perst 
r larger Practical situations are usually such that 
oth uw and ¢ are small, so that 
go/p! =~ O(1 
7 
singe=v, snu=u cosg=>cosz~p=>!] 
Then Eq 1) becomes 
_ Ago T l s oa An )? 
v/s)(Ary/r) — |8/(1 + B*)S,Ay 
yg 3 r)(goAn A? S 


is an additive correction to re 
ition (6) due to the hypersonic effect 3¢, > 0 (1). Itis 
that this correction vanishes if the 
Ax. Thus, for 


the evaluation of the pressure on the surface of the solid 


The last term in Eq. (8) 


nteresting to note 
ith line 21 is a streamline with ¢ Ar 
body, which is a streamline itself, Eqs. (6) and (8) be 
come identical and the entropy gradient term also drops 
out as Ayo; 0. While there is little difficulty in carry- 
ng these two terms along if we are interested in the 
flow properties at points not on the same streamline, 
the following development will be concerned mostly 
with the particular cases with the path line 2-1 either 
coinciding with or not deviating significantly from the 


streamline. For such cases, the equation 


[1 — (1/s) ]e(Ax Br) 


o/s Afro Yr) q) 


Ai Ag l - 


may be used throughout the supersonic-hypersonic 
range. 

Before considering these equations any further, we 
have as to the initial magnitude ot 


and its variation throughout the downstream field. 


to get some idea 


is the parameter representing the extent of the 
local axisymmetric correction on the planar contribu- 
tion Aw for the determination of the local change of 
flow speed. 
the axisymmetry and rotationality upstream of the 
given point as reflected in the local variation of stream 
deflection along the positive family Mach line. There- 
fore, 1 s must represent some integrated property of the 
upstream geometry and must vary from point to point. 
The supersonic flow field over a pointed, smooth con 
toured axisymmetric body at zero lift is well known. 
We can —O(In 7) O(In r) along the 
positive family Mach line from the solution with some 
We find that the 


evaluate l/s 


representative geometrical shape. 
linearized slender body theory gives 


ISCID SI 


It also represents the integrated effect of 


PERSONIC FLOW 1837 


| = {1 /s Oller (x — Br Oboe) < | 10 


If we take the Taylor-Maccoll solution over an unyawed 


cone and evaluate at the cone surface where ¢ we 


have 
l/s tan 2/Q2)[1/(1 + 6B tan Q 1] 
which ts valid for values of 2 and 8 not necessarily small 


In the limit of 3 
from slender body theory. 


this result agrees with that 


Calculations show that the 


value of 1 sis usually between 0.5 to 0.7, and decreases 
toward the shock (receding from the solid surface 
Furthermore, it may be shown that for a particular 


type of nonuniform inflow over a blunt section, the 


linearized (isentropic) supersonic flow theory will lead 


to a constant value of 1 1/2. On the other hand 
at the point where a streamline becomes parallel 
to the axis, s must vanish in a definite manner as 1s 
obvious from relation (2) if the flow field is locally 
regular. This last aspect certainly indicates that the 
functional form of (2) is inappropriate under this 
circumstance 

Consider (Fig. 2) two neighboring points 1-2 and 
denote the positive family Mach line through 1 and 
that through 2 by 1-3-5 and 2-4. Let 32 and 54 be the 
negative family Mach lines through 3 and 5, respec- 


tively. Then write down the characteristic relations 


along 24, 32, 35, and 54, with both ¢ and u sufficiently 
small to have their sines given approximately by their 


angles 
Ap t Ag T Oo A A y) l T y) 
Introduce the definitions 


A¢ ¢3/ S3)( Ars 7 
Ae —(wo Ss Ar } 


{ »y 


and eliminate v’s. Then take the limit that 4 ~ 2 and 


After dropping higher order small terms and 


- 9 
yo 2B. 


letting 


h l 2 y 14 
we obtain 
Ox OX), Ox Or ; UL ¢ Pu 
f ‘ t ‘ \ 
wet ¢)/2uley 4 I go)/2\(1/r) > 
) ») ») { = 
} b Yo —M rei a ~My ie) 


~~ +— ZZ, 4r 
1 '>~ ts Ppp “a 
9, 

‘ st 
4 3 

XN 1 

y 
4 

Fic. 2 








188 [TOURNAL OF THE ABRONAUTICAL SCIENCES 


where subscript p indicates differentiation along the 
path 2-1 and_ subscript indicates differentiation 
along a positive family Mach line. From this relation, 
we can infer that the fractional variation of « is of the 
order of uw or of ¢ if both of them are sufficiently small. 


Accordingly, a rough approximation will be taken that 
kK (1 (2/s)]¢ constant (16) 


This relation is not in contradiction with the pre 
vious observation that s 2 represents a stationary 
value of s in a particular linearized isentropic supersonic 
flow field. 

The simple relation (16) certainly cannot be a valid 
approximation in the immediate vicinity of ¢ 0 be 
cause the postulated variation of ¢ in the form of rela 
tion (2) loses its significance. s must vanish as ¢ ap 
proaches zero if the relation (2) implies a nonuniform 
stream. 
is at least qualitatively compatible in requiring s to 
In fact, if we 


In this respect, the rough approximation (16) 


vanish with ¢ if « is a nonzero constant. 
adopt for this region a more reasonable description of 
the » variation along the Mach line, in terms of a 


parameter k, as 
tal Yi T Ri \(r - Te ri} (17) 


we obtain a relation identical with Eq. (15) with « re 
placed by k. This indicates that it may be feasible to 
adopt the limiting definition 


lin [—(g/s)] lim « k (18) 


0 > 0 


and let relation (16) be approximately valid even when 
the streamline locally becomes parallel to the axis of 
symiunetry. 

It cannot be overemphasized here that Eq. (16) must 
be regarded more as an empirical relation selected for 
convenience rather than analytical consequence even if 
[O(In «)/Or|+ is, at its maximum, O(g). In practical 
applications, neither w nor ¢g is sufficiently small to 
warrant the order of magnitude. Consequently, 
reasonable numerical accuracy can be expected only 
when the axisymmetric contributions (involving ¢ s) 
may be considered as corrections to the planar contribu 
tion 

Eq. (9) may be written, for Ago; # 0 and ¢ ¥ O, as 


Avy = —Agu(l — (kpy/Br)} 1 — [1 + Be)/s]}) (19) 


where Ry is the local radius of curvature of the stream 
line through point 1, the second term being the ratio of 
the axisymmetric contribution to the planar contribu 
tion in the determination of Av. Regardless of the 


magnitude of By, the factor g[1 — (1 + Be/s)] is 
O(x), under the restriction (16), which is determined by 
the initial value of x. According to the approximate 
scheme proposed later for a pointed nose body, the 


tangent cone approximation will be used with 
Ky ¢go(Bogo — 1)/(1 + Boge) |< | ¢o 


Thus, kp go, 8¢o(”%/ go) determines the relative importance 


of the axisymmetric correction. It is then apparent 
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that, for a curved convex body where kp¢ 


usually less than unity, the axisymmetric correction js 


relatively unimportant when Sygp becomes larg: Chis 
is i agreement with the result of Eggers and Savin‘ that 
axisymmetric hypersonic flows are locally planar to th 
first approximation. The present result indicates 
however, that py 1 does not seem to be larg 


enough. Later calculations do indicate the importar 


of the axisymmetric correction. 


Over a conical or even cylindrical portion of a body 
of revolution the local planar contributions inish 
identically and the axisymmetric contributions pr 


dominate. Therefore, an approximate description 
the axisymmetric effect in terms of the variation of 


gross parameter is not likely to be a valid approximat 


scheme. For lack of anything better, it is propose 
that the following variation of k in relation (17) may 
serve as a rough approximation: 

k k exp i] —- (2n Or) . > 1/37 t A 
where ko ko 2 evaluated at the station .o where th 


cylinder portion begins. The constant 7 will be selecte; 
that the undisturbed free-stream pressure may be ap 
proximately reached at the far downstream end if th 
cylinder should extend to downstream infinity. In the 
case of a conical portion, the selection of 1 will be based 
on the presumption that in the limit of very large «x, the 
surface pressure must approach the Taylor-Maccoll 
value. 

Having selected the law of variation of the parameter 
and relation (16) or (20) for their respective range of ap 
plication, we may look upon Eq. (9) as a differential 
equation for v(x) if g(x) and 7(xv) are prescribed. 6 isa 
known function of vy. It is considered more consistent 
however, to adopt the following approximation 


where B is taken to be some average value of 3 between 
the initial station vj and that at the station 
For a pointed nose curved body, the following proc 
dure of approximate calculation is proposed 
1) The initial value of «x gol] 2 so)] will be 
evaluated at some convenient upstream station with the 


tangent cone approximation as 


K goll — (2/s gy —2tan ¢ l+ tan ¢ 92 


2) Estimate the value of 3 dim. (¥) at the station 
neglecting the axisymmetric contributions, and com 


pute 


Aa [2 | 2 tk PY 92 
p Li M2 dim. \*' a ed 


) 


(3) Then carry out the integrals as indicated in (21 
with the known function 7(v) and ¢(v) dr dx i ra 
is taken as a streamline. Then v(x) is obtained, from 
which we obtain the corrected value of B(x). A new 
average value of 8 may be obtained by replacing 
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in Eq. (23) with the corrected value of 3(4 
nd the process repeated, if circumstances so require. 
Sample calculations show that this is quite unnecessary. 

For arbitrary body shapes, the integrals in Eq. (21 
must be evaluated numerically. The following form of 


Eq. (21) is probably more convenient 


‘ go) + (1/2B)Je — 
2)I2 4 1 B) + xo] In [r(x)/r(x 


21a) 


vhere 


F dx : dx 
] and J y 
ve 743 2 ra 


ie purely the geometrical properties of the body shape 
that determine the axisymmetric corrections for flow 
properties on the body surface according to the present 
ipproximate scheme. 

For a cylindrical or conical portion of the body, rela 


tion (20) will be adopted, giving explicitly 


= }4k,/[(2n/6r) — 1]; X 
, Ry a vo) /Br]t) 24) 


where subscript 0 indicates the station where the conical 
or cvlinder portion begins. The value of 1 is deter 


mined by 
Vo — tk, [(2n Br) — 1] (25) 


with the value of y.. taken from the appropriate down- 
stream asymptotic value of v for an infinitely extended 
conical or cylinder body. 

As a final remark, the above specific selections of 
semiempirical rules are made primarily for the purpose 
of simplicity and convenience, not with the purpose of 
giving results in best agreement with the exact solu- 
tions available. The accuracy of the results obtained 
with the above scheme, the effects of the uncertainties, 
and the arbitrariness in the selection of the mean 


values are reported in the next section. 


53) ACCURACY OF THE INTEGRAL METHOD AND OTHER 
PRACTICAL CONSIDERATIONS 


We shall consider the complete characteristic calcu- 
lation of the inviscid flow of an axisymmetric body, in- 
cluding the effects of a curved shock wave at the nose, of 
the wave interactions at the shock, and of the entropy 
gradient across the streamlines as the ‘‘exact inviscid 
solution The deviations of the pressure ec: efficients, 

p,) p,, as evaluated from different approximate 
methods, from the values given by the “exact inviscid 
solution’’ will be considered as the errors of the ap 
proximate methods. The body shapes selected for de- 
tailed investigation are ogive cylinders. As previously 
noted, the results in the cylindrical portion will serve 
is a very critical test of the present approximation be- 
cause we expect analytically that the prsent approxima- 
tion is likely to be the least accurate. The present 
author fails to find any ‘“‘complete characteristic solu- 
tions’ for the flow field with negative flow deflection 
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(toward the axis of symmetry) such as about the down 
stream half of the body of revolution of a circular are 
i.e., biconvex body. 

Charts and tables in reference 6 are extensively used 
both for the conical shock solution and for the isentropic 
planar relations. Linear interpolation of the tabulated 
isentropic relations is found necessary 1n order to main 
tain the slide rule accuracy of three significant figures 
throughout the calculation. Accordingly, two decimals 
in the pressure coefficient calculated from the present 
integral method are considered as correct in the curves 
reported herein. The result of the characteristic caleu 
lation as is available to the present author are in the 
form of curves of similar scale as reproduced in the 
graphs of this report. Accordingly, no attempt has 
been made to improve the accuracy of the present cal 
culation under the integral scheme. 

In Fig. 3 the results of a sample calculation are sum 


marized. The body under consideration is an ogive 


cylinder of slenderness ratio / d 
diameter of the cylinder and / the axial length of the 
ogive nose. The flight Mach Number is /1/, 3. The 
ordinary hypersonic similarity parameter .J/,/(//d) is 
The ordinate in the curve is the pressure co 


3 where d is the 


unity. 
efficient expressed as the excess static pressure Dp p 
over the undisturbed free-stream static pressure p,. 
The abscissa is the axial distance downstream of the 
nose in per cent of nose length / of the ogive. Thus, 100 
on the abscissa indicates the station where the ogive 
joins with the cylinder with discontinuous radius of 
curvature, but with continuous slope. The angle mark 
ing each set of computed points 1s the angle representing 
the are length along the meridian section from the local 
point to the p unt of 100 per cent nose length where the 
cylinder begins. 

The triangles represent the pressure coeflicients as 
calculated from the tangent cone approximation. The 
curve joining these triangles should pass through the 
1.64, at the nose a 0 and 


point, Cp 
100 seyond this 


the point c, (p—p.)/p Oata 
point, x > 100, the tangent cone approximation, if ex 
trapolated, would give a zero pressure coefficient along 
the surface of the cylinder. Starting from the 14° sta 


) nose, where the 


tion—roughly 25 per cent from the 
tangent cone approximation is taken for the determina 
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tion of the initial values—the pressure coefficients at 
different downstream stations were also determined, 
based on the extreme assumption that the flow field 
downstream of the 14° station is planar, governed by 
the Prandtl-Meyer relation. 
sented by X up to the 0° station (v = 100). This two 
dimensional approximation gives a pressure coefficient 


These points are repre- 


over the entire length of the cylinder which is the same 
as that at 0° station. The dotted curve is taken from 
reference 7, which is the result of a characteristic cal 
culation without rotation term in the characteristic 
relation. The solid curve is taken from reference 5 
which is a result of a characteristic calculation, with 
rotation term in the characteristic relation taken into 
account. Superposed on the graphs are the two c, 
curves as determined from the uncorrected Newtonian 
approximation and from the shock expansion method. 
The Newtonian approximation results in an under- 
estimate of c, near the nose and gives a zero pressure co- 
efficient at the maximum diameter point (x 100) 
that is, a local overestimate of c,. The Newtonian c, 
curve crosses the characteristic c, curves in the up- 
stream region. Clearly, it cannot be expected to serve 
as a good approximation in the flow regions with negative 
deflection angles locally. The shock expansion method 
gives a very good approximation near the nose, slightly 
underestimating the c, consistently. But toward the 
downstream region, say near the maximum diameter 
point (x 100), the c, curve obtained from the shock 
expansion method agrees well with the locally planar 
curve but deviates further away from the ‘‘complete 
characteristic solution.’’ Therefore, we cannot ex- 
pect the shock expansion method to give a better than 
planar approximation for flow fields with negative de- 
flection angle, further downstream. 

The results of the present integral approximation with 
initial station 14° where the tangent cone approxima- 
tion is adopted, and with B evaluated as the arithmetic 
mean according to Eq. (23) arerepresented by ©. The 
cylinder portion x > 100 was calculated with Eq. (24) 
and » = 5/4, selected to have approximately the 
downstream asymptotic behavior of the zero pressure 
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coefficient along the cylinder surface based on Eq. (25 
No curve is drawn to joint these points because such g 
curve would be too close to and hardly distinguishab} 
from the solid curve representing the results of th 
“complete characteristic calculation.” 

It might seem that the good agreement could be dy 
to the overestimate of the initial pressure coefficient aj 
the initial station under the tangent cone approxima. 
tion. If the correct pressure coefficient from the ¢ harac 
teristic calculation were used as the initial value, the 
present approximate results would be displaced to lower 
values of c, all the way downstream. This is not so be 
cause, when the initial pressure is taken from th 
characteristic calculation, the value of 1/s, and ther 
fore ko, should also be evaluated from the results of g 
characteristic calculation locally, not from the tangent 
cone approximation. Whether a consistent use of the 
results of the characteristic calculation in evaluating the 
initial values will permit the present integral approxi- 
mation to give as good an agreement as is indicated her 
is open to question and is not to be resolved here. What 
seems to be the important question is the possibility 
that the good agreement might just be an accident of 
the particular choice of the initial station or other uncer 
tain factors in the present approximate scheme. We 
shall, therefore, look into these aspects. 

It is natural to expect that the errors committed by 
these uncertainties will be magnified and become more 
significant in regions further downstream of the + = 
100 point. To avoid the added complication and un- 
certainty of the cylinder portion, the pressure coefficient 
c, over the body of revolution of the circular arc asa 
continuation of the nose ogive is calculated with the 
present method even though there is no complete 
characteristic solution for such a biconvex body avail- 
able for comparison. The results of the present approxi- 
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ulation are plotted in Fig. 4 with a slightly en- 
The result of the present calculation 


curve to flatten out 


mate ¢ al 
larged ordinate. 
shows a strong tendency of the c, 
in the downstream half of the body because of the three 
jimensional effect. The differences among the two- 
dimensional approximation, the tangent cone approxt- 
mation, and the present result become large. As ex- 
nected on physical grounds, the recompression due to 
three-dimensional effects is milder in the present case 
{a continuous are than in the case of a cylinder. 

The effect of different initial stations in starting the 
present approximate solution is shown in Fig. 5 for the 
case of a body of revolution of a circular arc. Three 
sets of calculations were made with initial stations at 
16°, 14 and 12 are different for 


different initial stations, and the values of B are also 


The values of x 


different at a given value of x in the calculations with 


different initial stations. It is clear from Fig. 5 that 


the calculated results based on the present approximate 
method is not sensitive to the choice of the initial sta- 
tion. The divergence of the results with different initial 
stations is somewhat irregular, but the magnitude of the 
divergence lies mostly in the third significant figure of 
the calculation. This is considerably smaller than the 
difference between the present result and the tangent 
cone or the two-dimensional approximation, especially 
in the flow region with negative flow deflection where no 
results from characteristic calculation are available. 
For the ogive nose portion, the variation is within the 
two bounds of characteristic calculations, closer to the 
me including the rotation term. It is, therefore, in 
ferred that the effect of an arbitrary selection of the 
initial station in starting the present approximation 
within a little 


Clearly, one should not and will not select the station 


reasonable range is of consequence. 
too far downstream where the tangent cone approxt- 
mation is excessively off. In fact, if one should select 
the initial station at a 100, the present methods would 
ilso predict zero pressure coefficient all the way on the 
ensuing cylinder body. Even under this extreme case, 
the present method still gives, on the downstream half 
of revolution of a circular arc, a pressure 
“correct values’’ de- 


f the be dy 
coeflicient converging toward the 
termined with the proper selection of the initial station. 

Comparison of the three sets of results seems to in- 
dicate that better agreement with the rotation-charac- 
obtained the station 
The solution with the initial 


teristic solution is with initial 


taken further upstream. 
seems to be the closest to the result of 
characteristic calculation. (The half 
But it should be emphatically brought out that 


the difference is so small that such an inference cannot be 


Station of 16 
nose angle is 


18.9 


justified with the slide rule accuracy of the present cal- 
culation. Difficulties will certainly be encountered if 
one should try to start directly from the nose where r 

0. It seems that the upstream limit of selecting 
the initial station depends very much on the numerical 
the pur- 


poses, there appears to be no difficulty in selecting a 


accuracy of calculation. For engineering 


reasonable initial position. 
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Let us now investigate the effect of the different ways 
of determining B. The previous calculations for C 
over the body of revolution of the circular 
peated with all three initial stations when the mean 8 


arc are re 


is taken as the geometric mean of 8; and Be gim (A 


p [3 ; 


instead of the arithmetic mean given previously as Eq. 
(23). The difference of c, at a given station calculated 
with the two different ways of taking 8 is completely 
overshadowed by the numerical accuracy of a 12-in, 


slide rule for practically all the points calculated. 
There are several instances where deviations may ap 
pear in the third significant figure. The arithmetic 


mean of 8° defined as 


has also been tried with insignificant effect. The step 
by step integration with the immediately preceding 
point as the initial station has also been carried out. 
The values of c, so determined are still insignificant in 


the graph presented here. The difference seems to be 


bigger than that between the different averaging 
schemes and tend toward the characteristic solution 
without rotation terms. It 1s difficult to resolve 


whether this is due to the accumulation of inaccuracies 
of the calculation at successive points, or due to the 
self-compensation of different errors when average B 
is taken. The important thing, for the time being, 1s 
that any of these schemes will give substantially the 
Since the arithmetic mean can be handled 
may 


same result. 
more easily, the scheme as indicated by Eq. (25) 
be adopted in general. 

Let us consider the arbitrariness of the selection of 
for the calculation on the cylinder portion. The physi 
—Br(Oln«/Olnr)4 
Calcu- 


cal meaning of 7 is indicated by 1 
along the positive family characteristic line. 
lated results on the cylinder portion of the ogive cylinder 
shown in Fig. 3 are exemplified in Fig. 6 for three dif- 
ferent values of m = 0, 5/4, and 7 4, indicated by lines 
joining points marked DO, The 
case of nm = 0 represents a flow over the cylinder with 
the flow deflection taking a constant value over the 


>, and A, respectively. 














192 FOURNAL OF THE AERONAT 
CONE - CYLINDER 
8 - 
—e Me=3 — 
a 6} 
Aw 5 E ———— it 5S 44 
a | j 
.ct | 
3 | | - 
; SI 
S ot 9 
q | 
It 
} 
-J 
-4 — + 
oO 100 00 300 


Percent nose length 


FG. ( 


entire flow field downstream of the Mach cone at the 


LOO. 


initial condition at x 


shoulder x 
100 or the boundary condition 
on the cylindrical surface. Therefore, the calculated 
result over the cylinder with 7 0 cannot be expected 
as a satisfactory approximation in the overall picture. 
Some finite value of 7 must be selected. In the sample 
calculations for the ogive cylinder combination, as 
shown in Fig. 6, both the results with 7 5/4 and n 

7/4 give a satisfactory approximation to the character- 
istic solution with rotation terms up to one nose length 
downstream. The calculation was made with Eq. (24) 
and the asymptotic values of c, in the far downstream 
region are also shown in Fig. 6. The asymptotic value 
of c, with wm = 5/4is slightly positive and that with x 

7/4 is slightly negative, both fairly close to zero as they 
should be in the very far downstream region. The 
curve of c, as calculated with 7 7/4 as a function of 
x is nearly parallel to that of the characteristic calcula- 
t seems 


tion. Thus, for this particular example, 7 = 7 


to give a better approximation than 7 0/4, even 
though x 5/4 seems to give a smaller absolute devia- 
tion from the characteristic solution in the region with 
This observation certainly cannot 
be extrapolated to other body shapes. What 
portant is that either of the results with 7 5 


x < 200 per cent. 
is im- 
$ and 
t, as suggested by the rough approximation, rela- 
This 


aspect cannot be confidently extrapolated to cylinder 


n 7 
tion (25), may be adopted without serious error. 


bodies whose nose shapes are quite different from an 


ogive. For illustrative purposes, the cone-cylinder 


combination has been investigated. The results are pre- 
sented in Fig. 7. Within one nose length of the cylinder 
body, the calculated results using 7 QO and x 5/4 
may all be considered as fair approximations. The re- 
sult with 7 0 seems to give the least absolute devia- 
tion. It is, therefore, difficult even to make an intelli- 
gent guess as to what value of 7 should be taken in order 
to get reasonable approximations over the cylinders 
with arbitrary nose for common engineering problems. 
If one is willing to accept the cone-cylinder as an ex- 
static 


treme configuration, and tolerate an error in 





This situation contradicts either the 


T 
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pressure by 10 per cent, it really does not matter whay 
values of 7 one takes, especially if guided by a rule sy h 
as Eq. (25). 
As an additional example, th Me Ogir linder 
combination flying at .1/ 6 has also bee lyzed 


with the results given in Fig. 8S. The scale for the pres. 


sure coellicient c, in this plot is much coarse 


in Fig. 3, for example While graphically, the agree. 
ment between the present approximate solut 

the characteristic solution with rotation terms in the 
characteristic relations appear to be equally good, as js 
shown in Fig. 3 for .\/ 3, the deviation of the actu 
magnitude is considerably larger in the higher Mach 
Number case. There are several other aspects that are 
worth mentioning concerning this specific example 
First, the rotation term in the characteristic relatioy 
representing the effect of the entropy gradient across 


streamlines appears to be very important in the charac- 
teristic calculation. The error in neglecting this term 
than that of any of the 


he effect of 


is considerably bigger approxi- 


mate methods. This important part of t 
rotationality in determining the surface pressure is 
a general fashion in the present ap- 
Second, the 


accounted for in 
proximate method by the parameter «x. 
two-dimensional approximation becomes very good at 
Mach Number .J/, = 6 (Rk 
at M, 3 (k Lb). 
observations made previously in this report [see Eq 
(19) J. | 


2) but not nearly so good 


This is in agreement with the 


The foregoing numerical evidence of the accuracy of 
the present approximate method when applied to an 
ogive body is quite encouraging, especially in view of the 
insensitivity of the results of the present approximation 
toward uncertainties implicit in the method. Further 
numerical evidence as to the accuracy of the present 
method compared with the complete characteristic 
calculations when applied to other pointed bodies and 
for results in regions of negative flow deflection is ex- 
tremely desirable, if the applicability of the present ap- 
proximate scheme is to be carefully explored. The 
agreement as shown in Fig. 3 may not be obtained for 
other pointed bodies. Ordinary engineering practice 
however, will tolerate errors considerably larger than 
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those indicated in Fig. 5. Hence, some extrapolation in 


ive may be permitted. The advantage of the 
sent approximate method lies primarily in its sim- 
it) compared with the characteristic method, 
ie fact that the surface pressure can be deter 
ithout carrying along the flow field in its vicinity 
r even the pressure at intermediate points. Most im 
rtant of all, the previous investigation justifies the 
troduction of an overall parameter «x {1 2/s)le 
to represent the integrated effect of the upstream ge- 
metry in evaluating the axisymmetric corrections. 


[his makes it possible to study the downstream 


nortion of a flow field without referring to the details 
of the upstream configuration. Thus, the major object 


‘the present investigation 1s achieved.* 


CONCLUSIONS 


[he approximate method developed analytically and 
investigated numerically in this report is shown to be 
very satisfactory for pointed ogive bodies and may be 
ypplied with reasonable assurance to other curved, 
pointed nose convex bodies with axisymmetric inviscid 
flow fields, in flow regions with streamlines either diverg 
ing from or converging toward the axis of symmetry. 
The result of the approximate method is insensitive to 
the arbitrary factors implicit in the method. 

For an axisymmetric flow field over prolonged conical 
or cylindrical portions, an alternate formulation 1s pro- 
vided, which may be considered as a fair engineering 
ipproximation. 

Calculations indicate that a locally planar approxi- 
mation is good at relatively large values of the hypersonic 
similarity parameter, say J/g > 2, but that the axisym- 
metric correction is quite significant at J/g <1. Anin- 
The author was informed of reference 8 after the present 

nuscript had been completed. In this reference a second 
rder correction to the shock expansion method is introduced 


original shock expansion method, it will depend on the 


= 


details of the upstream geometry 


FRIC INVISC 


ving the present scheme to pointed bodies similar 
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ference from the analytic formulation of the present 


method is that the axisymmetric correction is likely to 
be more significant at the lower values of the hypersonic 


similarity parameters. 

The integrated effect ot the detailed configur 1uion up 
stream of a given station on the inviscid flow field in 
the downstream region may be represented by an overall 
parameter k 2 ¢. With the introduction of 


this parameter it 1s possible to study an isolated portion 
of an axisymmetric flow field in the downstream region 
without referring to the details of the upstream cor 


figurations 
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Incompressible Two-Dimensional Stagnation- 
Point Flow of an Electrically Conducting 
Viscous Fluid in the Presence of a Magnetic 


Field 


JOSEPH L. NEURINGER* ann WILLIAM McILROY** 


Republic Aviation Corporation 


SUMMARY 


When a blunt body moves through the atmosphere at high 


hypersonic speeds, there exists a 
shock and the nose of the body. It 
possible that, by impressing a magnetic field out of the nose of the 


cap of ionized gas between the 


detached was considered 


body and into the conducting gas, some beneficial effects could 


be derived from this region. In order to demonstrate this fact 
theoretically, consideration was given to the well-known problem 
of the two-dimensional stagnation-point flow on an infinite flat 
plate, modified in two ways: 1) the fluid was assumed incom 
impressed on the fluid was an 


From the 


pressible but conducting, and 
external parallel magnetic field normal to the plate 

resulting analysis of the hydromagnetic interaction, it 
that a skin 


with reasonable values of the applied magnetic intensity 


ippears 


considerable reduction in friction may be possible 


and 
electrical conductivity of the fluid 
The reduction in shear stress at the wall is found to be a func 


tion of a nondimensional parameter /?(0 defined in the text 


SYMBOLS 


Rational m. k. s. system is used throughout 
E* = electric intensity 
Pf = magnetic intensity 
B* magnetic induction 
F = current density 
0 fluid conductivity 
6 charge density 
7 = magnetic permeability 
€ dielectric constant 
q° = fluid velocity 
Pp" fluid pressure 
p fluid density 
be coefficient of viscosity 
i = coefficient of kinematic viscosity 
C* 9 coordinates 
b,b = parameters (dimensions sec 
1/pu,o 
n/t 
9 velocity components in x* and y* directions, 
respectively 
| Pigras & § magnetic intensity components in x* and y* 


directions, respectively 
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i(y), hy), ¢ ; 

Pi(y), Pay) dimensionless functions of the dimension- 
less y coordinate (see Analysis 

.* gradient 

i unit vectors 

The starred quantities are dimensional, the unstarr: 


dimensional 


INTRODUCTION 


| igen RE-ENTRY MISSILES, especially those in 
the intercontinental class, re-enter the earth's 
atmosphere at such tremendous speeds that the gas 
behind the shock may become highly dissociated and 
ionized. Jonization renders the gas electrically con- 
ducting and hence capable of being affected by electro- 
magnetic fields. The name given to the study of the 
interaction of an electrically conducting moving fluid 
in the presence of magnetic fields is magnetohydrody- 
namics. In its simplest terms, the interaction may be 
described as follows. The motion of the conducting 
fluid in the magnetic field induces electric currents in 
the fluid, thereby modifying the field; at the same time 
the flow in the magnetic field produces mechanical 
forces which in turn modify the motion. 

Up until the present, magnetohydrodynamics was 
virtually the sole concern of the geophysicist and _ the 
astrophysicist, since the fluid in the earth’s interior, and 
the gases in and around celestial bodies and in inter- 
stellar space, are already highly ionized and capable 
of interacting with the magnetic fields of these bodies. 
However, now that the behind the 
conducting, the subject becomes of interest to the missile 


gases shock are 
aerodynamicist. 

One of the chief problems associated with the re- 
entry of an intercontinental ballistic missile is that ol 
convective heating. However, since the gas near the 
stagnation region is conducting, the question may be 
raised as to whether the flow in this region can be suf- 
ficiently influenced magnetohydrodynamically so as to 
effectively reduce the heat transfer. In this report, 
we will not be concerned with the more complex prob- 
lem of heat transfer, but rather with skin friction. I, 
as is usually the case, reduced skin friction implies 
reduced heat transfer, then a marked reduction in skin 
friction hydromagnetically should also have the same 
favorable effect on heat transfer. 
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APPLIED MAGNETIC INTENSITY H, 


Stagnation flow of an ionized gas on an infinite plate in 


presence of a magnetic field 


In order to demonstrate the interaction, we choose to 
consider what perhaps is the simplest of the stagnation 
flow problems, the well-known two-dimensional stag- 
nation-point flow on an infinite plate, modified in two 
ways: (1) the fluid is assumed incompressible but 
conducting, and (2) impressed on the fluid is an exter- 
nal parallel magnetic field normal to the plate, Fig. 1. 
We wish to determine the flow in the vicinity of the 
stagnation point with particular emphasis on skin fric- 
tion. 


THE MAGNETOHYDRODYNAMIC EQUATIONS 


(a) The Electromagnetic Equations 


Assuming uw, and ¢ constant in the conducting fluid, 


Maxwell's equations in the rational m. k. s. system be- 


come 
v*xs" — (OB*/ Ot) (1) 

Vv" 2" 6/€ 2) 

V*-3" 0 (3) 

V*s J* + «€(0E*/Ot) (4) 

where B* u,H*. The most general expression for 


J* in a homogeneous isotropic medium is 


J* = cE* + o(q* x B*) + dq’ 5) 


where the terms on the right represent, respectively, 
the conduction current, the induction current, and the 
convection current. The last term of the right-hand 
side of Eq. (4) represents the displacement current. 
It can be shown,' by dimensional and order of magni- 
tude arguments, that the displacement current and 
convection currents are negligible compared with the 
conduction and induction currents, and they are there- 
fore neglected in the analysis. 

Taking the curl of Eqs. (5), using Eqs. (1) and (4), 


and remembering that B* = yu,H*, Eq. (5) becomes 
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— ou,(0H* Ot) + ou.V*x (q* x H*) 
(6) 


y*<i(v" <2") 


Expanding the vectors V* x (V* x A*) and V* x 
(q* x H*) and remembering that V*-q* = 0 (incompres- 


sibility) and V*-H* 0 (Eq. 3), Eq. (6) becomes, 
in steady state, 
(q*-V*)H* — H*.V*)q* nV **H* 7) 


where the Laplacian operator acts on each component 
of H* and n 1 ope. 

It is seen that, as far as the electromagnetic equations 
are concerned, the interaction with the moving con- 
ducting fluid is represented by the current induction 
term in Eq. (5) or the two terms on the left-hand side 


of Eq. (7). 


(b) Flow Equations 


The equations of continuity and momentum for an 
incompressible magnetohydrodynamically interacting 


fluid in steady motion are 
V*-g* = 0 8) 


pe)V *P* 1 
My po) (J* x H* + vV *-q* 


(g*-V*)q* ~(] 


va) 
where the interaction term u,(J* x H*), represents the 
ponderomotive force per unit volume exerted by the 
magnetic field on the fluid. 

Substituting for J* in Eq. (9a) from Eq. (4), we ob- 
tain 


po)V*P* 
M, po) (V* x H*) x H* t vv *-q" (Y) 


ANALYSIS 


(a) Stagnation Flow Equations 


Considering now the two-dimensional stagnation flow 
problem, introduce the dimensionless variables 


\ v*(b vy), y = y*(bv)'?, @ = g*, (bv) 
P*/(pobv), H (H* /H,*) (u,F1o*?/ pobv) 
where //,* represents the magnitude of the applied field, 
and } a parameter, of dimensions sec.~!, whose signifi- 
cance will be discussed later. 


In component form, Eqs. (9) and (7) become 


U(oU ox) + V(Ol' dy) 
—(OP dx) — H,{(Ol7, Ox) — (OH,/dy)] 4 
[(O°L" Ox?) 0°U’ oy 10) 
U(Ol ox) + V(Ol’ oy) 
—(OP/dv) + H,[(Ol1, Ox) — (OH,/ Oy 
[(O? 1" Ox?) + (0° V/Oy’) 11) 
U(OH, Ox) + V(0H, oy) — 
HOU Ox) —H,(0U/dy) = 
y[(O°H, Ox?) + (02H, /dy' 12) 
U(0H,/Ox) + V(OH, dy) — ,(01 ox) - 
IT,(OV/dOy) = y[(O°7T, Ox?) + (0°H,/Oy")] (13) 
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where 
q Ut Vj, A A,i+ H,j, 4 n/v 


Continuity and Eq. (3) yield the auxiliary equations 


(OU /ox) + (OV/oy) = 0; 
(O/T,/Ox) + (OH,/Oy) = 0 (14) 


Since, (1) the « components of both qg and H are 
antisymmetric with respect to v, (2) the vy components 
are symmetric, and (3) both g and H are divergenceless, 
we seek a solution of the system of Eqs. (10) to (14) 


of the form 
l u'(y)x, Tl h’(v)a 
|” —u(y), fH 
P Py — Pily) 


—h(y) (15) 
P(y)x? 


where the prime denotes differentiation with respect to 
y, and u(y), h(y), Pi(y), Po(y) are functions of y to be 
determined. /») represents the dimensionless stagna 
tion pressure. The odd powers of x are missing from 
the P expression because of symmetry with respect to x. 

Substituting Eqs. (15) into the system (10) to (15) 
and equating coefficients of like powers of ., we obtain 


from Eq. (10), 
u’? — un’! 2P.— hh" + u'" (16) 


from Eq. (11), 


uu’ Pr =a" (17) 
Ps h'h’’ (18) 

from Eq. (12), 
hu’’ — uh" yh’’’ (19) 

from Eq. (13), 
uh’ — uh —yh"’ (20) 


Integrating Eq. (18) we obtain 
P» (h’?/2) +C (21 
where C is the integration constant to be evaluated. 


(b) Boundary Conditions 
The velocity conditions at the plate require 
u(O) u’(O) 0 


The total magnetic field at the wall is the vector sum 
of the applied and induced fields, the latter resulting 
from the induced current distribution in the fluid. This 
distribution depends on the solution and is therefore 
unknown. 
effect of the interaction on the flow field, the induced 


However, for the purpose of calculating the 


field at the wall is neglected in comparison with the 
applied field. The magnetic field components at the 


wall thus become 
h(O) — (u,FIy*?/ pybv)' h’(0) 0 


h being obtained from the transformation to dimension- 


less variables. 


As in the field free case, the conditions for P; and | 
at the wall are determined by integrating the invise; 
momentum equation along the wall streamline 
comparing the result obtained with the expression for 
P assumed in Eq. (15). 

Returning to Eq. (9), neglecting the viscous tern 
and expanding (q*-V*)q*, we obtain 


v*[p* 4 


(1/2) poq” 


polg* x (V*xq’" u(V* x H*) x H’ 


Integrating along the wall streamline from 
to x” v*, and remembering that the velocity con 
ponents are zero at the stagnation point and the vertical 
component of g* is zero at the wall, the above expres 


sion becomes 


P*(x*, 0) Po*(0, 0) — po/2)U*?(x*, 0) 4 


0 } \q* x Vv" 2@")|-2 dx" 

ML, } iV *xh*) 2" |-tdx* (22 
The first integral on the right-hand side is identical 
zero. Expanding the last integral, we obtain 


P(e. O P,*(0, 0) — (g9/ 2) **{x",.'0) + 


fe } H,*[(0H,*/Oy*) — (OH,*/dx*) |dx* (23 


We assume the velocity components of the inviscid 
magnetohydrodynamic flow are given by L* by’ 
|"* = —by* since these forms are the simplest repre 
sentation satisfying both continuity and symmetry 
It is seen that the parameter, ), introduced in the di- 
mensionalization, has the same significance as in the 
field free case.*. In dimensionless form, these become 


U=x,V = -y. 


becomes 


In nondimensional form Eq. (23 


P(x, 0) P,(0, 0) — (1/2)x? }4F 


( IT, {(O71,/Ov) — (O/T, Ox) lda 24 
0 
\gain, the components of H are of the form assumed 


in Eq. (15)—4.e., 


IT h'(v)x, HH —h(y) 
Hence, Eq. (24) reduces to 
Pix. 6) Po(O, 0) — (1/2) h(O) h’’(O 
J 
or 
P(x, 0) P.(0, 0) (1/2)f1 h(O)h’’(0 


Now, in the inviscid flow, Eq. (20) is still valid, and 


since O and x’ 1 at the wall, yh’’(0) h(0 
The above expression for the pressure finally becomes 
P(x, 0) P\(O, O) — (1/2); 1 + [h2(0) /y]ixv 
When this is compared with the expression assumed for 


P in Eq. (15), the wall conditions 
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INCOMPRESSIBLE TWO 


TABLE | 





io> =| x 108 5x10° los 

nor | v0) | 7% | uO) [7% jut |%7, [uO |%7% 
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Fic. 2. Shear stress ratio versus parameter /(Q0) fer various 
values of + 


P,(0) 0, P.(0) (1 2)}1 + [h2(0) y]j 


are obtained. 

For the boundary condition at large distances from 
the wall, we require that the horizontal velocity com- 
ponent approach asymptotically the assumed inviscid 
form—i. e., U x. This requires wu’ l as y ap- 
proaches infinity. 

Returning to Eq. (21), and using the condition that 


h’ Q at the wall, the constant C (1,2))1 + 


Substituting for P, in Eq. (16), we obtain 
we’ + un!’ — uw’? +41 + [h2(0)/y]f + 
h’? — hh’’ = 0 (25) 


/ | 


and Eq. (20), yh’’ + uh’ — u’h = 0 (20) 


It is seen that the derivative equation of (20) yields 


Eq. (19); hence, a solution involving Eqs. (25) and 


20) automatically satisfies Eq. (19). Once uw and h 
are obtained from Eqs. (25) and (20), ?: and ?: can be 


obtained from Eqs. (17) and (18), respectively. 


DIMENSIONAL STAGNATION 
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(c) Recapitulation 
In summary, we need to solve the following system 


of equations: 


ul + UU — Hu + ,1 + [a-(0O)/¥Ii 
/ h/ 0) 25) 
yh"’ + uh uh 0 20 
sub }c ct to the conditions: 
At y (0), 
i u’ h’ 0 Hl utd* PyDv 
At y oD, u’ | 


It is noted that when the applied / 0, the equa 
tions uncouple and Eq. (25) reduces identically to the 


corresponding field free equation. 


SOLUTION AND DISCUSSION 
The pressure expressions at the wall with and without 
a field are, respectively, 
P,*(0, O) — (1° 2)pob?) 1 + [h2(0) y]hx 
P,*(0, 0) — (1 


pP*(s* 0) 
P*(s*. @) 


2) pudo*x" 


The assertion is now made that the nature of the 
hydromagnetic force (regarded as a strong magnetic 
“viscosity” tending to destroy the motion across the 
lines of force*) and the character of the interaction are 
such as to transform, by the application of the mag 
netic field, an initial infinite two-dimensional inviscid 
field free flow whose velocity components are charac 
terized by the parameter, ), to a flow whose velocity 
components are characterized by the parameter, ), the 
pressure distribution at the wall remaining essentially 
unaltered. Hence, comparing the above expressions, 
we have 


h by) 1 + h?(0) 7 ; ()) 
Now the shear stress, 7, at the wall is given by 


r = w(Ol* dy*),~o 


u(OU*/Oy),-9 (Ov Ov*) pyv'/*b3/*4""(O)x" 


Comparing wall shears, with and without the field, 
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Fic. 4 


yields 
(7/70) (20'’(O) /tt0’’(O) |(6/ be)?!” 

(u’’(O) wo’’(O) |f.1 + [h2(0), y]f -*/4 (27) 
where uy’’(O) 1.2326 


The solution of the system (25) and (20) was obtained 
on an analog computer and the results for «’’(O) for 
various values of 4(0) and + are shown in Table 1. 
Curves for the shear stress ratio, according to Eq. 
(27), are shown in Fig. 2. If the shear stress ratio is 
plotted as a function of the parameter, /°(0)/y, then 
the curves of Fig. 2 transform into the single universal 
curve represented in Fig. 3. The parameter )°(0)/y, 
in terms of the field free flow parameter, do, and dimen- 
sional quantities, is given by 

h?(O) 2 
¥ (1 + (4607 07/ Bo**a7) |'/? — 1 


kin” 


(2S) 


where By" 


The relation between /°?(0) y and the more fundamen- 
tal parameter, oB.*°, bop), is shown in Fig. 4. 

It is seen from Fig. 3 and Eq. (28) that the reduction 
in shear stress, for a given field free flow, depends on 
both the applied field and the fluid conductivity. This 
reduction depends much more strongly on the im- 
pressed field than on the fluid conductivity. 
magnetohydrodynamic 
At the wall, where the 
As the 
horizontal velocity component increases further away 


Physically the interaction 
may be described as follows. 


fluid velocity is zero, the induction is zero. 


from the wall, the induced currents increase and con- 
sequently the ponderomotive force, u,(J* x H*) in- 


creases. This force, whose direction is opposite to the 
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direction of the horizontal flow and therefore oppose; 
the motion, conspires in the steady state to give 

horizontal velocity component smaller in magnitug 
than would exist if the impressed field were absent 
The net effect is a reduction in the velocity gradient 
at the wall. 

To illustrate the orders of magnitude involved, con 
sider the case h?(0)/y = 1.6. This, according t 
Fig. 3, corresponds to a 30 per cent reduction in wall 
shear stress. Using Eq. (28), and the following values 
for py, d), and o which, for the re-entry problem, appear 
p 16 X 10-‘ke./m:.*, 65 = 5 X 19 
By* becomes & 0.18 webs 


reasonable 


sec.—!', o & 250 mhos/m.; 


m.° or 1,800 gauss. 


CONCLUSION 


By considering the two-dimensional stagnation-point 
flow of an infinite plate in a conducting fluid upon which 
is impressed an external parallel magnetic field normal 
to the plate, it has been demonstrated theoretically 
that a considerable reduction of wall shear stress can 
be effected. 
edly on the strength of the applied field and the con- 


The percentage reduction depends mark- 
ductivity of the fluid. It appears that it may be pos- 
sible to utilize the natural ionization existing in the 
shock layer of a high speed re-entry missile to reduce 
the skin friction magnetically. 

It is interesting to speculate on the effect of the mag- 
netohydrodynamic interaction on heat transfer. As 
was pointed out in the Introduction, if reduced skin 
friction implies reduced heat transfer, (Reynolds 
analogy), then a reduction in skin friction hydromag- 
netically should also have the same favorable effect on 
heat transfer. The situation here, however, is not as 
obvious as it seems, for it must be remarked that he- 
The Joulian 


heating, which increases with increasing distance from 


sides viscous there is Joulian dissipation. 


the wall, may more than offset some favorable thermo- 
mechanical analogy which might be in operation. 

Further, it is interesting to conjecture on the effect 
of the interaction on transition to turbulence. It 1 
well known* that the magnetic lines of force pervading 
the fluid give it a certain degree of rigidity so that the 
onset or growth of any random motions would be 11- 
hibited; also, the fact that the velocity gradient 
reduced may perhaps have some favorable effect 0 
transition. 

Finally, it must be emphasized that the potentialities 
of this interesting new field of magnetohydrodynamis 
and its practical application to hypersonies will require 
much careful experimentation in addition to theoretical 
studies similar to that of this paper. 
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Comments on ‘‘The Rotational Field Behind a 
Curved Shock Wave Calculated by the Method 
of Flux Analysis”’ 


E S L ve 
ipervisory Aeronautical Research Cngineer 
Langley Aeronautical Laboratory, NACA, Langle 7, V 
September 30, 1957 
CHIDA AND YASUHARA! have presented a means for calcu 
lating the rotational flow field behind a curved shock In 
their evaluation of the effect of rotation and the validity of the 
characteristic relation employed, application was first made to the 
flow about a semiconvex airfoil. Calculations were made of the 
shock shape and of the pressure distribution over the airfoil 
Figs. 10 and 12') and were compared with expcrimental results 
designated as being those of Ferri.? The almost exact agreement 


between the calculated and the experimei.tal results led to the 


conclusion! that the relation as developed ‘‘will be valid as a 
formula for giving corrections of effects caused by rotation as 
shown in the example of semiconvex airfoils.” Examination of 
this excellent agreement between the theory and experiment has, 


however, revealed some disturbing features which, as they stand, 
do not lend support to the above conclusion 
onvex, S&S 


Che airfoil considered two-dimensional semi 


thick 


was 


per cent with leading-edge and trailing-edge angles of 


20°, and was at zero angle of attack with .1/ 2.13. For these 
conditions, the pressure coc ficient C, at the I iding edge must be 
How 


leading 


about 0.6 and the initial shock angle must be about 49.6 


ever, reference 1 shows a theoretical value of C, at the 


edge of about 0.85 and an initial shock angle of about 58°; further, 


these theoretical values are shown to be in agreement wit! 


Ferri’s experimental results, thus implying that the experimental 





results contain errors that are identical with the deviations of the 
theoretical results from the correct values. Since this repre 
sents an unusual occurrence, a check was made of Ferri’s original 
Contributors Please Note 
Readers’ Forum items must be confined to the 
equivalent of one page in the Journal. Contribu- 


tions that exceed this limit will be returned to the 
authors for condensation and rewrite. To avoid incon- 
venience and delay, please adhere to the following 
specifications: 

(1) Five, double-spaced, typewritten, manu- 

script pages (8'/, by 11 in.), including 
wide margins, formulas, and headings 
equal one printed page. 
For every illustration, deduct at least 
one-half or as much as one manuscript 
page, depending on the size of the 
illustration. 


(2) 
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work where the pressures are koth tabulated? and plotted.* 3 
This check revealed that the experimental points contained in 
are incorrect. Ferri has confirmed this in a private 
the authors (of 
All the 


(experimental) values must be divided by 1.5 since the scale of 


Fig. 12 
communication in which he suspects that “ 
reference 1) misprinted the scale of the C,, coefficient 
qg in my original paper corresponds to 15 divisions in place of 10; 
this is probably the source of the error.”’” Thus, the theoretical 
results of reference 1 and the correct experimental results? do not 
agree. Moreover, division by 1.5 would still leave discrepancies 
between the experimental values in Fig. 12! and those obtained 
from the tabulations in reference 2 that are too large to be at 
tributed to plotting accuracy For example, the first experi 
mental point in Fig. 12! would be about 0.41 instead of the correct 
value of about 0.46. Further, there remains unexplained the 
almost 10° error in shock angle at the leading edge and its sup 
posed agreement with Ferri’s results, which in fact does not occur 

In a search for an explanation of these discrepancies, the writer 
noted that the values of pressure coefficient and the shock angle 
at the leading edge obtained in reference 1 would be correct if the 
airfoil were at a = H However, this explanation is unsatis 
factory since C,, at the 50 per cent chord would not be essentially 
zero at a = —4° 
for this airfoilat a = 4° and M = 2.13. 

Some clarification of the errors in the theoretical pressures and 


Further, reference 2 presents no pressures 


shock shape and of the errors introduced in the use of Ferri’s 
experimental results is needed, since these errors tend to cast 
doubt upon parts of the paper! which may otherwise be correct. 
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Behind 
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Authors’ Reply 


Shigeo Uchida 
Professor of Fluid Dynamics, Nagoya University, Nagoya, Japan 
November 7, 1957 


HE AUTHORS APPRECIATE Love’s comments on parts of refer- 
ence | because slight but clear mistakes were found in Fig. 
12? 

The application of the Flux-Analysis Method to the 8.8 per cent 
semiconvex airfoil was for the Mach Number J/; = 1.85, not 
VV, = 2.18 as given in reference 1. The pressure coefficient 
C,, at the leading edge is, therefore, about 0.815 and the initial 
shock angle is about 62.0°. These are almost satisfied by the 
original figure in reference 1. 

The calculation was compared with Ferri’s experiment given 
by Fig. 13 of reference 2, multiplying original values by 3/4 in 
stead of 2/3. Incorrectness of the description of the experi 
mental data was caused by this. In order to show the accurate 
result, a corrected figure for the original Fig. 12! is presented 
here. The full line is not altered from the original one excepting 
the part near the leading edge which was not exact in the tracer’s 
work. The chain line shows the Oth order approximation cal 
culated by the shock-expansion method with nonreflecting char- 
acteristics. The broken line is an analytical curve given by 
Ferri? which is coincident with the curve of shock expansion in 
the region of 0.56 S x/] S$ 1.0. The reason for the deviation of 
these curves in the foresection is not clear to the present writer. 

Since the effect of rotation diminishes the pressure coefficient 
as shown in the figure, better results of successive approxima- 
tions will be expected if the start of the present calculation is 


based on Ferri’s analysis 
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Rev ised . 


Errata for reference 1: 

(1) Change J, = 2.18 to JJ; = 1.85 on page 838 and in Figs 
10 and 12 

(2) Replace Fig. 12 with the revised Fig. 12 

(3) Change scale y,/d = 0, 2, 4, , 10, 12, 14 to y,/d = 


, 1.0, 1.2, 1.4 in Fig. 17 


Flows in Partly Dissociated Gases; 


Hans-Joachim Metzdorf* 

Institut fUr Theoretische Physik der Freien Universitat Berlin, Ber 
Germany 

October 17, 1957 


Bees SOLUTION of the laminar boundary-layer equations giver 
by L. Moore! for a compressible medium with variable pro} 
erties under the consideration of dissociation requires a corre¢ 
tion, since next to the possible influence of nonbalanced state of 
dissociation the diffusion effects have not been considered. As 
the dissociation energy Ey, in the range of temperature of partia 
dissociation is much greater than the internal and translationa 
energy of the molecules and atoms, one must expect that the dif- | 
fusion introduces a large influence on the heat-flow vector go | 

It is, therefore, suggested here that the following system ol 
equations be considered for treatment of the flows of partly dis 


sociated biatomic gases: 





dn,/dt = —m toe +n — dieo la 

dns/dt = —Noteg t r2 — dogo lb 

p(dv,/dt) = —D.» + (elle) .c - 
p(du/dt) + pio. = BUp,cllpe — Yo.o 


pb = (m + n)kT 


+t The research discussed in this paper was sponsored by The Glenn ! 
Martin Company, Baltimore 3, Md A detailed report was published by 
the company in May, 1956 

* I would like to express thanks to Prof. Dr. Ludwig for his constant en 
couragement and clear discussions during this work Further thanks 
and gratitude are owed to the Bundesversicherungsanstalt, the Firma 
AEG-Turbinenfabrik, and specially to Messrs. Giese, Luhnau, and Dr 


Nasitta. 
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Hoo - r T Vo.p 2/3)t pa 
ré QD 0/Ox,)o Um 
1/dtio@ 0/dt)d + Up 
Gmilar to the stvle of the relativity theory, we are to sum over 
double indices wherever they enter. The symbols used in 
ib equations are 
the three components (vy 1, 2, 3) of the flows 
velocities, 
components of the vector of the heat-flow density, 
components of the diffusion flow of atoms and 
molecules respectively, 
m.? number of atoms and molecules per unit volume, 
respectively, 
density, 
, change of m, and m» per unit time due to dissociation 
and recombination, 
m viscosity, 
\ heat conductivity of gas without considering dissoci 
ation and diffusion, 
hvdrostatic pressure, 
] absolute temperature, 
Boltzmann constant, 
é average internal energy of a molecule, 
E dissociation energy, 
by 1 for wh = vy, otherwise 0 
m,n masses of an atom and a molecule, respectively, 
k coefficient of thermal diffusion, 
kh, k coefficients of reaction rate, which can be calculated 
from the cross section of the molecules for disso 
ciation reaction, 
K temperature-dependent constant of law of mass 


action: ” Kn? (valid for equilibrium 
Weare to substitute 


1/p) ((3/2)kTn, 4+ 


For the heat flow we are to use approximately 


NT g + (5/2)RT (dig + dre) + 
pki 1/m)die — l/n dog +d o En 2 + d. of 


ie diffusion flow possesses, in the same approximation, th 


rm 


die —(mon?/p) Dyoldizg + Rr(log T),6 
dig mn?/p)Dys[dize + Rr(log T),¢ 


, m+n a + [mno(ms — m)/(m + N2)p| (log p).¢ 


If we treat the problem of a stationary laminar boundary layer 
s Moore has done along a tangentially extended infinite flat 
plate, under the simple assumption that the dissociation equilib 
rium appears only momentarily and that we can set kr = 0 
Without introducing large error, then the system (1) to (4) reduces 
to the usual boundary-layer equations for the flat plate. Only 
the coefficient \ (gg = —AT,¢) and the specific heat cp, possess 
igh maxima in the range of temperature of partial dissociation 
If Cp, is the specific heat of the mixture without consideration of 


the dissociation and cpp the part caused by dissociation, we get: 


« = Cou T+ Cyn 
ind N=A+ pD, Cpp 
With 
Cop = (n*/p2)m,[m,/(m, + me)|,r [En + (5/2RT — € 
rhe heat conductivity, viscosity, ete. for pure nitrogen subjected 


toa pressure of p = 3.105 g. cm.~! sec. ~? has been calculated under 
the consideration of experimental data and certain assumptions 


\s a result it came out that up shows a wave nature only by a 
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ary layer depending on 7 L/V/x | 'p v,y* )dy* without and 
e 0 
with consideration of diffusion, respectively The la and 2a 


curves represent the corresponding velocity profiles ¢' cm.se« 


small amount about an average value in a range of temperature 
from 3,500°K. to 7,500°K 


decreasing temperature and reaches, at 


It increases sharply, however, with 


1,Q00°K 


per cent of the average value. It appears that the assumption 


nearly 200 


up = constant, which had simplified the solution of the boundary 
layer equation, is not justified 

\ccording to a procedure similar to Dorodnizyn’s method the 
boundary-layer equations have been reduced to a system of 
ordinary differential equations. The coefficients concerned art 
given in table form 

With the help of an electronic digital computer, the equations 
for the case of an isolated plate have been solved for special 
boundary conditions, once with and then without the consider 
ation of the influence of diffusion on the heat conductivity In 
the second case, the boundary values have been chosen in such a 
way that the temperature of the wall and the nuclear flow take 


the same value as in the first case (see Fig. 1] 
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On Axisymmetric Thermal Stresses in Thin 
Shells of Revolution 


G. D. Galletly 
Mechanical and Electrical Engineering Department 
Shell Development Company, Emeryville, Calif 
October 21, 1957 
I THEIR RECENT BOOK on thermal stresses,!' Melan and Parkus 
give two coupled linear ordinary differential equations for cal 
culating the stresses in a thin shell of revolution due to axisym 
metric temperature distributions. Since these equations may 
be, or may have been, used for calculating thermal stresses in 
shells of varying wall thickness, the purpose of this note is to call 
attention to an apparent error in them. Only for shells of con 
stant wall thickness will their equations yield correct results 
The notation used by Melan and Parkus will also be employed 
herein except that a prime will be used to indicate differentiation 
with respect to @ 
With the radial displacement, w, measured positive outward, 
the angle of rotation of a tangent to a meridian, x, is given by 
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x = l/r v— w’ 1 
In Eq. (1) x is positive when the positive semitangent turns 


inward. The displacement along the meridian, v, is positive when 


directed in the direction of incre ising Dp 
The horizontal displacement, “, (considered positive outward 


is related tovand w by 
f w sind + vecos@ 2 
From Eqs. (1) and (2) we find 
XK l/r, sin ® “ucot gd — u colt @/? vcot d — 7 >) 
From Timoshenko,? we have 
col od r ro €00 ri €ge } 


The horizontal displacement, “4, is also given by 


re S11 @D €99 » 
Now it can be shown that 
, l ro/r,) Cot d 6 
Using Eqs. (5) and (6) it will be found that 
uCcOL @® = nu’ re F COS Meng — ? sin @ (€9¢ 7 


Hence, from Eqs. (3), (4) and (7), we obtain 


XK Ee — €99) COC dB — (Pfo/? €90 Ss 


Eq. (8) is the compatibility equation relating the middle surface 


strains, egg and egs, and the rotation of the meridional tangent, x 





Utilizing Eq. (8), and proceeding in the manner outlined in 


reference 1, vields the following differential equations 


x" 4 y y y y - COE @ 7 36/5 X 
' r/? cot o|? + ur) /? 1 — (36'/5) cot @ } Xx 
— (6"r,2/4A? O + a(l + p)r 36'/6)7r + 7’ q 
or t ry/? ro/? + cot d — (6'/6 - <= 
{ r/o) COL ® = (ors? 1 + (6/6) cot @ ; @) 
1F/5)n[(n/re) x ta T 10 


The left-hand sides of Eqs. (9) and (10) agree with those given 
in Timoshenko* for shells of revolution subjected to axisymmetric 
edge bending loads. The right-hand sides of the equations agree 
in essence with those given by Melan and Parkus; however, the 
above right-hand sides are simpler. The simplification occurs as 
a result of using Eq. (6 

It will be seen that the left-hand side of Eq. (10) differs from 
the corresponding expression given by Melan and Parkus. The 
author believes that they introduce an error here by using the 


following compatibility equation: 
x = €¢o r rijega| COt Gd — l/r ro€90 11 


The author believes Eq. (8) is the correct compatibility equa 
tion and not Eq. (11 

Denoting the value of x given by Eq. (8) by x, then Eq. (11 
gives 


X Xx — (2v/r,) cot? @ 12 


Usually, the tangential displacement, v, is small for shell prob 
lems and the error will probably not be important in many 
cases. However, near the apex of the shel! (nose of a missile 
cot @ becomes very large. The error may be of importance in this 
region; however, as v will be quite small near the apex, it is diffi 
cult to make definite statements even for this case. Resort 
will probably have to be made to numerical comparisons of spe 
cific examples. Whatever the result of this comparison, it is 
nevertheless theoretically desirable to use the correct and some 
what simpler equations. 

It is to be noted that in the special case of constant wall thick- 
ness (which is of importance in practice) the equations given 
herein coincide with those given by Melan and Parkus 
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On the Electric Tank Method for Study of 
Centrifugal Machines 


Syogo Matsunaga 


University of Osaka Prefecture, Sakai, Osaka, Japa 
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A * APPROXIMATE METHOD for study of nonaxisymmetric ff 
4 by the electric tank method is described Wit t 


method, the difference of circulation between blade nd 
radial thrust of runners will be roughly estimated fron 

tion of the streamlines around the blades Phe absolute fl 
runners is decomposed into three potential flows——viz., discharg 


flow caused by the centrifugal force, the displacement flow its¢ 
and circulation around the blade These flows are combi 
the electric tank so as to fulfill the boundary conditiv 





casing wall automatically 
The streamlines of the fluid flow through the blade 


caused by the centrifugal force acting on the fluid rhe vel 

is the velocity caused by the partial discharge quantity 0 
shown in Fig. 1. The total flow quantity 2 Q, distributed alor 
the circle around the eye of the pump gives rise 
sum Z. ( he streamlines caused by centrifugal force acting 
on the fluid must satisfy the boundary condition of the bl 


surface ,soa distribution of doublets on the blade surface must t 
added, 
} 


The velocity V, in Fig. 1 is the velocity caused b 


tributed doublets. The local strength of the doublet distributio 
must be chosen so as to satisfy the condition that the vector dire 


tion of vg be coincident with the tangent line of the blade cury 


a 


where vg is the vector sum of Vs; and 2 Qn. 


When the distribution of Q, 1s given, the distribution of d 


lets on the blade will be known, and then the velocity distributi 


















Distributed 
Sink and 
Source 
on the 
Blade 
Surface 
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raw-Hill p n the pump can be determined. The author considered the 
int P on the extension line of the blade near its edge (Fig. 1 

The velocity V is caused by the centrifugal force acting on the 
,and its component V, is perpendicular to the extension line 

blade (Fig. 3 Streamline through runner 

. he electric tank method for study of the fluid flow caused caused by the centrifugal force 
tudy of the rotating displacement at the point P is shown in Fig. 2, 
ere Itage difference added on the upper and lower surface 
each blade represents a distribution of doublets of approxi 
itely constant strength This approximation will be useful 
for the point near the blade edge. The streamlines caused by 
rotating displacement of the blade can be obtained as the 
a urves perpendicular to the equipotential lines. Using the nor 
Wit : nal velocity component rw-cos 8 normal to the blade surface 
HA > determine the relation between applied voltage and fluid 
— velocities, the velocity Vp at the point P can be determined by 

wn ea the stre unline, and hence the normal component, Von . 

ised The velocity v at the point P caused by the circulation is ob 

flow it oa tained by the method shown in the lower part of Fig. 3 rhe 
velocity is the velocity component of v normal to the extensioa 

mpimner —_ ‘ ‘ 
ine of the blade curve. The velocity v’ at the point P on the ex 
tension line of blade curve of one blade is caused by the circula 

. tion around another blade. v,’ is the component normal to the 
= xtension line of the blade curve 
ac a The fluid flow must go smoothly along the blade curve at the 
is exit of the bladed wheel, and this condition can be written as 
ROG aang follows: the velocity component normal to the extension curve ~— 
rise { f the blade must be zero at the point P. The relations can be 
— ritten at the point P of blade A as follows, Fic. 3 
bs te ValA 7 Vp 4+ (Un)A + (@rp: COS Bp) a T 
a Un’ )Ba—ra + (tn’)c—ra + (Un’) pa 4 = () The above equation for blade A contains the unknown vari 
— : ables I'4, I's, Tc, ete., and similar equations can be obtained for 
ss rhe sign of each term must be determined ” the flow direction, blades B, C, D etc., respectively ia acllcieeint iktamiltlaain 
iia ei 5 ee: pasting va ee oe as rs, Pc, oe., com Se ontained hy the culation of the group of 
aa tii ssid a £: : phen id equations mentioned above It will be noticed that the inter 
or direc- round the blade B on the point P of the blade A The values Maleaist . , : re 3 ; 

i ania fe), FE Nan pe ree eS RET HORM ferences of the circulations around the blade s are not considered 
a ee ae 7 j gh z because of the experimental difficulty of the electric tank method 
ly, Tz, Te, ete., having unknown magnitude This relation 

wa tween (v,)4 and Ty ete. can be obtained by the electric tank 
4 dout | 
ribution — Ms 
Lift and Rolling Moment on a Siender 
Wing-Body Combination With a Deflected 
Tip Control 
John F. Clarke 
English Electric Co. Ltd., Guided Weapons Division, Lutor 
Bedfordshir -, England 
October 24, 1957 

ed 

A 

e SUMMARY 

the Lift and rolling moment acting on a wing-body combination with a single 

ide deflected wing-tip control have been calculated on the basis of the slender 

rfate - body approximation 

7 INTRODUCTION 

"THE MUTUAL INTERFERENCE between the components of a 
wing-body tip-control configuration of small aspect ratio may 
equi-Potential be significant at Mach Numbers up to low supersonic values 
lines Some estimate of this interference can, in this case, be obtained 
Winns from a slender body theory We may take the configuration 

sketched in Fig. 1 as a fundamental element of such a system 
i.e., the outboard section (y < v2) of a flat plate wing on a 


circular cross-section body (radius a) is deflected through an 
angle a. The wing body is at zero angle of attack to a main 
stream of velocity / 

Ward! gives results for the wing body at angle of attack while 
Mirels? has dealt with the case of both wings deflected on a body 
at zero angle of attack. Superposition of these results on those 
to be given below will allow very general solutions to be built 
Fic. 2. up within the framework of the slender body approximation 











204 JOURNAL OF THE AERONAUTICAL SCIENCES MARCH, 1958 
>} 
wo Yy=0; Z=2aq t')dé 
J } 
VELOCITIES >} 
Vv Lia = ag Mé EG 
J} 
Ye where is dynamic pressure, pol*/2, and L(a) is t 
Y, moment due to control deflection a Phe function 
&( £° fa 2\( ¥2? rg 
: ‘ . Eq. (5) can be integrated in terms of simple function 
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Fic. 1. Physical and transformed plane configurations 


RESULTS 


Phe physical cross-flow plane (¢ plane) is transformed into the 


¢ plane by writing ¢ = ¢ + a’ see Fig. 1). The velocity field 


in the ¢ plane can be shown to be given generally by 


iw’ = 
>} 
( l/r ¥ } c?) Wwe’) [V1 § &’)] dé’ 
e } 

1) 
where ¥2 = yp a*y ¢=é in, and W(£) is the downwash 
distribution between ) é ) In the present case 
W(t) = 0; —} £<y;: yy =0 ) 
W(t) = aWU/2)|1 (t/¥ # $a?) 2 

} & Jo: n= a\ 
and }; = » + a*y, 


Ward! has shown that the complex normal force F acting on a 
slender body of any cross section is 


9 


F = 2rpoUa, + po U*(d/dx) [t,(x) S(x 3 


where px is undisturbed stream density, S is body cross-sectional 


area at station x and t, = y, + iz, is the (complex) co-ordinate 
of the body’s center of area at station x. The coefficient a; is 


derived from the expansion, ! 
= a(x) loget + d(x) +4 ay(x)t 


where ¢ is the velocity potential of the fluid motion created by the 


slender body. ao(x) = US'/27 for all body shapes (the dash 
indicating differentiation with respect to x) Using similar 


(momentum) methods we may show the rolling moment L to be 


L=@€8 2r pa Uae p 


U2(d/dx) [yyso5 + I 4) 


where J is the product moment of body cross section taken about 
axes Gy’, Ge’ parallel to Oy, Oz and passing through the center of 
areaG. For the body shape treated here J = 0; alsot, = 0 since 
the body is at zero angle of attack. We may find a, de directly 
from series expansion of Eq. (1) in inverse powers of ¢, using 


the relation ¢ = t + a?t 
that (N B. F = Y + 712), 


1 It is found from Eqs. (8) and (4) 


Z/aq = Y¥2?[(2r/2 sin Vi/3 
} Ja") sin V () } ! 
vo fy a 
y eg” ip 2 Vv} 

Noting that 7 sin~'(2a/y2) = 2 cot '(a/y2), Eq. (7) simy 
fies in the case of the single all-moving wing (}, = 2 1 giv 
a value precisely half that found by Mirels* for both wings d 
flected Presumably the term in parentheses multiplying ¢ 
a/ yz) in reference (2) should be squared. It can then 


that Dugan’s result! is identical with Mirels’ 
Eq 6) is rather more complicated but can be reduced t Lege 


dre elliptic integrals of the first and second kinds plu 


functions. Then 
SLia)/ag = () } 
) ) 7? R< vV } 1a? V/ } 
fa? V2[A(k F(0, k) | + Vox ¥22 — 8a2) | E(k E(6, 8 


where F(6, k) and A(k) are incomplete and complete elliptic i 


tegrals of the first kind, respectively, and similarly E(@, k 


E(k) are second kind integrals. The angle @is given by & sin 9 = 
l ta? Vy 2 and k is the modulus, k? = 1 ta*V.~?. Sor 
values of —L(a)/aqy2* are plotted versus (a in Fig. 2. 7 


rolling moment produced by deflecting an all-moving wing 


—~ ().3 


maximum when (a/ ys) 
The rolling moment L(p) due to rolling at angular velocity 


can be found by similar methods to those used for Lia I 
result is 
2( gp i Pp) = V34 cot ~! (a/v fra’ 
a( }.? Sa?) Y } 
Values of UL(p)/qpy2* are plotted in Fig. 2. Equilibriun 
roll rate is achieved when L(a@) = L(p) and Eqs. (8) and (9 
be used to define a rolling effectiveness parameter 8 = Q 
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Fic. 2. Rolling moment coefficients versus (body radius/sem! 


span) ratio 
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nation of nitrogen dioxide into tetroxide, 





ture is passed through an electric heater and fed into 


ture For smal 
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in be found from Fig. 2; and 8 = 0.5 when (a/y2) = 

Slender body results are often useful in the form of ratios ex 
essing the interference effect of the body Since extensive in 
kists on tip controls fixed to the wing alone (reference 

rexample) we may define the ratio L(a@) /L(a@)) where L(a)) = 

2 - y,*)*?, the rolling moment produced by an 

\ wing-control system without a body This ratio is 
greater than one for all (a > 0.9, but tends to zero as (a/y2) = 
Setting \ ™=daasaq--e jy results because Lia —> () as 

hile L(a@)o ~ Oas (4 a)*!2 The physical explana 


is that the induced loading on the body is unable to contribute 
the carry-over onto the wing with no body present 


i favorable influence on rolling moment 
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Experiments on Chemical Kinetics in a 
Supersonic Nozzle* 
Peter P. Wegener, Jack E. Marte, and Carl Thiele 
et Propulsion Laboratory, California Institute of Technology 
Pa adena, Calit 


MODERN TECHNIQUES in fluid dynamics have recently 
pplied to determine the rate and mechanism of fast 
reactions \lso, aerodynamicists have become inter 
ted in flows with chemical reactions and, in particular, those 
vs Where relaxation effects become apparent \ new method 
measuring fast, homogeneous gas reactions in a steady-state 
personie nozzle? will be described, and observations will be 
here a reaction is lagging with respect to chemical 
uillibrium throughout a steady-state flow field rhe present 
speriment concerns the measurement of the rate of recombi 
a reaction which was 
sen because the properties of both species are well-known in 
mical equilibrium The reacting species are carried in small 
neentration in dry nitrogen, and this mixture is passed through 
nie nozzle The process under investigation may be 


superso 


V. + NO, —— N. + 2 NO 


represent the rate constants for the dissociation 


ere ky and kp 
| recombination reaction, respectively Nitrogen is stored in 
eservoir to permit continuous operation of a 2 by 3 cm. nozzle 
9 for about 1 hour rhe reacting species are carefully 

pared, metered, and added to the nitrogen, and the gas mix 
l screened 


chamber at 2 atm. total pressure and 400°K. temper 


reactant concentration, the NO, is fully disso 


ed into VO. in the supply. The ensuing expansion through 


This work presents one phase of research carried out at the Jet Pro 
n Laboratory, California Institute of Technology, under joint sponsor 
Department of the Army, Ordnance Corps, and the Department 

Air Force 
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Fic. | Concentration of nitrogen dioxide mixed with nitrogen 


in a supersonic nozzle 


the nozzle to VJ 1.9 proceeds practically unaffected by the 
recombination reaction caused by the isentropic temperature 
drop Phis is because the heat of reaction released in the flow 
is a negligible fraction of the stagnation enthalpy of the mixture 
at such low reactant concentration Phis fact was experi 
mentally verified, and, from a pressure calibration of the nozzle, 
the Mach Number and all thermodynamic parameters in the 
flow field are known 

The concentration of nitrogen dioxide as function of Mach 


Number, if no reaction takes place, can be computed ft 


Cyo,/(¢ ,1 + — 1)/2] 173} 2 
where the subscript 0 denotes supply conditions, : 100 
and Cya Dyoo/ RI Conversely, assuming that the recom 





bination occurs infinitely fast, taking advantage of the fact that 
the stagnation enthalpy of the mixture remains unchanged during 
the expansion and making minor simplifying assumptions, it can 
be shown? that the concentration of VO. as a function of Mach 


Number is given by 


where / l — 1)/2 
po(ati supply pressure 
1 K supply temperature 
¢ = volumetric fraction of VO», in supply 
A 6747, B 20.72 (from experimental results on the 
equilibrium constant as function of temper 


ture 


Local concentrations of VO. determined from 
are shown in Fig. 1 as a function of distance along the continu 
ously diverging nozzle In Eq. (2) the experimentally ob 
served Mach Number was used, and this result is, of course, 


f supply concentration rhe theoretical concen 


independent ( 
tration for chemical eouilibrium computed from Eq. (3) is shown 
for the single admixture of 1 per cent VO, by volume in the 
supply \lso shown in the figure are observed VO. concen 


trations determined from several independent experiments 
The experimental concentration was measured by optical ab 
sorption using a beam of near monochromatic light centering 


around 4,400A, a wavelength which is strongly absorbed by the 
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brown .VO., while N.O; is transparent This light beam was 
passed through the nozzle center plane at right angles to the flow, 
and it was received by a battery-driven photomultiplier tube 
The optical system could be moved along the flow during the 
experiment, and the incident and transmitted light intensities 
were determined. The concentration was finally computed from 
Beer's law by the use of a molar extinction coefficient previously 
determined with a similar optical system 

Fig. 1 shows the interesting experimental fact that, in this 
experiment where the flow undergoes a cooling rate of the order 
of 108 < d7/dt(°K./sec 10°, the net recombination reaction 
is not rapid enough to permit chemical equilibrium to be estab 
lished during the entire expansion in contrast to Feldmann’s' 
observations on air. Since all flow parameters are well-known 
at every Mach Number and since it is possible to observe concen 
tration accurately in a steady-state system, it is evident that 
the rate of the reaction may now be computed \sa first attempt 
to determine kp, it will be assumed that for Eq. (1) the net rate 


of concentration of VO. is given by 
dCyo,/dt QWpCv.o, Cv, — Wer Cro,? Cy } 


Without any further assumptions, Re or ky may be computed 
from the experimental data and Eq. (4), remembering that A, = 


ky/ kp, Where A, is the known equilibrium constant It was 
found that for 220 7(°K.) < 300 the recombination rate con 
stant is independent of temperature and in the range 0.007 

¢ 0.013, it is independent of concentration within the meas 
uring accuracy. The average value obtained from 28 experi- 
ments is kr = 2 X 10!'* (cm.® mole~? sec , and it is known 


within about 70 per cent. Finally, in this temperature range 
general agreement is found with the dissociation rate measured 
by Carrington and Davidson® who used a shock-tube technique 
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Approximate Solution for Slightly 
Yawed n-Power Bodies at 
Hypersonic Speeds 


Robert W. Truitt 
Professor of Aeronautical Engineering, Virginia Polytechnic 
Institute, Blacksburg, Va 

October 24, 1957 

_ OF THE results of references 1 and 2 suggests that, 
neglecting centrifugal forces, a reasonable approximation 

for the pressure distribution over a general slender body at finite 

hypersonic speeds and small angle of attack may be written as 


j= 2f1 +A[(1 + N)/4]} (cos asin @ — sin acos 6 sin B)? (1 
In Eq. (1),\ = (y¥ —D/(y + DY, N = (2f/Mao)?/A, y = ratio 
of specific heats, f = fineness ratio, 1J.< = Mach Number, e@ = 


angle of attack, 6 = angle between local tangent to a meridian and 
the longitudinal body axis, and 8 = angular position of a point 


on the body surface.! For the special case of Ma /f—~ =(y—>1 


Eq. (1) reduces to the expression for C, derived in refer 


where centrifugal forces and real-gas effects were negk 
the other hand, when a — (0) Eq. (1) reduces to the 
C, value given in reference 2 for slender conical surfac 


hypersonic speeds ‘N 2 or 5 


SOLUTION FOR ”2-POWER Bopy 


Perhaps the most suitable tvpe body of revoluti 
sonic speeds is the m-power body, of fineness ratio / 


meridian section shape given by? 


y/(d/2) = (x/L 

where L = length and d = diameter of the body The ad 
tages of the power body are twofold: a With vit 
centrifugal force effects, Eq. (2) is a reasonable appro iti 
the correct minimum-foredrag shape (xv = 3/4 when ceuntrifug 
forces are neglected*® and nm = 2/3 when centrifugal force 
considered?:*) ; and (b) the blunt nose of the »-power b 5 
propriate from the point of view of aerodynamic heating 

It is in order to discuss the applicability of Eq. (1) in the regi 
of the blunt nose of the m-power body—1.e., at 0 = 7/2. Fort 
case when a ~ Oand \Va/f—m + — 1), Eq. (1) isa satisfactor 
ipproximation since it reduces to the Newtonian impact val 
for finite hypersonic speeds, however, Eq. (1) overpredicts t 
stagnation pressure coefficient Nevertheless, the error in 1 
drag at finite hypersonic speeds, due to this overpredictior 
the nose, is quite small since these high pressures act on 
tively small area of large slope near the nose For the case 
small angle of attack and ‘ 2 or 3, Eq. (1) should be a reasor 
ible approximation for the #-power body pressure distributior 

Following the method of reference 1, Eq. (1) was appli 
tangent-cone approximation) to the m-power body type I 


expressions thus derived for the normal- and chord-for 


ficients are, respectiy ely, 


. 
Cy = (2Kn/f?") sin 2a f px dx n/2f 
e 
and 
C. = KCl 3/2) sin? alin 
% 

x2"! dx n/2 x ; si 

. 
where A = 1 +A(1 + N)/4 and x is in per cent of the 


ce ( 


body ¢ 


ameter i.e., Cy and C, are based on the maximum cross-sectiol 


area. In Eqs. (3) and (4), the maximum body diame 
taken as unity and thus the fineness ratio f = L. Thes 


for Eqs. (8) and (4) can be written in the general forms 
Cy = Kk, sin2a 
ind C. = K [ke + ks sin? a 


Table 1 contains the results of the integration of Eqs. (5 


for a fineness ratio f = 3 and a series of values of » 


ter W 


olutions 


ind (¢ 


Since the results in Table 1 do not take into account centrifuga 


force effects, the minimum-drag body corresponds to n = 3/4, 
previously found in reference 3 
CENTRIFUGAL-FORCE CORRECTIONS 
Cole,? taking into account centrifugal forces, derived an es 
pression for the drag coefficient of n-power bodies for the spec! 
case of zero angle of attack and V = 0 \ comparison of Col 
rABLeE | 


= 3, Without Centrifugal-Force Effects 


n k k R 
0.6 0.975 0.0540 Q 
2/3 0. 992 0.0474 0) 
0.7 0, 997 0. 0464 ) 
3/4 0.996 0. 0457 0 
0.8 0.994 0. 0467 0 


0.9 0.990 0.0503 0 

















For t 


itisfactor 


Ibutior 
S appli 
pe. | 


rece coef 


a) 


be vy di 
CCUOI 
ter W 
olutions 


and (f 


rifugal 
3/4, a 


in ex 
speci ] 


Cole's 


TABLE 2 
= 3, With Centrifugal-Force Effects 


C ko Ck 
( 0.0400 0. 680 
2 0.0372 0.730 
a9 0. 0374 0.751 
} 0.0390 0.787 
S 0.0416 0.828 
q 0.0478 0.878 


trifugal-foree correction, C, in Eq. (6): 


Table 


s the results of Table 1 corrected for centrifugal force 


ere, for the range of n(0.6 <n < 1), C = 0.382 + 0.7n 





nall angle of attack, there is no corresponding centrifugal 
ree correction for the normal-force coefficient, Eq. (5 


For the special case of a = 0°, K = 1 (Cra Cp), Eq. (7 
lable 2 reduce, of course, to Cole’s original values for drag 


efficient In any case, based on the present approximate 
ilvsis, Eq. (7) and Table 2 show that for any Mach Number 
\ < 2 or 3) the minimum-drag body corresponds to n = 2/3 
The final expressions for the lift and drag coefficient for the n 
ver body, fineness ratio f = 3, are 
C, = K[ki sin 2 a cosa — sina (C kz + Chk; sin*? a s 
Cp K |k; sin 2 @ sin a + cosa (Ck: + Ck; sin? a o) 
For the case of a = 0 
Cp = K Che 10 


tis seen that when 1/2 /f— K — 1), the error in neglect 


g the centrifugal-force correction is appreciable—i.e., the error 


sl5 per cent for » = 3/4 and over 20 per cent form = 2/3. On 
other hand, for finite Mach Numbers, the product AC ap 


iches unity and, therefore, the drag coefficient predicted by 


pact theory (C, = k») would tend to be correct (for a = 0 

rexample, when f = 3, the product KC 1 for the following 
mbinations of Mach Number and n: (1) » = 3/4 and N = 3 
K = 1.167) and (2) nm = 0.8and N = 2(K = 1.125). It would 


ir, therefore, that one should expect Newtonian impact 
ry predictions of pressure-drag coefficients (for a — (0°) to 
gree rather well with experiment for certain combinations of 
body shape and Mach Number—i.e., when KC l. 
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Buckling of a Simply Supported Rectangular 
Plate Under Compression Reacted by Shear 


Victor Weingarten 
tructures Research Engineer, Northrop Aircraft, Inc 


Hawthorne, Calif 
Ietober 28, 1957 


SUMMARY 


rhe energy method was used in order to find the critical buckling load of 
simply supported rectangular plate under compression reacted by shear 
rhe plate here considered is assumed to follow the classical small-deflection 
ory. One point was checked by the method of finite differences. Re 


Sare presented graphically 


those in Table 1 suggests the use of an approximate 


. 
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Fic. 1 Variation of critical buckling factor with aspect ratio 
SYMBOLS 
x rectangular coordinates 
t lateral dimensions of plate (1 
u transverse deflection of plate 
h thickness of plate 
D Eh3/12 (1 yu bending stiffne of plat 
a b/a 
PN Ver b?/x?2D 
ue Poisson's ratio 


positive integers 


NV cr critical compressive buckling load 

| strain energy of bending 

1 work produced by external forces 

\ Vy normal forces per unit length of sections of a plate perpendicu 


lar to x- and y-directions, respectively 


Viry shear force in direction of y axis per unit length of section of a 


plate perpendicular to x axis 


SOLUTION OF THE PROBLEM 


A SIMPLY SUPPORTED rectangular plate (Fig. 1), is considered, 
along whose sides x = 0 and x = } shear reacts a uniformly 
distributed compressive force applied at y = a. The compres- 


sive force at any interior point is assumed by the equation 


N = N rl V/a ] 


v 
The deflection of the buckled plate simply supported on all 
sides is taken in the form of the double trigonometric series 
w= > YD a, sin (mmy/a) sin (nrx/b (2 


ri 
ra) — 


m=l1n=1 


The energy method will be used to calculate the critical value 
of the compressive force V,,. The strain energy of bending is 


given by the equation 


eo eh 
(1/2)D f, /, { [(O%w/dx?) + (0%w/dy?) |? 


(0*w /Ox? O*w/Oxdy)?]j dy dx (3 


AV = 


2 (1 bh) [(d%w/dx? 


Substituting Eq. (2) into Eq. (3) and integrating we obtain 


AV = (rtab/8)D D> do ann? [(m? 


m ln l 


The work done by external forces during buckling of the plate is 


eo eh 
AT = (1/2 / { N.r y/a)(Ow/Oy)? dy dx 5 
J0 J0 ’ 


Substituting Eq. (2) into Eq. (5) and noting that 


, 
f. y cos (r rv/b) cos (sry/b) dy = 67/4 for 
gt * 7 - i 


eh 
| y cos (rrv/h) cos (srv/b) dy = 0 lorr $= 
wv V 


r +s is an even number 


b 
y COS (’ryv b cos (S7rv/0 dy = 
0” . . 2 
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2+ sisan odd number. We obtain the following equation: 


AT =(Nwb/4a)( >> Dd a, 


2 (m/a)? (a?/4) 


n=1m=1 
> > ) AmnApn( Ma/a) (pa/a) X 
m=ln=1 —p 
$ (972 /_2 2 2 ,2 p2)2]) | . 
ylca*/n*) (Cm? + p*)/(m p { (6 


m ¥ p, m+ pis an odd number 
In order to simplify computation, the general expression, Eq 
(2), isreduced to 


w = sin (2x/b) B dy 


, sin (m7y/a) (7 
m=! 

which is equivalent to the assumption that the plate buckles into 

a half wave in the x direction. Letting the strain energy of bend- 

ing equal the work produced by the external forces in conjunction 


with Eq. (7) we obtain 


am | [(b?/a*) m? + 1)? Nor b4/2 a? 22D) m?} 
(2 N., b4/a? z4D) 5 Ap) mp |(m? + p*)/(m? p* 2] = (8) 
m ~*~ p, m+ pisan odd number 


By letting a = b/a and Xd = (N,, 6?/%?D) and writing Eq. (8) 


for three equations we obtain 


[(a? + 1)? — (a?/2) A] an + [(20a2/9%?)A Jan = 0 


[(20a2/92?)d Jan + [(4a? + 1)? — 2a*rA Jan 4 
[(60a?/9z2)\ ] an = 0 


[(60a2 Ox? Ala t [(9a2 + 1)? (9/2 ad | a =) 


Equating the determinate of these equations to zero, an equation 
for calculating the three term approximation is obtained which 
has been found to be sufficiently accurate for this problem 

A plot of @ versus J is found in Fig. 1 extrapolated to the limit 
ing cases: 


Case l: a— Approaches a plate under pure bending 


» 


A = 23.9 (reference 1) 
Case 2: a — 0 Approaches a plate under uniform compression 
\ = 4 (reference 1) 


For an aspect ratio of a = 2 a finite difference solution was 


performed using 21 interior points. The biharmonic Eq. (9 


(04g /Ox4) + 2 (046/Ox*0y?) + (046/04) = O (reference 2) (9) 


was solved in finite difference form using the appropriate bound- 
ary conditions in order to determine the stress distribution 


throughout the plate. The buckling equation,! 


(O4w/Oxt) + 2(04w/Ox*0y?) + (O4wW/ 044) = 
(1/D) [Nx(07w/dx?) + N,(0?2w/dy?) + 2 Nyz,(O%w/Oxdy | (10 


was then solved in finite difference form which determined the 
critical buckling load. A comparison of the finite difference 
solution with the energy solution showed a difference of 15 per 
cent. The finite difference solution has been found, in previous 
papers, to be a lower limit of the critical buckling load.*: 4 

The energy solution always gives an answer above the true 
buckling load. The true solution is therefore somewhere be 
tween the two. 

Since the energy solution assumes a stress distribution to equal 
N.- y/a and the finite difference techniques calculate the stress 
distribution, this might further account for the 15 per cent dif- 
ference of the two answers. This error will become small for 
lower values of @ since the stress distribution will approach the 


jinear variation assumption 
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A Note on the Expressions for the Local 
Heat-Transfer Coefficient 


Richard R. Gold 

Member of Technical Staff, Aerodynamics Department 
Hughes Aircraft Company, Culver City, Calif 

October 28, 1957 


1 


expressions for the local heat-transfer coefficients in lamin 


F yews DEAL Of effort has been put forth on establishing { 
and turbulent flow. Several surveys on the subject prest 
an extensive picture of the problem along with the assumpti 
made in the various theories for incompressible as well as cor 
pressible flow The remarks of this note will be restricted tot 
ease of an isothermal flat plate at zero angle of attack I 
standard transformations? provide the corresponding expressior 
for the cone 

For the case of laminar incompressible flow, the following gen. 


erally accepted result of Pohlhausen is used: 
Nu, = 0.332 Re!!2Pr! 


Frequently used results for laminar compressible flow are thos 
of Chapman and Rubesin, and Van Driest. 

Using Blasius’ equation for the turbulent skin friction 
Reynolds’ analogy, where S = Reynolds’ analogy factor P 
the local Nusselt number for turbulent incompressible flow 
given by 

Nuzy = 0.0296 Re*4Pr! 2 


\ widely used expression for compressible flow is due to \ 
Driest. 

Various authors have attempted to apply the incompressibl 
heat-transfer relations to compressible flow by evaluating t 
fluid properties at some reference temperature. The success 
this concept has been illustrated by Eckert? for both laminar 
turbulent flow 

In keeping with this idea consider* the variation of the lov 
heat-transfer coefficient as indicated through Eq l I 


equation reduces to 


hy = 0.332 (pu/Rx)s'? F(1 
where all the temperature-dependent quantities have been « 
lected into the single term 
F(T) = (k Pr*!?/V pT) (B.t.u./sec. ft. °F.) > 
it Ib. sec R 


The subscript 6 refers to conditions at the outer edge oi 
boundary layer while the subscripts 1 and 7 refer to laminar 





turbulent flow Using the results of the standard gas tables 


along with those of reference 4 to evaluate F(7) independent 
it can be shown that F(Z) is approximately constant in tl 
range 100 < T < 2000°R 


0.000256, in which case Eq. (3) reduces to 


This approximate value is F(T) = 


hy, = 2.0579 X 10-6 (pu/x)s! B.t.u./ft.2 °F. sec 


This result is within the accuracy maintained in the analysis 
of Eckert, for example, where the validity of using the constant- 
property solution for the compressible case is shown. As a ft 
sult, one can observe that, within the accuracy of our knowledge 


wnt 


the variation of the fluid properties with temperature, the assum] 


* This approach was used by A. E. Bryson, C. P. Stensgaard, Jr., 4 


J. M. Bozajian as early as 1952 based on reference 5 
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Journa n that #7, is independent of the choice of reference temperature, Eckert, E. R. G., Survey on Heat Transfer at High Sp WADC TR 
p. 384 I . ‘ ; 54-7 ) 95 
ndeed that it isactually independent of temperature, isreasonable. 1-70, April, 1954. 
vs . : ai : : ; Decoursin, D. G., Bradfield, W. S., and Sheppard, J. J ferodynamt 
25S For ease of use, the resultant simplified expression for h;, given by Kistius aut Gant Tranter Phsomes t Mach Numt 
57-164 Eq. (5), can be reduced to the form of a nomograph by standard WADC TR 53-379, February, 1954 
jods. It is interesting to note that if the frequently used Tribus, M., and Boelter, L. M. K., An [nvestigats 1 aft Heaters II 
cpressions for the variation of viscosity and conductivity with Properties of Gases, NACA, A.R.R. W-9, 1942 
P . ; - Johnson, H. A., Rubesin, M. W., Sauer, F. M., Slack, E. G., Possner 
nperature, namely Pr = constant, 4 ~ A ~ T, are applied to Re i aa Er ee ee: ‘iu Thermal Characiosinix af Wa 
the function F one would obtain: F(7 = constant Ex Speed Ai aft, AAF TR 5632, September, 1947 
I rimental verification of the above has been obtained from ref 
nce 3 where the laminar expressions of Pohlhausen, Chapman 
> - . ; « 
ind Rubesin, and Eq. (5) are compared with test data. Further 
erification has been obtained from existing flight-test data on 
nes of various sizes covering a wide range of Mach Number, On the Generalized Prandtl Relation 
ids’ Number, and ratio of wall temperature to temperature 
utside the boundary layer. The results of this correlation indi Rudi S. On 
‘ pee udi > dng 
that Eq. (5) predicts the test data reasonably well Phe > 
: is I sate Assistant Professor, Department of Aeronautical Engineering 
ae results ol Eq 5) are in excellent agreement with the values pre University of Michigan. Ann Arbor. Mich 
ishing { : oe 


ed by Pohlhausen’s expression, as would be expected, and are October 98, 1957 
in} lamin ‘ 


in good agreement with the values predicted by the expression 
1g ' I - ] 


eae f Ch —— Rubesin, oon Van I rriest The data includ co R A FOUR-DIMENSIONAL Euclidean space-time mani 
mytantes: 1 uppeectatne amount of scatter at times in addition to several fold, E, The time catlable will be deacted tw &. and th 
eltend ) mats lying in the transition region, which accounts for some of space variables by x*, A 1,2,3. In particular, we assume that 
ck. TI disparity with theory \s a result of the above discussion, the space variables are Cartesian and define Euclidean three 


pressior tn be stated that A, is independent of temperature for all spaces E; in E, 

| ae practic urp ‘ he result: si ified expression for ; ; . ‘ , 

practical purposes and the resultant simplified expression for The tensor formulation of the hydrodynamical equations in 

can be co : repres d by asi omogri . . . 

— in be « nveniently represented Dy a imple nom gt iph Gina scactaie wechein tide de Geren 
wii \ similar reasoning can be applied to the case of turbulent 


mpressible flow in which case the results serve to justify the 0/ot 














pv» + Va pvr», + po, ] 
ilvsis of Eckert and present a calculation simplifying expres 
‘ : ; yh Op/Ot) + Vy pr" 0) 2 
sion based on his equations We will use the analysis of Eckert 
ire t mae = om ss 2 > /é ‘ i 
utlined in reference 2 Phe experimental calculations examined 0/Ot) } pl (g?/2) + e]§ + Vyuipr"*l(g?/2) + h)j () } 
) ithe present investigation have generally been in the range Reg F 
ion an ‘ ; ; 6” is the usual notation for the Kronecker delta, and p, ¢ et 
Pp < 40 X 10° in which case the constant value S = 1.18 suggested a ee : ' ” ikl idiots 
a 1 the literature is used \ theoretical correlation with experi and Af denote the density, component of velocity, magnitude of 
ow aa : ' 3 . ' : i | 
mental data using Blasius’ incompressible expression gave good the velocity, internal energy, pressure, and enthalpy, respectively 
results in this range. Asa result, this expression is used instead Eq. (1) is equivalent to the Eulerian equations of motion for 
2 f the several expressions suggested by Eckert for turbulent compressible nonviscous fluids in the absence of body forces; Eq 
67 incompressible flow In order to simplify the caleulation of the 2) is the equation of continuity; and Eq. (3) is the energy rela 
o \ - ; 
it-transfer quantities from Eckert’s analysis, the temperature tion 
as ependent terms have once again been collected into a single The hypersurface (or lower-dimensional manifold) in space 
ing t term in the expression for the local turbulent heat-transfer time along which discontinuities occur will be denoted by 
ne efficient. The resultant expression is given by 
one ; : @ (t, x” k where the /& are constants } 
nar hy = 0.0349 [(pux)s°8/xR°8|G(T* / 
, : . " - = : rE ‘onsists of o on ation (7 ) het lis 
1 k Pr/(pl Btu. /sec. it. °F 6 If Eq. (4) consists of only one equati , then the di 
lo ar) «© 1 j ifo 1 Ss Ss i ] 1 oO 
1 soca ft.2/lb. sec. °R 3} continuity manifold defines a hypersurface If the m con 
This sists of more than one equation, then the equations define a mani 
vhere 7* is Eckert’s reference temperature. Using the afore fold of lower dimension. In any case, the vector fields for the 
mentioned variation as a guide, namely, 4» ~ A ~ 7, one ob _ j j — ' 
, “0 . . . various types of 7, 0¢/O0/, Om@/Ox", determine vectors normal to 
tains G(7 ~ 7 6 and applving the previously mentioned = ‘I cif: ’ wh 3 
, btaj ; E the discontinuity manifold Furthermore, the unit normal 
en lata one obtains, with adequate accuracy, ages ; ° 
vectors of this manifold are determined by 
G(T*) =z 0.113 7" 6 150 < T < 2000°R, 
n, = (Op/Ol Og /Ot)? + g™#*(0~/Ox") (0G/ Ox" 
|} suchthat 5 
: ~ ny, = (Q@/Odx° Op /Ot)? + g*#(0p/Ox*) (0G /Ox# 
hr = 10.25 XK 107*[(pux)5°*/x7T™- ( 
of the a ; where g®? is the metric tensor 
Eq. (7) can also be conveniently represented by a nomograph a ; 
ir é pe : rae : In the case where the discontinuity manifold is a shock wave 
ve Various experimental results of flight-test data which were ag . : 
tables ; the ’ : a - A the values of p, vy, ¢, p, and # and their space-time derivatives are 
lescribed above indicate that Eq. (7) predicts the test data gen . - . 
lent ‘ : , ‘ assumed to be continuous, while in crossing the discontinuity 
erally within 10 per cent, and, in several instances, this figure ; 
n the a age Sans manifold these values are discontinuous 
| rises to 15-20 per cent It is of interest to note that the results 
= > ae = : e . ’ Ty } 1] of o7'. across the sho lis 
f Eq. (7) are conservative. For comparison purposes the well Let pv, denote the jump in the value of pz) act the shock di 
known ¢ xpression of Van Driest was used to predict the same test continuity, pv4“7, the jump in the value of pr“zy, et Phen by 
5 lata. The results of this correlation were similar to those ob means of the integral forms of Eqs. (1), (2), and (3), the following 
l tained using Eq. (7) except that the results of Van Driest were jump relations corresponding to Eqs. (1), (2 may be 
alysis 2 - = . : 
nconservative. Eq. (7) generally overpredicted the test data obtained 
stal : . : . ; ° 
Van Driest’s expression gave values which were generally 
re : \ . . ° —— 4 \f lk 4 a yoy + 84H ( ‘ 
1 si \s a result, the analysis of Eckert seems justified Ob/ OF) pt 09 / Cx " f 
ge 
ump- : 06/d1)5 + (d/dx") pr 0 - 
, REFERENCES 
Kaye, J Surve Fricti Coe fictent Re ery Fa ind Heat Od Ol pg? /2 + pe + Od Ox” pq*7 /2 - pri’ U Ss 
- 1 ( fi { Su pe ic Flow, Journal of the Aeronautical Sci 


ences. Vol. 21, No. 2, pp. 117-129, February, 1954 Let the subscripts 1 and 2 indicate values on the sides of the 
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shock wave, then 
aB = (aB)y — (@B)i = (a2 — ay) (Bs — Bi) F 


a) P2 T @ 3; = 2a8; 


Upon simplifying the right hand side of the above, the jump may 


be expressed in the form 


Hence, pUx pur~n T PrTK TT Vp Pp 
vy,pu" ete 
Expanding Eq. (6) in this manner the following equation is 
obtained: 
(Og/Ol pv’, tT pilr tT PU t 
\ 

Od / Ox" pv", 4+ pv" ivy + pv* vy ro. Od Ox 0 
Assume now that p = co where c? is the ‘velocity of sound.” 
, the above equation may be reduced to 


Then using Eq. (7 


pity |(Om/OL) + 7,4(O0G/Ox*) + c7p(Od ox” = 


Solving for #, 
p= —} pivrl (0d O/) + 74(0p/Ox*)] /c2(O@ ox” } 9) 
Expanding Eq. (7) the following equation is obtained: 


O¢/Ox") (pv* + pyv* + pv") = 0 


(Od ol ” ae a 


or pl (dp ol + (p* + 4 Od ox ] = piv" Od Ox’ = () 
Substituting the value of from Eq. (9) into the last equation 
vields 


pib| (Od Ot) + v,4(0d/Ox*)|1(06/Ot) + (F* + 7,4) (0G/Ox4 


c? (0¢/0x*) 
pid#(Ob/Oa*) = O 
or Py O¢/ Ol + vv Op Ox™) |, (O@/Ot) + rtd Od Ox? asi 
cv! (Op/Ox") (06/Ox*) = 0 (10) 


Since the tangential component of velocity is continuous across 


a shock discontinuity, 


F, = F4(Oh/ON#) (OG /ON>) /g%4 (OH/OK®) (OG/OKF 
Thus Ea. (10) becomes 
By} [(Op/Ot) + v*(OG/Ox® O¢/Odt) + v8(0¢/dx*)]. — 


t 


0¢/dx’)} = 0 


2,08 a 
c#g®" (Op /Ox 


Since by hypothesis 3, # 0, 


[(O0p/Ct) + 7%(0p/Ox® 0¢/Ot) + v79(0b/Ox?)]2 — 
c2g (Op ox” (Od Ox?) () 11 


Expressing this equation in terms of m;, my, as defined by Eqs 
(5), vields the following equation: 


(ny + vny), (Me + vYNy)2 = Cc? 


This is the generalized Prandtl relation. 
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Further Comments on Diffuser Instability 


W. C. Randels 

Missiles Systems Division, Lockheed Aircraft Corp., 
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I HAVE FOLLOWED, with much pleasure, the spirited polemics of 
Dailey! and Trimpi? on the subject of diffuser ‘‘buzz.” 


ABRONAUVTICAL 


SCIENCES MARCH, 1958 


In view of the obvious satisfaction derived by the antagonist 
it seems almost cruel to point out that the whole thing cap }, 
explained by a classical theorem on compressor surge which mg 
be found in standard textbooks*~ on rotating machinery 
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A Test of the Uniqueness of Solutions for 
Problems of Nonsteady Flow Under Given 
Boundary Conditions* 


J. Altenhoft 
Research Staff, General Motors Corp., 
November 4, 1957 


Detroit, Mich 


2 STEADY STATE Which is approached by a flow system, fron 
given space-like initial states and under given time-like 
periodic boundary conditions, is generally postulated, for lack of 
a rigorous proof,’ to be uniquely determined by the boundan 
conditions and independent of the initial state—i.e., independent 
of the transient processes through which it is established. Even 
in the absence of viscosity, when the approach to periodicity can 
be expected to be asymptotic, the existence of a uniquely deter 
mined final steady state, regarded in this case as a limiting state 
is usually taken for granted 

These postulates stem from physical intuition and have als 
derived at least qualitative support from the experimental ob 
servation of the behavior of periodic-flow devices under con 
trolled boundary conditions 

This note concerns an investigation, the results of which pri 
vide further evidence of the uniqueness of the final steady stat 
under given boundary conditions. The object of this investi 
gation was to determine, by the method of characteristics,? the 
final steady state or limiting state which would be approached by 
a particularly simple, inviscid, and isentropic flow system under 


given boundary conditions from various randomly chosen initial 


states. 
The system chosen for the analysis consisted of a tube filled 


with and surrounded by air, open at one end and closed at the 
other end by a piston oscillating with simple harmonic motion 





* The author wishes to thank Dr. J. V. Foa, Rensselaer Polytechnic 
Institute and Dr. R. Davies, Research Staff, General Motors Corp. for their 


guidance and assistance in preparing this note 
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System used for analysis. 
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ntagonists ol2 TABLE 1 
ing can hy : er 
wl Path Cycle at Which Steady 
= tae O10 Run State Was Reached 
, l 1] 
- 14 
008 3 16 
\erona yO -O—0 -0-0-0- = = = -0 4 11 
o i> 











imn 
Pe 004 nificant changes in this distribution from one cvele to the next was 
Pitma; interpreted in each case as indicating the establishment of the 
cadets 002 final steady state _— validity of this interpretation Was veri 
fied in each case by carrying the computation further, for several 
ae 000 eveles, to check for the possibility of a subsequent divergence of 
(@) 2 o 6 8 lO 12 4 16 18 the solution (which was, however, never encountered lable 1 
Ww York CYCLE NUMBER gives the number of cycles required to reach the steady-state 
re 2. Values of A at 0.1 at the end of each cycle (Run | condition 
‘ One can illustrate the way in which the solution converges to 
Fig. 1 fhe piston oscillation frequency, the ratio of piston — Steady-state value by means of a curve as shown in Fig. 2 
stroke to tube length, and the temperature of the air in the tube which is for the initial conditions of Run | The curve shows 
sere so chosen as to avoid the formation of shock waves, and the the value of A at the end of each cycle and represents the condi 
tube diameter was assumed to be large enough to make it per tions at the location ¢ 0.1 along the tubs cms waMue OF 
1s for missible to neglect boundary-layer effects. The open end of Berke — the variation of 1 was greatest at the piston 
Given e tube was assumed to be flared so that inflow would be nearly = . “ ; _ ' ; eae 
‘seutropi For the “standard” initial conditions (Run 1), the m iximum 
In all the cases analyzed, the piston oscillation frequency was eiierenct senna ote Aiope se oo ve ee 
ssumed to be 4000 eyeles per min. and the ratio of piston stroke It was arbitrarily decided that the solution had come sufficiently 
tube length was assumed to be 0.25. The following initial close to the steady-state value when the change during one cycle 
sill tas aati oidnaliliewaal was less than 2 per cent of O.OLO77 his value was then used as 
the criterion for convergence for all runs 
1) A 1.0, C v for 05€3 10 The value of A along the length of the tube at the end of vari 
n, fron 2) A 101, V=0 for 0OS€510 ous numbers of evecles is shown in Figs. 3 and 4 Fig. 3 illus 
me-like 3 0.989, 0 for0 S&S 1.0 trates the convergence for Run 1 in which the steady-state con 
lack of 4 1.0 for o 3 $ =| 0, l = 0.001 ie s dition was reached by the eleventh cvcle Fig t shows the con 
undar\ 5) A 1.01 for 0.5 — = 1.0, A = 0.989 for 0 S & 0.5, vergence for Run 6, the case that showed the most erratic vari 
endent Ofor0 S§S 1.0 : ation in A Although the solutions for these two sets of initial 
Even 6) Same A as (5) but U = 0.025 for 0 S & S 1.0 conditions differ greatly during the first few cvcles, they do 
ty can where £ v/Lo, and A and Ul’ are the ratios of the local speed of eventually converge to the same steady-state value Likewise 
deter sound and of the local particle velocity, respectively, to a refer for the other cases the first few cycles may differ consider 
State ence speed of sound ably, but eventually the solution approaches a steady-state value 
The calculations were carried out on an IBM model 704 that is independent of the initial conditions 
e als ligital computer In summary, within the accuracy of these computations, the 
al ob The distribution of state and flow parameters along the tube final steady state is the same for all cases considered, independent 
~ Con was calculated at the end of each cvcle, and the absence of sig of the transient processes through which it is established 
| pro 
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Fic. 3. Variation of A along the tube at the end of various Fic. 4. Variation of A along the tube at the end of 
cycles (Run 1) various cycles (Run 6 
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REFERENCES where E is the elastic modulus and the skin is taken the ¢ 
R. Courant and Kk. O. Friedrichs, Supersonic Flow and Shock Wave pression member ; 
p. 367; Interscience Publishers, Inc., New York, 1948 For given values of C, in Eq. (3) the skin stress f,/ F.,, ca 
2G. Rudinger, Wave Diagrams for Nonsteady Flow in Ducts; D. Van graphed against the elastic skin stress (f F 2); Fig. 1 ows? 
Nostrand Company, New York, 1955 aan Mie ae ease Ok a 00 if C, is larger than valu rat 
the roles of the skin and stringer can be reversed to give 
; = 
+ Cy 
Note that the stress-strain curve in Eq. (2) is represente 
. . . Fig. 1 by C, = 3/7 with (f/F = eE/F If during the apy 
Inelastic Combined Thermal and Applied eee ee ee oe — Py 
Stresses in Skin-Stringer Aircraft Structure cation of the temperatures (J.) should decrease, then be 
most materials unload elastically the unloading curves Fig = 
SE Guewood will be parallel to the C; = Ocurve. 
° J we . ° ° . ° , . 
Research Professor, Air Force Institute of Technolog For the case in which the skin-stringer element is subjected t 
, ! Cc CCU SY, = x am P 
Wright-Patterson AFB, Ohio tension or compression force Q, Eq. (1) becomes y 
er 6,19 . or / 
November 6, 1957 f, = [—E,(a,T, — a7;) + (E,/E;) (Q/A: : 
1 + (A,E,/A,E, 
N CHAPTER 7 of reference 1, procedures involving strain equa f, = (O/A, (A./A,)f 
tions and trial-and-error calculations are used to determine ;. 
inelastic thermal stresses. If the stress-strain curves are repre By using Eq. (2) the inelastic stresses for Eq. (4) can be writter 4 
sented by approximate equations, it is possible to obtain the in the form 
elastic thermal stresses directly from charts, and the procedures f iat, ~ mt, + BAB/R) C/A 
for combined applied and thermal stresses can be considerably : = — : eae eae 
a " F F.,, (1+ (A,E,/A,Ey)] (+ Gf. /Feys)” | 
simplified 
" ° ° ° ° ° f , = , sie , , f , 
Consider a skin-stringer element in which the skin of area A, o/ Fy Q/ AF iy A,/A I F iy Py 2 
is at the uniform temperature 7’, above datum and the stringer By = 1+ 6/7) (4/F,3)" 
of area A, is at uniform temperature 7, above datum. From Wb A.E. (f./F,..)%7! 
Eq 3-1) of reference is the thermal stresses in the skin and 1 + (A.E./A,E, A.E, (f./F 
stringer are (the skin and stringer restrain each other and there 
is no bowing or buckling where m, and m, are the values of » for the skin and stringer 
: Tee , respectively 
j = ~ EsasT, — aT, + (A,E./A,E,)) 5 1 Eq. (5) can be solved for the inelastic stresses by a successiy 
_ f . . . s+ a Py 
Ay/A, approximation process using Fig. | rhe procedure ts as foll L 
where a@ is coefficient of thermal expansion, E is secant modulus, l Calculate (f,/F,, and (f,/F;,,)e. from Eq. (4 
and f is thermal stress 2) Caleulate B,; and C; from Eq. (5) using the elastic values 
\n approximate equation for the stress-strain curves for vari 3) Calculate a corrected (f,/F,,,)e using B; but taking ( 
ous materials is the Ramberg-Osgood equation (see section 6-3 of in Eq. (5 
reference | for discussion that » = 10 gives a good approximation +) Use this corrected (f,/F,,,.)e. and the calculated C, to re 
at various temperatures for various materials the inelastic f,/F,,,o0n Fig. 1. This assumes that » 10. Fora . 
PEO toe ne different #,, curves similar to Fig. 1 can be constructed fron \ 
cE/F f/F 1 + (3/7) (f/Fo,)” 2 t7 (2 | 
/ / : y Eq. (3 Swit 
where F,, is the compressive yield stress and E is the elastic 5) With this value of f,/F,,, calculate f,/ F,,, from Eq. (5 
modulus. Assume that the stress-strain curve in tension has 6) Use these inelastic stresses to calculate a new B, and ¢ 
the same shape as in compression and take the value of » for the 7) Repeat steps 3, 4, and 5 nd 1 
skin and stringer stress-strain curves to be the same. Then by (8) Repeat the process until the stresses repeat within t vith 
using the secant modulus from Eq. (2), Eq. (1) becomes desired accuracy rl 
. Note that this process can be used for thermal stresses 1 
f /F f /F a Pe a when m, + n, by taking B 0 Note also that if unloading 
) Fae f l ( I : 
f /F th le T. = wal occurs during the temperature and applied load cycles, B l Prof 
F...{1 + (A.E,/A,Bs)]! 3 for the stringer unloading and C, = 0 for the skin unloading eralt: 
C 3/7) [1 + (EB,/E,) (A,/A P/F ; 
e 1E/A.E REFERENCES 
Gatewood, B. E., Thermal Stresse McGraw-Hill Book Company, In 
New York, 1957 
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10 
It 
‘ad A Special Case of Swirling Viscous Flow 
L, . yp! 
“ye Dennis G. Mabey 
» ine 
06 Royal Aircraft Establishment, Bedford, England 
November 7, 1957 < 
OF 1 
/ SYMBOLS 
QZ f+ } } } } j } } } r,9,@ spherical polar co-ordinates 
U,V,u velocities in 7, 6, @ directions 
F mn cos 6 I 
ot I kinematic viscosity the 
0 OF 08 /2 4. 20 oF es 32 R characteristic Reynolds Number l ! constant refe 
Be ) 
| tous ” | Subscript 





a magnitude on jet axis 
Fic. 1 Inelastic thermal stresses The prime indicates derivatives with respect to u 
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PHEORY 


N EXACT SOLUTION has been found to the equations of incom 
A pressible viscous flow for axially symmetric motion with 
Che swirl velocity W is of the form 


W = —of(u)/r(l — p?)!/? l 


nd the stream function is ¥ yrf(m) as in the earlier solutions 


vithout swirl.? 
then give 


rhe three equations of motion (reference 5, Eq. 45 


ff? = pe fr 9 
Prof. H. B. Squire remarked to the author that 7 may be gen 
eralized in the form W’ vh(p)/r(l — p*)'* Phe equations 
f motion then give 
/ 1 — ph 
Z2i(1 — pw — i l— yp _ 1 — p 
The simple case considered in Eq. (2) when A(p f(u) corre 


sponds to We = V 

It can be shown from the variation of the pitot pressure that 
}. (2) can only describe flows away from the origin, for which 

<0. No wall boundary condition (1 V W 0) can be 
pplied, except in the trivial case when the fluid is quiescent,' 
since the solution represents a jet mixing process. The relevant 
boundary condition is the characteristic Reynolds Number of the 
et R, and, to illustrate the two types of solution, the appropriate 


undary-layer solutions are now given 


TYPES OF SOLUTION 


Formally the Type 1 solution represents a source vortex along 


the 1 axis and a force applied at the origin (as discussed in 
reference 2 The boundary conditions are 

f(—1) = 0, f’(—1) = (ur/s , —R 
When R >> 1, then f tends toa limit when w= — 1. The bound 
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irv-layer solution is 
_ 1k ] T pw ‘+ f | T pw 


and the potential flow solution is 1, which give 
streamlines in the inflow 

The Type 2 solution represents a radial jet in the plane u 
combined with a pair of source-vortices along the u +] axes 


The boundary conditions are 


f(O 0, f’(O ur/1 = 
When R > 1, then f tends to a limit when p< 1 Phe boun 
ary-layer solution is 
—(2R)'/* tanh [(R/2)'/* p 

Both types of solution can be obtained numerically for R 
without making the boundary-layer approximations Note that 
for an air jet with a speed of 300 ft./sec., 1 ft. from the orifics 
R 2.108 

ryPES OF FLow 
The Type 1 and Type 2 solutions are radically different. The 


Type 1 solution represents a small jet along the pu —1 axis witl 


swirl in the inflow about the jet axis Fig. 1 shows the stream 
lines of the motion for R 2.10 Phe angular width of a jet 
may be specified as the angle @ between the 
O.5U,r,/v, and the jet axis For 


1/2 


half-velocity cone 
mw = mw, on which L’(p))re/s 
the Type 1 solution @ 1.82/(R 

normal to the jet axis, the maximum swirl velocity is 0.71»(R)!/ 


\long a line ry 


rq (see Fig. 1 
The Type 2 solution represents a radial jet in the plane yu 
+1) which is equal and opposite 


with swirl (about the axes yu 
on both sides of the jet Fig. 2 shows the 


For the Type 2 solution the half-velocity 


streamlines of the 


motion for R = 2.10% 
= 1.25/(R)'/* 


ixis, the swirl velocity 


cone is 6 \long a line r(1 — y?)!/? ry normal 
to the jet 


y(2R)1/- r, (see Fig. 2 


attains a constant value of 
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Phe streamlines close to the jet axis diverge for the Type | 


solution and are parallel for the Type 2 solution as shown in 
Figs. | and 2 
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Subsonic Lift Response to Penetration of a 
Sharp-Edged Gust Moving at Supersonic 
Speed * 


E. Krasnoff 


4 


Staff Engineer, Allied Research Associates, Inc., Boston, Mass 


November 20, 1957 


INTRODUCTION 


of a 


been stimulated be 


pee REST in the problem of the penetration (by a wing 


moving sharp-edged gust has recently 
cause of its application to the blast wave penetration problem 
The relation of a moving sharp-edged gust to a blast wave is 
discussed in reference 1 which includes caleulations of lift re 
sponse functions for the following cases: airfoils in both incom 
pressible and supersonic flow; wide delta and rectangular wings 
in supersonic flow; and very narrow delta wings in both com 


pressible and incompressible flow. However, the only lift re- 


sponse functions that have been calculated for subsonic, com 
pressible, two-dimensional flow are those associated with a 
fhe author wishes to express his appreciation to the Aircraft Labo 


ratory, Wright Air Development Center, USAF, for permission to publish 


this note 
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Fic. 1. Representation of the airfoil-gust penetration problem 


in the x, ¢ plane 
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Fic. 2. Typical lift response functions for an airfoil M 


Number of 0.8 


sudden change in airfoil angle of attack (infinite gust f1 


locity) Results of these calculations are included in referer 


2 through 6 

The purpose of this note is to present lift response functior 
for an airfoil moving at subsonic speed corresponding to penetr 
tion of finite supersonic spr 


a sharp-edged gust moving at 
Linearized airfoil theory is employed and calculations are earrit 
out for periods of time somewhat exceeding the time required f 
the airfoil to be completely immersed in the uniform gust fi 

Both 


orientations are considered and the general results are spt 


head-on and tail-on (gust front overtaking the airfoi 


to the case of infinite gust front velocity 


DISCUSSION 


The representation of an airfoil moving at subsonic spec 


being overtaken from the rear by a sharp-edged gust moving 
supersonic speed is depicted in the 
Here the 
nondimensionalized by the chord of the airfoil, ¢ 
T, by ao/c, is the sound speed in the homentropic fl 

field That is, X/c = x and aol /e t Vand J/ re t 

Mach Numbers of the airfoil and gust front respectively, taker 


plane of the atr 
m8 s beet 


, and the tin 


v, ¢ plane 
foil) shown in Fig. 1 space dimension, 


where a 


positive in the positive x direction 
Calculation of the Lift Response Functions 


It is well-known (see, e.g., reference 7) that the boundary v 


problem posed by compressible nonstationary airfoil theory 


analogous to that of supersonic finite wing theory, the governing 
partial differential equation in each case being the wave equa 
tion. Thus the techniques of supersonic wing theory (in particu 
lar those of references 8 and 9) were employed to compute t 
pressure coefficients in the five regions indicated in Fig. 1. TI 
lift coefficient was then obtained as a function of time by inte 
grating the pressure coefficients over the chord of the airfoil 
computations are discussed in some detail in reference 7 


The following notation is introduced: 


where Cr, = 2ra,/V 1 72, the steady-state valve 


lift coefficient, C,, attained as ¢ approaches infinity 
sults are presented below 


For the tail-on case 


0<t<S 1M, —- MW: 6 =8t VO i 2 
ind, during the analogous time interval for the head-on cas 
7, « l 

iY <2 < l V7 a VV 





Indic 
Tl 





onse functions during the remaining time inter\ 
mputations have been performed are presented be 
I pper (lower) sign refers to the tail-on (head-on 
In specifying the time interval, the upper (lowe: 


to the tail-on (head-on) orientation 


A UV (iM VW) (M, + 1)/|M,| 1 + M 65 


Indicial Lift Response Functions 


rhe lift response functions corresponding to a sudden change 
gle of attack (infinite gust front velocity) can be obtained 
Eq 1) through (6) by performing a limiting process in a 
raightforward manner. In view of Eq. (2) or (3), the first timé 


| shrinks to zero and the lift instantly becomes 


Se¢ From Eqs. (4), (5), and (6) there results 


1/1 + M @ = 6,(1 l VW)t 9 





Numerical Results 
Representative numerical results are presented in Fig. 2 for 
urfoil Mach Number of 0.8 and a wide range of gust front 
locities. The abscissa is specified in terms of chord lengths 


led, s, which is determined from the time by 
s= Mt = UT/ 1] 


plete set of numerical results is presented in reference 7 for 
irfoil Mach Number range from 0.5 to 0.9. Reference 7 also 


neludes a discussion of methods of approximating the lift re 


sponse functions for large values of time These methods are 


Dased on the premise that the asymptotic behavior of the indicial 


READERS’ 


FORUM 21 





lift response function in compressible flow is the same t 
incompressible flow rhe validity of this premise is disctussec 
reference 2 where it is also shown that the distribution of pressurt 
wer the airfoil chord attains the steady-state distril 
only several chord lengths traveled for the 
front velocity 
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A Remark on the Optimum Design of Round 
Tubing 


William L. Poesch 
Graduate Student, College of Engineering, University of Califorr 
/ iid 


Los Angeles, Cal 


December 2, 1957 


FMIS NOTE presents experimental verification of a previously 
I published theory on the optimum design of round tubular 
columns in axial compression 

This investigation was directed toward the analytical deter 
mination of the optimum diameter /thickness ratio of a pin-ended 
tubular column using an experimentally obtained optimum struc 


tural index his ratio was compared with the actual ratio of 
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Fic. 1 Optimum (maximum) stress for round tubular colum1 
of 2024-T3 aluminum 
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diameter and thickness of the specimens to obtain a measure of 
the accuracy of the theory 

Twenty-one test specimens with lengths varying from 4.0 to 
89.0 inches were utilized in the tests. The material, 2024-T3 
aluminum alloy, and a D/t ratio of 70.5 were used throughout the 
experiment 

For convenience, two relations from reference | are shown b« 


low: 
= 2[k.2E/a(P 
kowE, 
Iq. (1) gives the optimum cross section for a given material and 
structural index. Eq. (2) was used to plot the optimum stress 
versus structural index curve shown in Fig. 1. Note that to plot 
this curve, it is only necessary to have the material properties 
These were obtained from compression tests on the column ma 

terial 

After data on the buckling or crippling load had been obtained 
for each column, it was possible to plot the test data as shown in 
As the length was varied, the test points approached the 
The optimum point 


Fig. | 
optimum curve and then broke away again 
was obtained as (P?/L?)'/? = 6.9 
Substituting the experimentally obtained optimum structural 
index into Eq. (1), an optimum D/t ratio of 74.2 was obtained 
An error of 5.25 per cent resulted when this value was compared 
with the actual D/f ratio of the specimens. A value of 0.84 was 
on the basis of experimental data 
California at Los 


used for kz in Eqs. (1) and (2), 


previously obtained at the University of 


Angeles, Calif.* 
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On Optimum Nose Shapes for Missiles in the 
Superaerodynamic Region 


David H. Dennis 

Aeronautical Research Scientist, Ames Aeronautical Laboratory, 
NACA, Moffett Field, Calif 

November 97, 1957 


W J. CARTER! and i. D. Chang? have calculated the shape 
© of the body of revolution of given length and base diam 
eter, having minimum drag in hypervelocity free-molecule flow, 
With this 


assumption, the expression defining surface pressure coefficient 


assuming specular reflection from a smooth surface. 


is the same as that for Newtonian flow except for a constant 
factor of 2. 
the same shape, including a small flat nose, as that derived by 


Accordingly, the minimum drag body should have 


Newton.® The parametric equations defining this shape are 


In y 


where y; is the nose radius.‘ 

The flat nose results from the restriction on body length. — If 
the length is not fixed, the minimum drag shapes have sharp 
noses. For example, the cone is the minimum drag shape of given 
surface area and base diameter .4 

As a final point it may be noted that since the local drag on a 
body in the flows considered here depends only on local slopes 
and body radii, the first integral to the Euler equation may be 
written immediately and the body coordinates in terms of body 


slopes are easily deduced 
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Structural and Impact Loads 


for the Flexible Airplane. ., 


Continued from page 170) 


airplane was obtained from experimental measurements 
over a range of impact conditions. From these meas- 
urements the hydrodynamic force and moment on the 
airplane were obtained. The most critical parameter 
in measuring the effects of landing impact was found 
to be the rate of vertical descent. 

A comparison of the experimental impact data with 
the results of the analysis showed that the analytical re- 
sults generally were in good agreement with the experi- 
This includes both time 
variation and magnitude. It that the 
hydrodynamic step impact and the dynamic response 
of the flexible seaplane may be accurately given by the 


agreement 
is concluded 


mental results. 


method of analysis presented in this paper. 

The response of the flexible airplane was found to be 
a function of both the magnitude and shape of the fore- 
ing function. The maximum inertial loading for some 
structures was not obtained for the conditions of maxi- 


mum hydrodynamic loading. In order to obtain the 


maximum loads for the flexible airplane it was neces- 
sary to study a range of landing impacts. 

The inclusion of airplane flexibility in the equations 
of motion tended to reduce the hydrodynamic foree 
obtained for the rigid airplane. In this instance the 
maximum reduction was of the order of 7 per cent. 
However, for a more flexible airplane this reduction 


may be of a large magnitude. 
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